Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,328,182 A→T 6.2% intergenic (+99/+57) 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,1820AT 6.2% 441.3 / 14.6 129intergenic (+99/+57)yhbY/greARNA binding protein associated with pre‑50S ribosomal subunits/transcript cleavage factor
Reads supporting (aligned to +/- strand):  ref base A (56/64);  new base T (6/2);  total (62/67)
Fisher's exact test for biased strand distribution p-value = 1.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 131 total aligned reads displayed.

TGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAG  >  NC_000913/3328034‑3328328
                                                                                                                                                    |                                                                                                                                                   
tGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCa                                                                                                                                                     >  4:759697/1‑149 (MQ=255)
tGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCa                                                                                                                                                     >  4:504732/1‑149 (MQ=255)
   tGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑g                                                                                                                                                 >  4:10632/1‑149 (MQ=255)
            gCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCttt                                                                                                                                         <  1:408604/148‑1 (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)
                   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)
                     aaGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCctttctt                                                                                                                               <  3:738469/149‑1 (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                                                                                                                     <  3:916171/149‑1 (MQ=255)
                               aaTCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACAt                                                                                                                     >  4:784096/1‑149 (MQ=255)
                                  ctcGCTGCGACGCTAAGATTATCCTAGAGTTACACTCGTTTGTTGTGTAAATCGAGCGGTTTTCCGCAGGCAGGAGGGCAAAACGCTACGCTCTGGGCGTTAATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACGATACATCaa                                                                                                                  <  4:162690/149‑1 (MQ=255)
                                  ctcGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCaa                                                                                                                  <  4:55724/149‑1 (MQ=255)
                                   tcGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAAc                                                                                                                 <  3:200992/149‑1 (MQ=255)
                                   tcGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAAc                                                                                                                 >  4:933164/1‑149 (MQ=255)
                                       tGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATc                                                                                                              <  1:12612/148‑1 (MQ=255)
                                        gCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCtt                                                                                                            <  3:442714/149‑1 (MQ=255)
                                         ccACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCtt                                                                                                            <  1:284493/148‑1 (MQ=255)
                                          cACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTg                                                                                                           <  1:483253/148‑1 (MQ=255)
                                          cACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTg                                                                                                           <  2:653591/148‑1 (MQ=255)
                                           aCGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAg                                                                                                         <  3:66375/149‑1 (MQ=255)
                                             gCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTa                                                                                                       >  4:456407/1‑149 (MQ=255)
                                              cTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTa                                                                                                       >  1:121061/1‑148 (MQ=255)
                                              cTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTa                                                                                                       >  2:610833/1‑148 (MQ=255)
                                               tAAGATTATCCTAAAGTTACACACTTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTAtt                                                                                                     <  3:116131/149‑1 (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                                                                                                    >  2:590279/1‑148 (MQ=255)
                                                 aGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTg                                                                                                    <  1:889848/148‑1 (MQ=255)
                                                    ttATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGt                                                                                                 >  1:50982/1‑148 (MQ=255)
                                                      aTCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAAt                                                                                              <  4:67857/149‑1 (MQ=255)
                                                         cTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTc                                                                                            >  1:427521/1‑148 (MQ=255)
                                                          tAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCt                                                                                           >  2:785529/1‑148 (MQ=255)
                                                          tAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCt                                                                                           >  2:887819/1‑148 (MQ=255)
                                                          tAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCGGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCtt                                                                                          >  3:573135/1‑149 (MQ=255)
                                                           aaaGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCtt                                                                                          <  2:427521/148‑1 (MQ=255)
                                                           aaaGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCtt                                                                                          >  1:203753/1‑148 (MQ=255)
                                                            aaGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTAc                                                                                        <  4:291089/149‑1 (MQ=255)
                                                                   acacaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTAtt                                                                                 >  3:930253/1‑149 (MQ=255)
                                                                   acacaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGC‑TATGCGGCCtttta                                                                                                                                       <  4:556325/95‑2 (MQ=255)
                                                                   acacaTTTGCTCTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCTAT‑GCGGCCtttta                                                                                                                                       >  3:556325/1‑94 (MQ=255)
                                                                    cacaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTAtt                                                                                 >  2:948250/1‑148 (MQ=255)
                                                                     acaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACaa                                                                                                                          >  4:304843/1‑106 (MQ=255)
                                                                     acaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACaa                                                                                                                          <  3:304843/106‑1 (MQ=255)
                                                                     acaTTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTcc                                                                               <  3:604457/149‑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)
                                                                                         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)
                                                                                                          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)
                                                                                                                   aaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtgatga                                 >  4:314282/1‑149 (MQ=255)
                                                                                                                   aaTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtgatg                                  >  1:69508/1‑148 (MQ=255)
                                                                                                                    aTGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGAc                                >  3:836486/1‑149 (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                     <  1:91079/148‑1 (MQ=255)
                                                                                                                                gTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAAcataat                     <  1:247854/148‑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:10632/149‑1 (MQ=255)
                                                                                                                                          aaaaGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTTAGTATTGGGTAATTCTTACAGGTATTTCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCttt           >  1:255192/1‑148 (MQ=255)
                                                                                                                                          aaaaGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCgccgcc                                               >  2:137889/1‑112 (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                                     <  4:881580/119‑1 (MQ=255)
                                                                                                                                             aGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTtg                                     >  3:881580/1‑119 (MQ=255)
                                                                                                                                                   cATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAg  >  1:630372/1‑148 (MQ=255)
                                                                                                                                                   cATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAg  >  2:347366/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
TGGTGCTTTATCGCCCAACTAAAGAACGTAAAATCTCGCTGCCACGCTAAGATTATCCTAAAGTTACACACATTTGCTGTGTAAAACGAGGGGTTTTCCGCAGGCAGGAGAGCAAAATGCCACGCTCTGTTCGTTGATAAAAGGCCGCATA‑GCGGCCTTTTTCCTTTCTTTACAATACATCAACATCTTGAGTATTGGGTAATTCTTACAGGTATTCCACCTTAATTACTTCAAATTCTACTTCGCCGCCCGGCGTTTTGATGACCACAACATCATCTTCTTCTTTGCCGATCAG  >  NC_000913/3328034‑3328328

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