Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,328,184 A→T 6.2% intergenic (+101/+55) yhbY → / ← greA RNA binding protein associated with pre‑50S ribosomal subunits/transcript cleavage factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,328,1840AT 6.2% 452.1 / 15.9 130intergenic (+101/+55)yhbY/greARNA binding protein associated with pre‑50S ribosomal subunits/transcript cleavage factor
Reads supporting (aligned to +/- strand):  ref base A (56/66);  new base T (6/2);  total (62/68)
Fisher's exact test for biased strand distribution p-value = 1.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

Only 100 of 130 total aligned reads displayed.

TGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGC  >  NC_000913/3328037‑3328330
                                                                                                                                                   |                                                                                                                                                   
tGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑g                                                                                                                                                   >  3:527277/1‑149 (MQ=255)
tGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑g                                                                                                                                                   >  4:10632/1‑149 (MQ=255)
         gCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCttt                                                                                                                                           <  1:408604/148‑1 (MQ=255)
           ccAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCttttt                                                                                                                                         <  2:175785/148‑1 (MQ=255)
           ccAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTc                                                                                                                                        >  4:379125/1‑149 (MQ=255)
           ccAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTc                                                                                                                                        <  3:759697/149‑1 (MQ=255)
              aCTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCct                                                                                                                                      >  1:743/1‑148 (MQ=255)
               cTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCcttt                                                                                                                                    <  4:523550/149‑1 (MQ=255)
                tAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCTAT‑GCGGCCtt                                                                                                                                            >  3:941484/1‑140 (MQ=255)
                tAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCctttc                                                                                                                                   <  3:691940/149‑1 (MQ=255)
                tAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCcttt                                                                                                                                    <  2:510050/148‑1 (MQ=255)
                tAAAGAACCAAAAATCCCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTGCcttt                                                                                                                                    <  2:471399/148‑1 (MQ=255)
                 aaaGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCctttc                                                                                                                                   >  2:425523/1‑148 (MQ=255)
                    gAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTa                                                                                                                               >  3:729343/1‑149 (MQ=255)
                     aaCGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCctttct                                                                                                                                  >  1:516446/1‑145 (MQ=255)
                     aaCGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCctttct                                                                                                                                  <  2:516446/145‑1 (MQ=255)
                        gTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑Gc                                                                                                                                                  >  2:607406/1‑126 (MQ=255)
                         tAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAAt                                                                                                                           >  2:767337/1‑148 (MQ=255)
                         tAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATa                                                                                                                          <  3:12975/149‑1 (MQ=255)
                            aaTCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACAt                                                                                                                       >  4:784096/1‑149 (MQ=255)
                               ctcGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCaa                                                                                                                    <  4:55724/149‑1 (MQ=255)
                                tcGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAAc                                                                                                                   >  4:933164/1‑149 (MQ=255)
                                    tGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATc                                                                                                                <  1:12612/148‑1 (MQ=255)
                                      ccACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCtt                                                                                                              <  1:284493/148‑1 (MQ=255)
                                       cACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTg                                                                                                             <  2:653591/148‑1 (MQ=255)
                                        aCGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAg                                                                                                           <  3:66375/149‑1 (MQ=255)
                                           cTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTa                                                                                                         >  1:121061/1‑148 (MQ=255)
                                           cTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTa                                                                                                         >  2:610833/1‑148 (MQ=255)
                                            tAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTAtt                                                                                                       >  4:725388/1‑149 (MQ=255)
                                             aaGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTAtt                                                                                                       >  1:321639/1‑148 (MQ=255)
                                             aaGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTAtt                                                                                                       >  1:722505/1‑148 (MQ=255)
                                             aaGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTg                                                                                                      <  3:886333/149‑1 (MQ=255)
                                              aGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTg                                                                                                      <  2:179780/148‑1 (MQ=255)
                                              aGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTg                                                                                                      <  1:610833/148‑1 (MQ=255)
                                              aGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTg                                                                                                      <  1:889848/148‑1 (MQ=255)
                                                 ttATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGt                                                                                                   >  1:50982/1‑148 (MQ=255)
                                                      cTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTc                                                                                              >  1:427521/1‑148 (MQ=255)
                                                       tAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCt                                                                                             >  2:887819/1‑148 (MQ=255)
                                                       tAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCt                                                                                             >  2:785529/1‑148 (MQ=255)
                                                       tAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCGGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCtt                                                                                            >  3:573135/1‑149 (MQ=255)
                                                        aaaGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCtt                                                                                            >  1:203753/1‑148 (MQ=255)
                                                        aaaGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCtt                                                                                            <  2:427521/148‑1 (MQ=255)
                                                        aaaGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTa                                                                                           <  4:527277/149‑1 (MQ=255)
                                                        aaaGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTa                                                                                           <  3:725388/149‑1 (MQ=255)
                                                         aaGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTAc                                                                                          <  4:291089/149‑1 (MQ=255)
                                                         aaGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTAc                                                                                          <  4:177346/149‑1 (MQ=255)
                                                                acacaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTAtt                                                                                   >  3:930253/1‑149 (MQ=255)
                                                                acacaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGC‑TATGCGGCCtttta                                                                                                                                         <  4:556325/95‑2 (MQ=255)
                                                                 cacaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTAtt                                                                                   >  2:948250/1‑148 (MQ=255)
                                                                  acaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACaa                                                                                                                            <  3:304843/106‑1 (MQ=255)
                                                                      ttGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCt                                                                             <  4:729343/149‑1 (MQ=255)
                                                                          tgtgTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAAt                                                                         <  3:784096/149‑1 (MQ=255)
                                                                                aaaCGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACtt                                                                    <  2:300032/148‑1 (MQ=255)
                                                                                 aaCGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTc                                                                   >  2:247854/1‑148 (MQ=255)
                                                                                   cGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCaaa                                                                <  3:143379/149‑1 (MQ=255)
                                                                                      ggggTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAAtt                                                              >  1:433905/1‑148 (MQ=255)
                                                                                      ggggTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTc                                                             >  4:336961/1‑149 (MQ=255)
                                                                                       gggTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCt                                                            <  3:350408/149‑1 (MQ=255)
                                                                                       gggTCTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCGGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACAGCTTGAGTATTGGCTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCt                                                            >  4:16868/1‑149 (MQ=255)
                                                                                        ggTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTa                                                           <  3:504732/149‑1 (MQ=255)
                                                                                        ggTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCAAA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTa                                                           <  3:933164/149‑1 (MQ=255)
                                                                                          ttttCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACt                                                         >  3:426474/1‑149 (MQ=255)
                                                                                          ttttCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACt                                                         >  3:449030/1‑149 (MQ=255)
                                                                                            ttCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACtt                                                        <  1:873637/148‑1 (MQ=255)
                                                                                                gcaggcagGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCt                                                                                                               >  1:473900/1‑89 (MQ=255)
                                                                                                    gcaggAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCg                                               >  3:110608/1‑149 (MQ=255)
                                                                                                       ggAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCg                                            <  3:290227/149‑1 (MQ=255)
                                                                                                         agagCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCAGGCGt                                           >  1:820816/1‑148 (MQ=255)
                                                                                                          gagCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGtt                                          <  2:722505/148‑1 (MQ=255)
                                                                                                           agCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGttt                                         <  2:328704/148‑1 (MQ=255)
                                                                                                             ccaaaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtga                                      <  3:571768/147‑1 (MQ=255)
                                                                                                             cAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtga                                      <  4:873751/149‑1 (MQ=255)
                                                                                                              aaaaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtga                                      <  1:531376/148‑1 (MQ=255)
                                                                                                               aaaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtgatg                                    >  3:437411/1‑149 (MQ=255)
                                                                                                                aaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtgatga                                   >  4:314282/1‑149 (MQ=255)
                                                                                                                aaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtgatg                                    >  1:69508/1‑148 (MQ=255)
                                                                                                                   gCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACca                                <  3:314282/149‑1 (MQ=255)
                                                                                                                      aCGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACcaca                              >  1:768458/1‑148 (MQ=255)
                                                                                                                          tctGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCAcatcat                          <  1:357611/148‑1 (MQ=255)
                                                                                                                             gTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAAcatcat                       >  2:52886/1‑148 (MQ=255)
                                                                                                                             gTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAAcatcat                       <  1:91079/148‑1 (MQ=255)
                                                                                                                             gTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAAcataat                       <  1:247854/148‑1 (MQ=255)
                                                                                                                             gTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACatcatc                      <  3:898056/149‑1 (MQ=255)
                                                                                                                                cGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttc                   <  3:379125/149‑1 (MQ=255)
                                                                                                                                   tGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCAtcttctt                 >  1:795053/1‑148 (MQ=255)
                                                                                                                                    ggaaaaaGGCCGCTAT‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACtt                                                                    >  1:370027/4‑96 (MQ=39)
                                                                                                                                    ggaaaaaGGCCGCTAT‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACtt                                                                    <  2:370027/93‑1 (MQ=38)
                                                                                                                                      tAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCttt             >  3:814224/1‑149 (MQ=255)
                                                                                                                                      tAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCttt             >  3:581179/1‑149 (MQ=255)
                                                                                                                                      tAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCttt             <  3:10632/149‑1 (MQ=255)
                                                                                                                                       aaaaGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTTAGTATTGGGTAATTCTTACAGGTATTTCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCttt             >  1:255192/1‑148 (MQ=255)
                                                                                                                                       aaaaGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCgccgcc                                                 <  1:137889/112‑1 (MQ=255)
                                                                                                                                        aaaGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTg            >  2:494518/1‑148 (MQ=255)
                                                                                                                                         aaGGCCGCTAT‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGcc          >  3:121395/1‑149 (MQ=255)
                                                                                                                                         aaGGCCGCTAT‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGc           >  1:629686/1‑148 (MQ=255)
                                                                                                                                          aGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtg                                       >  3:881580/1‑119 (MQ=255)
                                                                                                                                          aGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtg                                       <  4:881580/119‑1 (MQ=255)
                                                                                                                                                cATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAg    >  1:630372/1‑148 (MQ=255)
                                                                                                                                                cATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAg    >  2:347366/1‑148 (MQ=255)
                                                                                                                                                  tA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGc  <  1:753639/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAGGC  >  NC_000913/3328037‑3328330

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑