Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,991,593 A→C 13.7% S9R (AGT→AGG yecA ← UPF0149 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,991,5930AC13.7% 185.2 / 2.6 73S9R (AGT→AGGyecAUPF0149 family protein
Reads supporting (aligned to +/- strand):  ref base A (43/18);  new base C (3/9);  total (46/27)
Fisher's exact test for biased strand distribution p-value = 6.49e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.79e-01

TCTGCGGTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGT  >  NC_000913/1991357‑1991818
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   
tctGCGGTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGAGGTCAGTGGTCTATTTGGTCCGAACATCGTCCAGCCTTTCCACCCCACt                                                                                                                                                                                                                                   >  1:446518/1‑239 (MQ=255)
    cGGTATCGGCCATATGTTGAAACGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGCATGCCGTGGTCAACGTTGTATTTGGTCAGAATATCGTCCAGCCCTTACAACTCACTTTCg                                                                                                                                                                                                                               >  1:463797/1‑239 (MQ=255)
         tCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCAATTCCAAACAACTTTCGTTTaa                                                                                                                                                                                                                          >  1:673885/1‑239 (MQ=255)
                       aaaaGCCCAGTTTCATAAAGCCGGTCACTCTTTTCCAGGCCCCCCGGGCCGGCACATAGCAGCCCCACCCCCACACGCCACCCAGCACCTGTCCCGGTTCAATCTCTTGCGGAAACCCCACCACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTc                                                                                                                                                                                                             <  1:264226/239‑1 (MQ=255)
                       aaaaGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAAccttcctcttatac                                                                                                                                                                                                            >  1:87862/1‑226 (MQ=255)
                                 aTTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCttt                                                                                                                                                                                                    >  1:129036/1‑238 (MQ=255)
                                      aaacgccccttgcccttctctcTGACCCCCCGCGCGTCACATAGCCGCCCCCCCCCCCCCCGCAACCCGCCCTCGGTCCCGTTACAACGCCGCGGGGAGAACTCAACACGGCCGTCAATAACCCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTtg                                                                                                                                                                                              <  2:129036/221‑1 (MQ=255)
                                             cgTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGCAGAACTCAACACCGCCGTCAATAACCCGTCCAGCTCCGCCACATCAAGGATGGCGTGGCCAGGGTTGTCTTTGTTCAGAATATCGTCCAGCCATTCCAACTCCATTTCGTTTTACGGTCCCGTTTTCATCAGATTTTCCTtgtgtcacc                                                                                                                                                                                       >  2:785089/1‑234 (MQ=255)
                                                        tctcTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGGCAGAATATCGTTCAGCCCTTCCCACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGGGGATCTCAACTCTCCac                                                                                                                                                                            >  2:126758/1‑238 (MQ=255)
                                                            tGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGGAACCAGCAACTGTTCCGGTTCAATCTCGTGCGGAGAACTCAACACCGCCGTCAATAAACAGTCCAGCTCAGCAACACCAAGGAAGACGTGGTCAGTTTTGCAATTGTTCAGAATACCGTCAATCCATTCAACCGCACTTTCGTTTAACGGTCCCGCTTCATTAACCTTTTCATTATGGATCTaaaatgaccacaccc                                                                                                                                                                        >  2:741643/1‑227 (MQ=255)
                                                            tGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCCTTCCCACTCCCTTTCGTTTACCGGTGCCGTTTTCATACGCATTTCCTTGTGGAACTCAACCCCCCAGCAcc                                                                                                                                                                        >  2:357009/1‑239 (MQ=255)
                                                                   cccgcgcagcccatctccccccccccccccccCGGCACCCGCCCCTTGTCCGGCTTCAATCCCGTGGGGAAAAACCACCCCCGCCTGAAATAACCCGCCCGCCCCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCtt                                                                                                                                                                <  2:446518/219‑1 (MQ=255)
                                                                      cgcgGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTCCAt                                                                                                                                                            >  1:147176/1‑239 (MQ=255)
                                                                                 aGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGt                                                                                                                                                  >  1:272226/1‑239 (MQ=255)
                                                                                      gCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCCCTTTCGTTTAACGGTCCCGGTTTCATACCCTTTTCCTTATGGATCCCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCag                                                                                                                                             >  1:666686/1‑239 (MQ=255)
                                                                                      gCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCag                                                                                                                                             >  1:222751/1‑239 (MQ=255)
                                                                                      gCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGGGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCGTTTCCTTGTGGATCTCAATTCCCCAGCACCTA‑‑TCTTACATGCCAGTCCGTATCac                                                                                                                                             >  2:587059/1‑238 (MQ=255)
                                                                                          cccccccacCCGGCAACCAGCCACTTTTCGGGTTCAATCGTGTGGGGAAAACTAAACCCCGCCGTAAATAAACCGGCCGGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCCCTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAgaga                                                                                                                                          <  1:412022/236‑1 (MQ=255)
                                                                                           ccccccacACGGCCACCCCGCCCTTTTTCCGTTTAAATTCTTGGGGGAAACCTAAACCCCGCCGTCAATAACCCGTCCAGCCCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAgaga                                                                                                                                          <  1:73832/237‑1 (MQ=255)
                                                                                                  cacGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCGGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACAAGCCGGTCCGAATCAAAGAAACttttt                                                                                                                                  >  2:351423/1‑239 (MQ=255)
                                                                                                  cacGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGCGATACttttt                                                                                                                                  >  1:537838/1‑239 (MQ=255)
                                                                                                  cacGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCCAAACTCCCCAGCACCTA‑‑TCTTACATGCCGGGCCCTACAAAAGAAACttttg                                                                                                                                  >  2:412022/1‑238 (MQ=255)
                                                                                                      gCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTT‑‑TCTTACATCCCGGTCCGTATCAGACATCCTTTTTGatt                                                                                                                              >  1:843705/1‑237 (MQ=255)
                                                                                                                        ccGTTTCAATCTCTGGCGGAAAACTCAACGCCGCCGTCAATAAACCTTCCAGCCCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCTAGCCATTCCAACTCCCTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGAtt                                                                                                              <  1:250688/239‑1 (MQ=255)
                                                                                                                         ggttttCAACCTCTGGCAGAAACCCCAACCCCGCCGTCAATAACCCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGAtt                                                                                                              <  1:537572/235‑1 (MQ=255)
                                                                                                                          cggttcaactcttggggggtgaACTCACCACCGCGTGAAATAACCGGCCGACCCCCCCCACATCAAGGATGGCGTGGTCAGTTTTTTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCCCTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGAtt                                                                                                              <  2:105150/219‑1 (MQ=255)
                                                                                                                           gTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGACTGGCTTTCATGGCGATTaaa                                                                                                           >  1:105150/1‑239 (MQ=255)
                                                                                                                                      ttgggggaaactccacgcccgcggTAATTAATCCGTCCGGCCCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAagga                                                                                                  <  2:222751/221‑1 (MQ=255)
                                                                                                                                      ctggggagaagctaaccccgtccgTAATTAACCCGGCCAGCTCCGCCACAGCAAGGATGGCGTGGTCAGTGTTGTATTTGGTAAGAATATGGTCCAGCCATTCCAACTCCCTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAagga                                                                                                  <  1:587059/218‑1 (MQ=255)
                                                                                                                                      ctggggaagaaccaaccacCCCCGTCAATAACCCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAagga                                                                                                  <  1:81801/223‑1 (MQ=255)
                                                                                                                                       tGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTAAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGAACTCAACTAGCCAGACCCTA‑‑TCTTCCATGCCGGTCGGTATCAGCGATACTTTTTGAGTGGATTTGTTGGTGTTTAACCATAAaggagg                                                                                                >  1:829880/1‑239 (MQ=255)
                                                                                                                                                 cTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCtt                                                                                                                                                                >  1:202119/1‑166 (MQ=255)
                                                                                                                                                           ccgTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATTCCGGCCCGTATCAAGGATACTTTTTGAGTGGCTTTGCTGGTGTTTAAAAATTATGGGGGGGGTAAAGCAAAATTGCAg                                                                             >  2:792654/1‑239 (MQ=255)
                                                                                                                                                           ccgTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGGGTGGCTTTGCTGGGGATTAAAAATTAAGAAGGGGGGAAAGAAAAATTGCAg                                                                             >  2:59468/1‑239 (MQ=255)
                                                                                                                                                           ccgTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGAAAAGTTGCAg                                                                             >  1:629337/1‑239 (MQ=255)
                                                                                                                                                              tCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAACATTAAGGGGGGGTGAACGACAAGTTGCAGGca                                                                          >  1:795496/1‑239 (MQ=255)
                                                                                                                                                                        gTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGAc                                                                                      >  2:582253/1‑216 (MQ=255)
                                                                                                                                                                        gTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGAc                                                                                      <  1:582253/216‑1 (MQ=255)
                                                                                                                                                                         tCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGGGTAACGCCAAGTTGCAGGCACAAAAAAacc                                                                >  1:669075/1‑238 (MQ=255)
                                                                                                                                                                         tCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGGGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAAacc                                                                >  1:719974/1‑238 (MQ=255)
                                                                                                                                                                                   gccaactaaagtgtGGGGTGGTCAGTGTTGTTTTTGGTCAAAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTAATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCTGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAAACCACCCGAAggt                                                      <  2:147176/226‑1 (MQ=255)
                                                                                                                                                                                   aagtttttttaaggtgaatGGTTAGTTTTTTTTTTTGGCCGAATTTCCTCCGTCCATTCCCACCCCCTTTCGTTTACCCGTCCCGTTTTCATCCCCTTTTCCTGTTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAAACCACCCGAAggt                                                      <  2:51690/221‑1 (MQ=255)
                                                                                                                                                                                    ccACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTTCCTGGGATTAAAAATTAAGGAGGGGGTCACGCCAAGTTGCAGGCCCAAAAAAACCAcccacagggg                                                     >  2:630638/1‑232 (MQ=255)
                                                                                                                                                                                    ccACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGTTTAAAAATTAAGGAGGGTGTAACg                                                                                        >  2:582089/1‑204 (MQ=255)
                                                                                                                                                                                    ccACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACg                                                                                        <  1:582089/204‑1 (MQ=255)
                                                                                                                                                                                          cAAGGTTGGCGTTGTTGTTTTGGTTTTGGGTCAGAATTTCGTCCAGCCATCCCAACTCCCTTTCGTTTACCGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAAGTGGCCAGCGCCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAAACCACCCGAAGGTGGTTTCa                                               <  2:843705/239‑1 (MQ=255)
                                                                                                                                                                                          cAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAAACCACCCGAAGGTGGTTTCa                                               <  1:351423/239‑1 (MQ=255)
                                                                                                                                                                                            aGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGCTTAAAAATTAAGGAGGGTGTACAGCCAAGTTGCAGGCACAAAAAAACCACCCAAAGGGGGTTTAAca                                             >  2:239248/1‑238 (MQ=255)
                                                                                                                                                                                                tgtgttGGTCTGTGTTTTTTTTGTGCAGAATATCGGCAAGCCTGTCCAACCTCTTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATTTCACCTCGCCAGCCCCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAAACCACCCGAAGGTGGTTTCACGaca                                          <  2:795496/233‑1 (MQ=255)
                                                                                                                                                                                                tggtgtgttccttgttttTTTTGGTCAAATTTTCGGCCAGCCATTCCAATCCCATTTCGTTTAACGGTCCCGTTTTCATACGCTTTGCCTTGTGGAGCTCAACTCGCCAGCCCCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGAGGCTTTGCTGGTGAATAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGTCACAAAAAAACCACCCGAAGGTGGTTTCACGaca                                          <  2:402533/226‑1 (MQ=255)
                                                                                                                                                                                                tGGCGTGGTCAGTTTTTTTTTTGGTCAAAATACCGTCCAGCCATTCCATCTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTGGCCTGGGGGAACTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAAACCACCCGAAGGTGGTTTCACGaca                                          <  2:272226/238‑1 (MQ=255)
                                                                                                                                                                                                tGGCGTGGTCAGTGTTGTATTTGGTCAGAATTGCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTGTCCATGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAAACCACCCGAAGGTGGTTTCACGaca                                          <  1:126758/238‑1 (MQ=255)
                                                                                                                                                                                                   cGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGCTTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAAACCACCCGAAGGTGGTTTCACGACACt                                        >  1:83006/1‑237 (MQ=255)
                                                                                                                                                                                                      ggTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTATCTCTTACATGCCGGTCCGTATCAgaga                                                                                                                                          <  1:483434/138‑1 (MQ=255)
                                                                                                                                                                                                      ggTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTATCTCTTACATGCCGGTCCGTATCAgaga                                                                                                                                          >  2:483434/1‑138 (MQ=255)
                                                                                                                                                                                                        tCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAaggag                                                                                                 >  1:616145/1‑175 (MQ=255)
                                                                                                                                                                                                        tCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAaggag                                                                                                 <  2:616145/175‑1 (MQ=255)
                                                                                                                                                                                                            tgttgtATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGGGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGCCAAGTTGCCGGCCCAAAAAAACCACCCGAAggggggtcaacgatatgcattcttgttt                             >  2:888480/1‑216 (MQ=255)
                                                                                                                                                                                                             tttttttttGGTCAGAATTTCGGCCAGCTATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACCGTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCtt                             <  2:669075/232‑1 (MQ=255)
                                                                                                                                                                                                                 tATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAAACCACCCGAAGAGGGTTTCACGa                                            >  1:51690/1‑219 (MQ=255)
                                                                                                                                                                                                                      ggTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTtg                                                                                                                                                                                              >  2:275487/1‑70 (MQ=255)
                                                                                                                                                                                                                      ggTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTtg                                                                                                                                                                                              <  1:275487/70‑1 (MQ=255)
                                                                                                                                                                                                                      ggTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTC‑‑TCTTCCATGCCGGTCCGTATCAGCGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGCAAAGTTGCAGGCACAAAAAAACCACCCGAAGGTGGTTTCACGAAATTGCTTCTTGCTTTGATTTTAtt                   >  1:436889/1‑239 (MQ=255)
                                                                                                                                                                                                                      ggTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCaga                                                                                                                                            >  2:554163/1‑118 (MQ=255)
                                                                                                                                                                                                                      ggTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCaga                                                                                                                                            <  1:554163/118‑1 (MQ=255)
                                                                                                                                                                                                                      ggTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGGGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGCCAAGTTGCAGGCACAAAAAAACCACCCCAAGGGGGTTTCACGAACATGCTTATTGCTTTACTTTTAtt                   >  2:461882/1‑239 (MQ=255)
                                                                                                                                                                                                                        tCAGAATATCGTCCAGCCATTCCAACTCCCTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGAtt                                                                                                              <  2:801538/146‑1 (MQ=255)
                                                                                                                                                                                                                        tCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGAtt                                                                                                              >  1:801538/1‑146 (MQ=255)
                                                                                                                                                                                                                           gAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACCAAAAAACCACCCGAAGGTGGTTTCACGaca                                          >  1:174690/1‑211 (MQ=255)
                                                                                                                                                                                                                                  gTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTTTTGCTTTGATTTTATTCTTATTTTTccc       >  2:665157/1‑239 (MQ=255)
                                                                                                                                                                                                                                      aGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAGAAATTAAGTAGGGCGTAACGACACGTTGCATGCACAACAAAACCACCCCAAGGTGGTTTCACGACACTGCTTATTGCTTTGGTTTTATCTTTATCTTTCCCTTgg   >  1:572817/1‑239 (MQ=255)
                                                                                                                                                                                                                                      aGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTCTTTGAGTGGCTTTGCTCGTGATTAGAACTTAAGGAGGGTGTAACGACAAGTTGCAGGCACAACAAAACCACCCGAAGGTGGTTTCACGACACCGCTTATGGCTTTGATTTTATTCTTATCTTTCCCATgg   >  1:460632/1‑239 (MQ=255)
                                                                                                                                                                                                                                        ccATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGGTTTTATTCTTATCTTTCCCATGGt  >  1:122846/1‑238 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                   
TCTGCGGTATCGGCCATATGTTGAAAAGCCAGATTCATAAAGCGCGTCATCTCTTTCTCTGACGCCCAGCGCGGCACATAGTCAGCCCCACCCCACACGGCAACCAGCCACTGTTCCGGTTCAATCTCTTGCGGAGAACTCAACACCGCCGTCAATAAACCGTCCAGCTCCGCCACATCAAGGATGGCGTGGTCAGTGTTGTATTTGGTCAGAATATCGTCCAGCCATTCCAACTCACTTTCGTTTAACGGTCCCGTTTTCATACGCTTTTCCTTGTGGATCTCAACTCGCCAGCACCTA‑‑TCTTACATGCCGGTCCGTATCAGAGATACTTTTTGAGTGGCTTTGCTGGTGATTAAAAATTAAGGAGGGTGTAACGACAAGTTGCAGGCACAAAAAAACCACCCGAAGGTGGTTTCACGACACTGCTTATTGCTTTGATTTTATTCTTATCTTTCCCATGGT  >  NC_000913/1991357‑1991818

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: