Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 361,668 C→T 19.6% E65K (GAA→AAA)  lacA ← thiogalactoside acetyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913361,6680CT19.6% 52.0 / 4.6 44E65K (GAA→AAA) lacAthiogalactoside acetyltransferase
Reads supporting (aligned to +/- strand):  ref base C (13/22);  new base T (7/2);  total (20/24)
Fisher's exact test for biased strand distribution p-value = 5.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.29e-01

ATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGG  >  NC_000913/361453‑361887
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                  
aTTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGAGGTTGGAACCGTAAGAGAAATAGAAAGGCGGGTTCTACCAAgcggt                                                                                                                                                                                                                         >  2:427954/1‑224 (MQ=255)
     cAATCGTTAACGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTccc                                                                                                                                                                                                                  <  2:528851/228‑1 (MQ=255)
                   aaaGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGCTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGGTTTTCCCCAGGGTGTTTccccacccgggaaaaa                                                                                                                                                                                                     >  1:93202/1‑216 (MQ=255)
                                 gCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTCCTACCCACACGGTTTACCCCACCGTGGCAAAAAATTTTTTATTCAg                                                                                                                                                                                     >  1:195589/1‑228 (MQ=255)
                                 gCCGTTTTTTCTCAATTCATGGAGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGGGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAACATTGCGGCCTATCTGGATGTTGGAACCGTAAGACAAATAAACATACGGGCCTACACTGGCGATTTTCCCGCCAGTTGGAAACAAATCTTTTAttcg                                                                                                                                                                                     >  1:451238/1‑225 (MQ=255)
                                                 tCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTc                                                                                                                                                                     <  2:359479/227‑1 (MQ=255)
                                                          aCAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGATCTAACCCGGGGTTTTCCCCTACCGTGGGAAACAATTCTTTAATCAGGCTTTTTCTTTTTTCAACTTCTgat                                                                                                                                                           >  2:476/1‑228 (MQ=255)
                                                                 tgcccGTAAAGGAAAGAGGAACGTTGGGTGCAATTAGTACGTTATCACCGATTGTTACCGTGGAGTCATTGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCg                                                                                                                                                    <  1:427954/226‑1 (MQ=255)
                                                                 gtCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTTTACCCAAGCGTTTTTCCCTTCACTGGCAAACAATTTTTTTATACAGGTTTCtttttttttcattttttttggtggc                                                                                                                                                     >  2:50131/1‑210 (MQ=255)
                                                                 gtCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCCGGCGTTTTTCCCTACCGTGGCAAACAATTCTTTAATTAGGGATTCTCTTTTTTTAACTTCTtatgggggtc                                                                                                                                                    >  2:96077/1‑223 (MQ=255)
                                                                   cccTTAACGGAAAGAGTAACGTTTGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGa                                                                                                                                                   <  1:606118/227‑1 (MQ=255)
                                                                   cccGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCCTACCAAGCCGTTTTTCCTACCCTGGGAAACAATTTTTTTTTAAGACGTTTCTTTTTTTTCACCtgtccttattcacca                                                                                                                                                   >  1:378357/1‑214 (MQ=255)
                                                                    ccGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTtc                                                                                                                                                                              <  2:378357/199‑1 (MQ=255)
                                                                       tAACGGAAAGAGTAAAGTTGGGTGCAATCAGTTAGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGAt                                                                                                                                              <  1:198194/228‑1 (MQ=255)
                                                                          cGGAAAGAGTAACGGTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGAt                                                                                                                                              <  2:534052/225‑1 (MQ=255)
                                                                              aaGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGATAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTc                                                                                                                                       <  2:174545/228‑1 (MQ=255)
                                                                                    agggTTGGGTTGAATTAGTACGTTAACCTCGGTTGTTACCGTGGAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCAACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACAt                                                                                                                                <  1:80172/225‑1 (MQ=255)
                                                                                      cGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCCGGCGTTTTTCCCTACCGTGGCAAACATTTTTTTTATCTGACTTTCTCTTTTTTCAACTTCTGATGGGTGTGAGTGATTATAATCTTACAtt                                                                                                                               >  2:19252/1‑227 (MQ=255)
                                                                                                 aTCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGAGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGCTTTTCCCCGGGGGTTTTCCCCACCCTGGGAAACATTTCTTTTAGTAGGCTTTCTCTTTTTTCCACTTCTGATGGGTTGGAGTGATTTAACTCATACATTAAAAGTTTTcct                                                                                                                  >  2:80172/1‑227 (MQ=255)
                                                                                                     gTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAAGCGGTTTTCCCTACCGGGGCAAACATTTTCTTAAACAGGCATTCTCTTTTTTTTACTTCTGATGGAGGGCAATGATTTAACTCATACCTTAAACATTTTcccccaa                                                                                                              >  2:570792/1‑224 (MQ=255)
                                                                                                           tATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTTTATCCAGGCGTTTTCCCCTACCACGGGAAAAAATTCTTTAATCAGGACTTCTCTTTTTTTAACATCCGGTGGCTACGGTTGGTTAAACTCAAAAAATTAACGTTTCCCAccgcatttt                                                                                                        >  1:512835/1‑219 (MQ=255)
                                                                                                           tATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTTCCCTACCGTGGCCAACATTTCTTTAATTAAGCTTTCTCTTTTTTCAACTTCTGATGGTTGCGAGTGATTAAAATCCTACATTTACGTTTTTCCCCCAAGTCtt                                                                                                       >  2:545830/1‑227 (MQ=255)
                                                                                                            aTCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAACAATTGCGGCCTATATGGATGTTGGAACCGGAAGAGAAATAGACAGGCGGTTTTTACCAGAGGTTTTTCCCTACCACGGTGAAGATTTTTTTAATACGGCACTTTCTTTTTTCTACTTTTATTTGATTCAAGGTATTAAAACCATACATTTACACTTTTCCACgtaatgttt                                                                                                       >  2:574704/1‑219 (MQ=255)
                                                                                                                ccGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGGGCTTACCAGACGGTTTTCCCCAACCTGGGAAAAAATTTTTTTATCAAGACTTTTCTTTTTTTACCTTTTCGTGGGAGTGTGTGCTTAAACTCATACATTAACGTTTTTCCCACAAGAttttttttt                                                                                                  >  1:637974/1‑227 (MQ=255)
                                                                                                                               tAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGGTCTTCCCAGGAGTTTTTTCCTACCCTGGGAAACATTTCTTTTATCAAGCTTTTTCTTTTTTTAACTTCTTTTGTAGTCGAGCGATTTAAACCATCAACTTACTATTTTCCACCAAGACtttttttcccttacacattcaac                                                                                  >  1:587427/1‑211 (MQ=255)
                                                                                                                                            tGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTTTAACCAAGCGTTTTTCCCTACCCTGGGAAACATTTTTTTTATTAGGCTTTTCCTTTTTTTAACTTTTCATGGGGGCGACTGTGTTAACTCATACATTTACCTTTTTCCCCTAAGGTCTTTTTTCTGAAAACACTCGCACATAACTGTAAAt                                                                      >  2:302557/1‑227 (MQ=255)
                                                                                                                                            tGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACAGTAAGAGAAATAGACAGGCGGTTTTACCCAGGCGGTTTCCCCTACCCAGGGAAACATTTTTTTAATCAAGGATTTTCTTTTTTCACCTTCTGATGGGTGCGGGTGATTAAAACTATACAATAACAGTTTTCCACTAAGGCTTTTTTTCTGGAAACATTTGTACATATCTGTGAATa                                                                     >  1:91959/1‑228 (MQ=255)
                                                                                                                                                          tttGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGctct                                                        <  2:426883/227‑1 (MQ=255)
                                                                                                                                                                   aaaTTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTTACCAGGCGGTTTTCCCTATCCTGGCAAAAAATTCTTTAATACGGCTTTCTCTTTTTTTTACTTTTGTTGTATGCCAGTGGTTAAACTCATACATTTAAGTTTTTCCCCGAAGTGTTTTTTCTGGTATACCTTCTCACATATCTGTAAATAGCTTTCCTGCTCCTTTTATTTct                                              >  2:651506/1‑227 (MQ=255)
                                                                                                                                                                       tgtgtcttattTGGAGTTTGGAACCCGAAAAGAAATATACAGGCTGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCa                                          <  1:476/221‑1 (MQ=255)
                                                                                                                                                                            tctattggggtgtttgagacGTTAGGGAAAAAGACAGGCGTTTTTACCCTGTCGTTTTCCCCTTCCGTGGCAAACCTTTCTTTAATCAGGCTTTCTCTTTGTTCAACTTCTGATGGATGCGAGTGATTAAACTCATCCATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGc                                     <  1:31183/209‑1 (MQ=255)
                                                                                                                                                                                          tgAACCTTAAGAGAAATAGACAGGCTGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGAt                       <  1:89881/227‑1 (MQ=255)
                                                                                                                                                                                          ggAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGGTTTTCCCCTCCGTGGGAAACATTTCTTTAATCAGGCGTTCTCTTTTTTCAACTTCTGATGGATGCGCGGGAGTAAAGTACTAACTTAACGGTTTTTCACGAAAGTTTTTTTCTGTGAACACTTTCCCCCTAACTGTAAATAAATTGCCTGGTCTTGTACTTTCTGTCAtttgaaagtgtaatacgat                       >  1:178034/1‑211 (MQ=255)
                                                                                                                                                                                              ccGTAAGAGAAATAGACCGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCTCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                   <  2:639334/228‑1 (MQ=255)
                                                                                                                                                                                              ccGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTTGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                   <  2:195589/228‑1 (MQ=255)
                                                                                                                                                                                              ccGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                   <  2:93202/228‑1 (MQ=255)
                                                                                                                                                                                              ccGTAAGAAAAATAGAAAGGCTGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGTTCATTGGCATGTTCAATGCGACCACt                   <  2:637974/228‑1 (MQ=255)
                                                                                                                                                                                               ccgtaagagaaaattgtggggggTTCTACCCAGGCTTTTTCCCCTACCGTGGCAAACATTTCCTTAATCAGGCTTTCTATTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                   <  2:316380/207‑1 (MQ=255)
                                                                                                                                                                                                gTAAGAGAAATAGACCGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACt                   <  2:248454/226‑1 (MQ=255)
                                                                                                                                                                                                 taaggaaaaTAGAACGGCGGTTCTACCCAAGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                <  2:91959/222‑1 (MQ=255)
                                                                                                                                                                                                 tAAGAGAAATAGACCGGCGGTTTGACCCATGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                <  1:387081/228‑1 (MQ=255)
                                                                                                                                                                                                  aaggaaattAGACAGGCGGTTCTACCCCTGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTGCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                <  1:96077/219‑1 (MQ=255)
                                                                                                                                                                                                  aaGGGAAATAGACAGGCGGTTCTTCCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCATTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCg                <  2:512835/227‑1 (MQ=255)
                                                                                                                                                                                                                 gCGGTTCTAGCCAGGCTTTTTCCCCTACCGTGGCAAAAATTTCTTTAAACAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCTAGTGATTAAACTCATTCATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTgg  <  1:50131/226‑1 (MQ=255)
                                                                                                                                                                                                                          |                                                                                                                                                                                                                                  
ATTGCCAATCGTTATCGGAAAAGAGTACATCTCGCCGTTTTTTCTCAATTCATGGTGTACAGGGTGTCCCGTAACGGAAAGAGTAACGTTGGGTGCAATCAGTACGTTATCACCGATTGTTACCGTGTAGTCATCGACAATGGTTAAATTGAAATTTGCATAAAAATTGCGGCCTATATGGATGTTGGAACCGTAAGAGAAATAGACAGGCGGTTCTACCCAGGCGTTTTCCCCTACCGTGGCAAACATTTCTTTAATCAGGCTTTCTCTTTTTTCAACTTCTGATGGATGCGAGTGATTAAACTCATACATTAACGTTTTCCCACGAAGTCTTTTTTCCGGTAAGCCTTCGCACATATCGGTAAATAGCTTGCCTGCTCTTATTCTTTCGGTCATTGGCATGTTCAATGCGATCACTCCGTTATGATATGTTGG  >  NC_000913/361453‑361887

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