Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,781,049 C→T intergenic (+32/‑166) lldD → / → trmL L‑lactate dehydrogenase/tRNA (cytidine/uridine‑2'‑O)‑ribose methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,0490CT97.3% 233.8 / ‑7.2 73intergenic (+32/‑166)lldD/trmLL‑lactate dehydrogenase/tRNA (cytidine/uridine‑2'‑O)‑ribose methyltransferase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base T (43/28);  total (45/28)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCG  >  NC_000913/3780807‑3781281
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        
gctgGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGga                                                                                                                                                                                                                                   <  2:36865/250‑1 (MQ=255)
gctgGGTCGTGCTGTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAggg                                                                                                                                                                                                                                    <  2:385554/249‑1 (MQ=255)
   gggTCGTGCTTTCTTGTATGCGATGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCCCTGGCTCCCCTGGCGGAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAgg                                                                                                                                                                                                                               >  1:101369/1‑251 (MQ=255)
             ttCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGGGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTGGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTc                                                                                                                                                                                                                      <  2:281948/250‑1 (MQ=255)
                           ggCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGAGGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCGCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCtgtatgta                                                                                                                                                                                                         <  1:34860/249‑1 (MQ=255)
                             cAACAGCGGTCCAGGCGGGTGTATCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGGAAAGAGTGGCCTTCGGCACGGGCTCCCATGGCGAAAGGTAATGCGGCATAGTCGTTTGCGCCCCCTCTCCCGAACCCTCTCCTTCATTGAGGGGGGACCGGTCGGCGCAGTATGTAcgcc                                                                                                                                                                                                     <  2:34240/251‑4 (MQ=255)
                             cAACAGCGGGCCAGGCGGGTGGAGCTAACCCTCTAAATCTGAACGGAAAAGAGATGAAAGGGGGGATGACGCTGACTGGCGCGAAATCGATCCGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTTGGGAAAGAGGTGGCTGCGGCCCTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCACCTCACCCTCTCCTTCAGGGGGAGGGGGCCGTTCGGCGCTTGATGTACTcc                                                                                                                                                                                                     >  2:471282/1‑251 (MQ=255)
                                 aGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCa                                                                                                                                                                                                 <  1:432401/251‑1 (MQ=255)
                                  gCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCAc                                                                                                                                                                                                <  1:229628/251‑1 (MQ=255)
                                          ggCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCCAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAAc                                                                                                                                                                                        >  1:35122/1‑251 (MQ=255)
                                              ggcgtAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGTGGATGCGGCATAGACGTTGGCGCCCCCTCACCCTAAACCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGCACTCCCTCACTCTGACACGaa                                                                                                                                                                                     >  2:127305/4‑249 (MQ=255)
                                              ggTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGaca                                                                                                                                                                                    <  1:412931/251‑1 (MQ=255)
                                               gtgtAACTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGACTGCGGCATAGTCGTATGCGACCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTgacacgacac                                                                                                                                                                                   >  1:90268/1‑251 (MQ=255)
                                                 gtAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCg                                                                                                                                                                                 <  1:19775/251‑1 (MQ=255)
                                                       aaCCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGACGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGTTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTACCCCTCTCCTTCAGGGAGAGGGTACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCActc                                                                                                                                                                            <  1:519422/250‑1 (MQ=255)
                                                         ccTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttt                                                                                                                                                                         <  1:181247/251‑1 (MQ=255)
                                                            gCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCtttt                                                                                                                                                                        <  1:506278/249‑1 (MQ=255)
                                                             cTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttttt                                                                                                                                                                     >  1:205379/1‑251 (MQ=255)
                                                               aaaTCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCttttttt                                                                                                                                                                     <  1:419884/249‑1 (MQ=255)
                                                                 aTCTGATCGAAAAAGCGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGgcg                                                                                                                                                                                                                  >  2:250022/1‑202 (MQ=255)
                                                                 aTCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGgcg                                                                                                                                                                                                                  <  1:250022/202‑1 (MQ=255)
                                                                  tCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTCCTCCCTCACTCTGAAACGCCCCCGCACTCTTTTTTTCTccc                                                                                                                                                                >  2:399861/1‑251 (MQ=255)
                                                                   cTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGACAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCt                                                                                                                                                               >  2:332936/1‑251 (MQ=255)
                                                                     gATCGAAAAAGAGATGAAAGTGGCGATGTCGCTGACTGGCGCGAAATCGATCAGCGAAAGTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCGCCTAACCCGCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCGCACTCTGAATCCACAGCGCGCTCGTTCTTTCTCCCTCg                                                                                                                                                             <  2:73130/251‑1 (MQ=255)
                                                                      aTCGAAAAAGAGATGAAAGTGGCGCTGACGCTGACTGGCGCGACATCGATCAGCGAAATTAAGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCCTGGCGAAAGGGAATGCGGCATAGTCGATTGCGCCCCCTCCCCCTAACCCTCTCCTTCAGGGCGAGGGGACCGTTCGGCGCTGTATGTCCTCACTAACTCTGAAACGACACCGCCCTCTTTTTTTCTCCCTCGc                                                                                                                                                            >  2:353979/1‑251 (MQ=255)
                                                                      aTCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtct                                                                                                                                                                   >  1:113406/1‑244 (MQ=255)
                                                                      aTCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtct                                                                                                                                                                   <  2:113406/244‑1 (MQ=255)
                                                                           aaaaTAGATGAAAGGGTCGATGACGCTGGGTGGGGCGAAATCGATCAACGAAATATCGCAAGATTCGCTGGTACAGGGGCTGGGTAAAGAGTTGCCTGCGGAACTGGCTCCCATGGCGAAAGGGAATGCGGGATAGTCGTTTGCGCCCCCTCACCCTAACCCTCACCTTCAGGGAGAGGGGACCGTTCGGCGTTGTATGTGCTCCCTCACTGTGAAACGACTACGCACACTTTATTTCTCCCTCGCCCCTc                                                                                                                                                       <  2:341220/251‑1 (MQ=255)
                                                                           aaaaGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTGTCTCCCGCGCCCCTc                                                                                                                                                       <  2:224565/251‑1 (MQ=255)
                                                                              agagATGAAAGTGGCGATGAAGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTAGTTTGCGCCCCCTCACCCTAAACCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAACCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgg                                                                                                                                                    >  1:255736/1‑251 (MQ=255)
                                                                                agaTGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTGGCGCCTCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgggg                                                                                                                                                  <  1:257055/251‑1 (MQ=255)
                                                                                   tGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCgggg                                                                                                                                                  >  1:10184/1‑248 (MQ=255)
                                                                                    gAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACg                                                                                                                                                                                       <  2:111752/210‑1 (MQ=255)
                                                                                    gAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACg                                                                                                                                                                                       >  1:111752/1‑210 (MQ=255)
                                                                                     aaaGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGTGTAAAGAGTTGCCTGCGGCGCTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGTGGACCGTGCGGCGCTGTATGTACTCCCTCGCTCTGAAGCGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAgg                                                                                                                                             <  1:259361/251‑1 (MQ=255)
                                                                                     aaaGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACGCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAgg                                                                                                                                             <  2:35122/251‑1 (MQ=255)
                                                                                         tGGCGATGACGCTGACTTGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGGGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCACCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCACTCGCCCCTCCGGGGAGAGGGCCg                                                                                                                                         >  2:331216/1‑251 (MQ=255)
                                                                                          ggCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                                                                                                        <  1:221035/251‑1 (MQ=255)
                                                                                          ggCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCgg                                                                                                                                        >  1:103897/1‑251 (MQ=255)
                                                                                          ggCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCGGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACGCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGTGCCgg                                                                                                                                        <  2:389372/251‑1 (MQ=255)
                                                                                                   gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTCTTCTCCCTCGCCCCTCCGGGGAGAGGGc                                                                                                                                           >  2:199879/1‑239 (MQ=255)
                                                                                                   gCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTCTTCTCCCTCGCCCCTCCGGGGAGAGGGc                                                                                                                                           <  1:199879/239‑1 (MQ=255)
                                                                                                   gCTGACGGGCGCGAAATAGAGCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGGGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCTTAGGCGTTTGCGCCCCCTCTCCCTACCCCTCCCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACCCTGAAACGACACCGCGCTCTTTTTTTCTCCCTCGCCCCTACGGGGAGAGGGCCGGGGTGAgggg                                                                                                                               <  2:9073/251‑1 (MQ=255)
                                                                                                     tGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGa                                                                                                                              >  1:221350/1‑250 (MQ=255)
                                                                                                      gACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGAAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCa                                                                                                                                                                                                                                       >  1:220818/1‑144 (MQ=255)
                                                                                                      gACTGGCGCGAAATCGATCAGAGAAATTCCGCAAGATCAGCTGGTGCAGGGGCTGGGTAAAGAGATCCCAGAGACACTGGCTCCCATGGCGAAAGGGAATGCGGCCTAGTCGTTTGAGCCCCCTCACCCTAACCCTCACCTTCAGGGAGAGGGGACCGTTCGGCGCCGTATGTACTCCCTCCCTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAAGGCAGGGGTGAGGGGaa                                                                                                                             >  2:563869/1‑250 (MQ=255)
                                                                                                                        cagcgaagttaagaaagttccgcgggTGCACGGGATGGGCACAGAGTTGCGTGAGGCACTGCATCCGATAGAGCAAGGCAATGCGGCATAGTCGTTTGCGCACCCTAACCCTAACCCTCTCCTCCAGGGACAGGGGACCGTTCGGCGATGTATGTACTCCCTCACTCTGAGACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGACACGGCCGCACTGCTCAca                                                                                                           >  2:243841/25‑249 (MQ=255)
                                                                                                                           cGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTTCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAAtt                                                                                                       >  1:422802/1‑251 (MQ=255)
                                                                                                                                    gCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGCGGGGAAAAGGCCGCACTGCTCACCACTTTTCTGCTAt                                                                                              >  1:313361/1‑251 (MQ=255)
                                                                                                                                        gATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCATGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                                                          >  1:539739/1‑251 (MQ=255)
                                                                                                                                         aTTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGc                                                                                         >  1:24039/1‑251 (MQ=255)
                                                                                                                                          ttCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGACACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAACGGCCGCACTGCTCACCAATTTTCTGCTATCATGcc                                                                                        >  1:233159/1‑251 (MQ=255)
                                                                                                                                             gCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGCCCCCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGccccc                                                                                     >  1:394128/1‑251 (MQ=255)
                                                                                                                                              cTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCg                                                                                    >  2:12629/1‑251 (MQ=255)
                                                                                                                                               tGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGcccc                                                                                                                                                         <  1:4732/181‑1 (MQ=255)
                                                                                                                                               tGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGcccc                                                                                                                                                         >  2:4732/1‑181 (MQ=255)
                                                                                                                                               tGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGc                                                                                   >  1:164986/1‑251 (MQ=255)
                                                                                                                                                  tGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACt                                                                                >  2:321707/1‑251 (MQ=255)
                                                                                                                                                              gTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGt                                                                     >  1:51357/1‑250 (MQ=255)
                                                                                                                                                                          cTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCg                                                        >  1:224662/1‑251 (MQ=255)
                                                                                                                                                                                cACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt                                                  >  2:280484/1‑251 (MQ=255)
                                                                                                                                                                                 aCTGGCTCCCCTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACttt                                                  >  2:106169/1‑250 (MQ=255)
                                                                                                                                                                                         ccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTCAGGGGGCAGTATGCTAACCATCGTACTTTACGAAcc                                           >  2:410330/1‑249 (MQ=255)
                                                                                                                                                                                          cATGGCGAAAGGGAATGCGGCATAGACGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGAGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTCCGAAcc                                           >  2:354980/1‑248 (MQ=255)
                                                                                                                                                                                                   ggggAAGGCGGAATAGTGGTTTGCGCCCCCCCCCCCTACCCCCCTCCTACGGGGGGAGGGGCCCGTCCGGCGCTGTATGCACCCCCCCACCCGGAAACGACACCGCACTTTTTTTTTCCCCCCCGCCCCCCCGGGGAGAGGGCGGGGGTGGGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGcc                               <  2:103897/250‑1 (MQ=255)
                                                                                                                                                                                                    gggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGACCCAGAAATTCCGcc                               >  1:505600/1‑250 (MQ=255)
                                                                                                                                                                                                     ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGCGGGGAAAAGGCCGCACTGCTCACCCATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGCACCAGAAATTCCGCCa                              >  1:113408/1‑250 (MQ=255)
                                                                                                                                                                                                       aaTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGcaacat                           >  1:27494/1‑246 (MQ=255)
                                                                                                                                                                                                           cGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTg                       >  1:388276/1‑251 (MQ=255)
                                                                                                                                                                                                                           gcCCCCTCACCCTACCCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTCCTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGCGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGAACTACGGGGGCAGTATGCTAAACATCGTACTTTACGCACCAGACATTCCGCCAAATACTGGCACAATCATCCGTCt       >  2:366840/1‑251 (MQ=255)
                                                                                                                                                                                                                            cccccTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCtt      >  2:186042/1‑251 (MQ=255)
                                                                                                                                                                                                                              cccTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTg    >  1:313494/1‑251 (MQ=255)
                                                                                                                                                                                                                                cTCACCCTAACCCTCTCCTTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGCGAGGGCCGGGGAGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTgcg  >  2:342590/1‑251 (MQ=255)
                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                        
GCTGGGTCGTGCTTTCTTGTATGCGCTGGCAACAGCGGGCCAGGCGGGTGTAGCTAACCTGCTAAATCTGATCGAAAAAGAGATGAAAGTGGCGATGACGCTGACTGGCGCGAAATCGATCAGCGAAATTACGCAAGATTCGCTGGTGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATTCCGCCAAATACTGGCAACATCATCCGTCTTTGCG  >  NC_000913/3780807‑3781281

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