Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,781,101 A→C 18.6% intergenic (+84/‑114) lldD → / → trmL L‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,781,1010AC18.6% 127.2 / 2.5 90intergenic (+84/‑114)lldD/trmLL‑lactate dehydrogenase, FMN‑linked/tRNA Leu mC34,mU34 2'‑O‑methyltransferase, SAM‑dependent
Reads supporting (aligned to +/- strand):  ref base A (36/33);  new base C (8/10);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATT  >  NC_000913/3780953‑3781247
                                                                                                                                                     |                                                                                                                                                  
tGCAGGGGCCGTGTAAAGAGTTGCCTGCGGCAATGGCTCCCATGGCGAAAGGGAATGCGGAATATTCTTTTGCGCCCCCTCAACTTAACCCTCTCCGTCAGGGAGAGGTGACCGTTCGGAGCTGTATGTACTCCCTCACTCTGAAACGa                                                                                                                                                    <  2:647611/149‑1 (MQ=255)
ggCAGGGGCTGGGTAAAGAGTGGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGa                                                                                                                                                    <  2:996758/148‑1 (MQ=255)
 gCAGGGGCTGGGTAAAGAGTTGCCTGCGGAACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCCTCGCTGTATGTACTCCCTCACTCTGAACCGac                                                                                                                                                   >  2:602701/1‑149 (MQ=255)
  cAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGaca                                                                                                                                                  <  1:827034/149‑1 (MQ=255)
   aGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGacac                                                                                                                                                 <  1:775272/149‑1 (MQ=255)
    ggggCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTGGGCGCTGTATGTACTCCCTCACTCTGAAACGCCAcc                                                                                                                                                >  2:442951/1‑149 (MQ=255)
               aaGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATCCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACGCttttt                                                                                                                                     <  1:211543/149‑1 (MQ=255)
               aaGAGGTGCCTGAGGCACTGGCCCCCATGGCGAAAGGGAATGCGGCATCGTCGTTTGCGCCCGCTATCCGTACCCGTCTCCCTCCGGGAGTGGGGACCGTTCGGCGCTGTATGTCCTCCCTCACGCTGAAACGACACCGCACTCttttt                                                                                                                                     <  1:851617/149‑1 (MQ=255)
                  tgTTTCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATACGCGTTTGCGCCCCCTCAGCCTAACCCTCTCCCTCAGGGAGAGGGGAACGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtc                                                                                                                                  <  2:174668/148‑1 (MQ=255)
                   ggtGCCGGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtct                                                                                                                                 <  2:564301/147‑1 (MQ=255)
                    ttGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCACTCACTCTGAAACGACACCGCACTCTTTTTTtatc                                                                                                                                >  2:138601/1‑149 (MQ=255)
                    tggCCTGCGGCACTGGCTCGCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCACGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtctc                                                                                                                                <  1:538678/147‑1 (MQ=255)
                    gtGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTAGGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTtctc                                                                                                                                <  1:543089/148‑1 (MQ=255)
                     tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGCGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTcc                                                                                                                               <  2:797983/149‑1 (MQ=255)
                     tGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTcc                                                                                                                               <  2:154058/149‑1 (MQ=255)
                          gCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCGCCCTCGGGGAGAGGGGCCCGTTCGGGGCGGGCTGTACACCCGCACTCTGCAACGGCACCGGCCGCTTTTGTTCGCCCgcgc                                                                                                                          >  1:668999/1‑149 (MQ=255)
                          gCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGc                                                                                                                          <  2:91650/149‑1 (MQ=255)
                          gCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAACCCACTCCCTCAGGGAGAGCGGAACGATCGGCGCTGTATGTACAGCCGCCCTCTGAAACGACACCCCCCTCTATTTTACACCCACGc                                                                                                                          >  2:65287/1‑149 (MQ=255)
                               aCTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAAACCTATCCCTCAGGGAGAGGGGTCCGGTCGGCGCTGTATGTACTGCATCACTCAGCACCGCCCCCGCACTCGTTTTTTCTCCCTCGcccgac                                                                                                                     >  2:56917/1‑146 (MQ=255)
                                   gcacccATGGCGATAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTCACCCTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGATACGACACCGCACTCTTTTTATCTCCCTCGCCCCTCCggg                                                                                                                 <  1:233548/146‑1 (MQ=255)
                                   gCTGCCTTGGCGAACGGGAATGGGGCATAGTCGTTTGCGCGCCCTCTCCCTAACCCTCTCCCTCATGTAGAGGGGAGCGGTCGGCGCTGTAAGGACTCCCTCACTCTGTAACGAAACGGCTCTCTTTTTTTCTCCCTCGCCCCTCCgtg                                                                                                                 <  1:818653/149‑3 (MQ=255)
                                   gCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCCTCAGGGAGAGGGGGCCGTTGGGGGGTGTTTGTTCTGCCGCACAGGGATACGCCCCCGCACTCTTTTTTTGGCCCGCGCCCCTACggg                                                                                                                 >  2:217121/1‑149 (MQ=255)
                                      cccATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAACCCTCACCCTAAGGGTGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGAACCCGCACTATTTTTTTATGCCTCGCCCCACCGGGgag                                                                                                              >  2:425448/1‑149 (MQ=255)
                                        cATGGCGAAAGGGAATGCGGCATAGTCGATTGCGCCCCCTCACCCGAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGagag                                                                                                            <  1:563304/149‑1 (MQ=255)
                                         aTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGGTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCTCCCCGGAGGGGCGAgg                                                                                                           <  1:483687/149‑1 (MQ=255)
                                         aTGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCAAACCCAAACCCACTCCCACAGGGAGGGGGGCCCGGTGGGCGCGGTTTGGACTGCCAGACTCAGAAACGGCACCCGACGGGTTTTTTAGACCAAGCCCCTCCggggcgggg                                                                                                           >  2:805433/1‑144 (MQ=255)
                                          tGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAggg                                                                                                          <  2:622511/149‑1 (MQ=255)
                                          tGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCAAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTAATCCCTCACTCTGCAACGACACCGCCCGCTTTATTTATCCCTCGCCCCGCCGGGGAGAggg                                                                                                          >  1:1006231/1‑149 (MQ=255)
                                           ggagAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCACCCTAAAGGAGAGGGGCCCGTTCGGCCCTGTATGTACTCCCTCAAACTGAAACGACCCCGCACGATATTTTTCTCCCTCCCCACACCGGGGAGAGGGc                                                                                                         >  2:621909/4‑149 (MQ=255)
                                           ggCGAACGGGACGGCGGCAGAGGCGTTTGCGCCCCCTCACCCTAACCCACGCCCTCAGGGAGAGGGGACCGTTCGGCGCTGGATGTACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGc                                                                                                         <  1:722676/149‑1 (MQ=255)
                                             cGAAAGGGAATGCGGCTTAGTCGGTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCGTCACTCTGAAACGACACCGCACTGTTTTTATCGCCCACGCCCCTCCGGGGAGAGGGCCg                                                                                                       <  1:904046/149‑1 (MQ=255)
                                               aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCGGGATGTACTCCCTCACGGTGAAACGGCACCGCACGCGTTTTTTCTCCCTCGCCCCGCCGGGGGGAGGGCCggg                                                                                                     >  1:13489/1‑149 (MQ=255)
                                               aaaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                                                     <  1:150385/149‑1 (MQ=255)
                                               aaaGGCAAGGCGGCATAGACGTTGGCGCCCCCTCACCCACCCCCTCTCCCACAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCggg                                                                                                     <  1:570438/149‑1 (MQ=255)
                                                aaGGGAATGCGGCATAGTCGTTTGCGCCCCCTCCCCCAAACCCACTCCCACAGGGTGAGGGGGGCGGTGGGCGGGGTTTGGACGGCCTCACACATCAACGACACCGCCCACGTTTTTTATCCCTACACCCTCCCGGGGGcggggcgggg                                                                                                    >  2:827190/1‑149 (MQ=255)
                                                   ggAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCCCGCTTTTTTTCTCCCACGCCCCGCCGGGGAGAGGGCCGGGGTgc                                                                                                 >  2:629606/1‑148 (MQ=255)
                                                    gAATGCGGCATTGGCGTTTGCGCCCCCTCACCCTACCCCACTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAg                                                                                                <  1:942344/149‑1 (MQ=255)
                                                     aaTGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCACCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                               <  1:928844/149‑1 (MQ=255)
                                                     aaTGCGGCATAGGCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCTGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAgg                                                                                               <  2:621534/149‑1 (MQ=255)
                                                          ggCATAGTCGTTTGCGCCCCCTCCCCCGAACCCTCTCCCTCAGGGAGAGGGGCCCGTTCGGCGCTGGATGTACTCCCTCACTCTGAAACGACACCGCAGTCGTTTTTTCTCCCTCGCCCCTCCGGGGGGAGGGCCCGGGGGAGGGGaaa                                                                                          >  2:625183/1‑149 (MQ=255)
                                                           gCATAGTCGTTTGCGCCCCCTCACCCTAACCCGCTCCCTCAGGGAGAGGGGACCGTGCGGCGCTGTATGTACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGaaaa                                                                                         <  2:345335/149‑1 (MQ=255)
                                                             aTAGTCGTTGGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCGCCGGGGAGAGGGCCGGGGTGAGGGGAAAAgg                                                                                       <  2:388489/149‑1 (MQ=255)
                                                                 tCGTTTGCGCCCCCTCACCCTAACCCTCTCCCGCAGGGAGAGGGGACCGTTCGGGGCGGTATGTACTCCCTCCCGCTGAAACGCCGCCGCACTCGTTTTTTCACCCTCGCCCCGCCGGGGAGCGGGCCGGGGTGGGGGGTAAAGGCCGc                                                                                   >  2:252869/1‑149 (MQ=255)
                                                                         gcCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCa                                                                           >  2:26158/1‑149 (MQ=255)
                                                                         ccccccTCCCGCTCAGCCACCCCGTCGGGGCGAGGGGACGGTTCGTCGCTGTATTTACTCCATCACTCTGAAAGCACACCGCATTCTTTTTTTCTCCCTCGCCCCTCCGGGGTTAGGGCCGGGGTTCGGGGAAAAGGCCGCACTGCTCa                                                                           <  2:199192/148‑1 (MQ=255)
                                                                          cccccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACCCCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGGACTGCTCAc                                                                          >  1:24109/1‑149 (MQ=255)
                                                                          cccccTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCAc                                                                          >  1:20439/1‑149 (MQ=255)
                                                                               tCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtt                                                                     >  1:964262/1‑149 (MQ=255)
                                                                                 aCCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAAtttt                                                                   >  2:305218/1‑149 (MQ=255)
                                                                                 aCCCTAACCCTCACCCTCCGGGAGTGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGTAAAAGGCCGCACTGCTCACCAAtttt                                                                   <  2:964262/149‑1 (MQ=255)
                                                                                  cccTCACCCTCTCCCCCCGGGAGAGGGGTCCGTGCGGCGCGGTATGTACACCCACCCTCTCAACCCACACCCCCCTCTTTTTTTCTCCCTCCCCCCGCCGGGGAGCGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTc                                                                  <  2:74884/149‑1 (MQ=255)
                                                                                   ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACCCCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                 >  1:963818/1‑149 (MQ=255)
                                                                                   ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCGCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGCGGGCCGGGGTGAGGGGAAAACGCCGCACTGCTCACCCATTTTCt                                                                 >  1:871616/1‑149 (MQ=255)
                                                                                   ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCCCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCt                                                                 >  1:178209/1‑149 (MQ=255)
                                                                                   ccTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGATGTATGTACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTATCTCCCTCGCCCATCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCCCGGCGCACCAATTTTCt                                                                 >  2:212064/1‑149 (MQ=255)
                                                                                    cacaCCCGCCCCCTCCGGGAGCGGGGACCGGTCGGCGCTGTATGTACACCCTCACTCGGACACCCCCCCGCCCTCTTTTTTTCCCCCACCCCCCTCGGGGGAGAGGGCCGGGGGGCGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                                <  1:998289/146‑1 (MQ=255)
                                                                                    caaccccTCTCCCTCCGGGTTTGGGGACCGTTCGGGGCTTTTTGTACTCCCTCCCACTGAAACGCCACCTCACTTTTTTTTTCTCCCTCGCCCCTCCGGGGAGTGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                                <  2:57967/145‑1 (MQ=255)
                                                                                    caaaCCCTCCCCCACAGGCTGAGGGGACCGTCCGCCGCTGTTTGACCTCCCACACTCTGAACCGACCCCCCACTCTTTTTTTCACCCTCGCCCCGCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTg                                                                <  2:949679/147‑1 (MQ=255)
                                                                                     tAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGTAAAGGCCGCACTGCTCCCCAATTTTCTGc                                                               >  2:237421/1‑149 (MQ=255)
                                                                                        cccTCTCCCTCCGGGAGAGGGGACCGTTCGGCGCTGTATCTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTAt                                                            <  1:731480/149‑1 (MQ=255)
                                                                                         ccTCTCCCTCAGGGTGTGGGGACCGTTCGGCGCTGTATGAACTCCCTCACTCTTAAACTACACCGAACTCTTTTTTTATCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGTTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                           <  2:20439/149‑1 (MQ=255)
                                                                                         ccTCTCCCTCAGGGTGAGGGGTCGTTTCGGGGCTGTATGTACTCCCTAAATCTTAAACGACACCGAAATCTTTTTTTCTCCCTCGCCCCTCCGGGGAGTGGGCCGGGGTGTGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATc                                                           <  2:24109/149‑1 (MQ=255)
                                                                                         ccTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACCCCGCCCTCTGTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGCGGATCGGCCGCACTGGTCACCAAGTTTCTGCTATc                                                           >  1:409990/1‑149 (MQ=255)
                                                                                          ctctCGCTCAGGGAGAGGGGACCGTTCGGCGCTGTAGGTACTCCCTCACTCTGAAACGCGCCCGCCCTGTTTTTTTCTACCTCGCCCCTCCGGGGCGAGGGCCGGGGTGAGGGGAAAGGGCCGCCCTGCTCACCAAGTTTCTGCTATcc                                                          >  1:114076/1‑149 (MQ=255)
                                                                                            cacccTCCGGGCGTGGGGACCGTTGGGGGCTGTATTTACTCCCTCACTTTTAAACCACACCGCACTATTTTTTTCTCCCACGCCCCTCCGGGGAGTGGGCCGGGGTGCGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTg                                                        <  2:453361/147‑1 (MQ=255)
                                                                                                   ctggtgagGGGACCGGTCGGCGCTGATGCTCCTCCCTCTCCCTTACACGCCCCCGCACTGTTTTTAGCTCCCGCGCCCCCCCGGGGTGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGc                                                 <  1:953428/144‑1 (MQ=255)
                                                                                                     ggAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGCCACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGCGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCAc                                               >  1:471016/1‑149 (MQ=255)
                                                                                                           gtgACCGTTCGGCGCTGTATGAACTCCATAACTTTCAAACTACACCGCACATTTTTTTACTCCCTCGCCCCCCCGGGGAGTGGGCCGGGGTTAGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                         <  1:138601/147‑1 (MQ=255)
                                                                                                           gggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGCCACCGCACTCTTATTTTCTCCCTCGCCCCTCCGGGGATAGGGCCGTGGGGAGGGGTAAAGGCCGCAATGCTCACCAATTTTGTGCTATCCTGCACCCGCACAAAggg                                         >  2:668836/1‑149 (MQ=255)
                                                                                                           gggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                         >  1:210227/1‑149 (MQ=255)
                                                                                                           gggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggg                                         >  1:217526/1‑149 (MQ=255)
                                                                                                            ggACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGAAACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGTGAAAACGCCGCACTCATCGCCAATTTTCTGATATCCTGCCCCACCACTAAgggg                                        >  1:283494/1‑149 (MQ=255)
                                                                                                            ggACCGTGCGGCGCTGTATGAACTCCCTCCCTCTGAAACGACCCCGCACCCTTTTTTTCTCCCGCGCCCCGCCGGGGCGAGGGCCGGGGGGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAgggg                                        <  1:631354/149‑1 (MQ=255)
                                                                                                             gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTATggggg                                       >  2:69355/1‑149 (MQ=255)
                                                                                                             gACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAggggg                                       >  2:870328/1‑149 (MQ=255)
                                                                                                              aCCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGc                                      >  2:657219/1‑149 (MQ=255)
                                                                                                               ccGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTTCTCACAAATTTTCTGCTATCCTGCCCCCGCACGAAGGGGGCa                                     >  2:957499/1‑149 (MQ=255)
                                                                                                                cGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAg                                    >  2:740267/1‑149 (MQ=255)
                                                                                                                    cggcggGGTATGTCCTCCCTCCCACTCAAACCCCCCCGCACACTTTTTTTCTCCCGCCCCCCGCCGGGGAGCGGGCCGGGGTGAGGGGACAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATg                                <  2:994136/149‑1 (MQ=255)
                                                                                                                      gcgcTGTATGTACTCCCTCACTCTGAAACGACACCGCCCTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCt                              <  2:178209/149‑1 (MQ=255)
                                                                                                                        gcTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCACTCCGGGGAGAGGGCCGGGTTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTATGCTATCCTGCCCACGCAGTAAGGGGGCAGTATGCTaa                            >  2:838419/1‑149 (MQ=255)
                                                                                                                        gcTGTAAGTACTCCCTCACGCTGAAACGACACCGCACTCGTTATTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGCGTGAGGGGAACCGGCCGCCCTCCTCCCCAATTTTCTGCTATCCTGCCCCCGCTCTTAGGGGGCAGGATGCTaa                            >  2:808152/1‑149 (MQ=255)
                                                                                                                         cTGTATGTACTCCCTCACTCTCAACCCCCCCCCCACTCTTTTTTTCCCCCTCCCCCCCCGGGGGAGTGGGCGGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTaaa                           <  1:981647/149‑1 (MQ=255)
                                                                                                                                 tccccccTCACTCTGAAACGACACCCCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTa                    <  1:442951/145‑1 (MQ=255)
                                                                                                                                   ccccccCACTCTGAACGGCCACCGCACTCTTTTTTGCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGTTGTGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACt                  <  2:975393/147‑1 (MQ=255)
                                                                                                                                    tCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACtt                 >  1:620128/1‑149 (MQ=255)
                                                                                                                                    tCCCTCACTCTGAAACGACACCGCACTCGTTTTTTCTCACTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTACCCAATTTTCTGATATCATGCCCACGCACTAAGGTGGCAGTATGCTATACATCGGACtt                 >  2:923228/1‑149 (MQ=255)
                                                                                                                                      cccccctctGAAAGTCCCCCGCACTCTTTTTTTCTCCCGCGCCCCCCGGGGTTGTGGGCCGGGGGGCGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTa               <  2:529792/144‑1 (MQ=255)
                                                                                                                                      ccTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCTGGGAGAGGCCCGGGGTGAGGGGAAAAGGCCGCACTCCTCCCCAAGTTTCTGCTATCCTGCCCCCGCACGAAGGGGGCAGTATGCTAAACATCGTACttaa               >  1:120294/1‑147 (MQ=255)
                                                                                                                                       cTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGTAAAGGCCGCACTGCTGACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTAGGCTAAACCTGGTACTTTAc              >  2:232239/1‑149 (MQ=255)
                                                                                                                                        agaCTCTGAAACGCCACCGCCCTGTTTTTTGCACCCTCGCCCCTCCGGGGCGAGGGCGGGGGTTCGGGTAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACg             <  1:957499/147‑1 (MQ=255)
                                                                                                                                           ctctAAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAAc          >  2:397947/1‑149 (MQ=255)
                                                                                                                                              tGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGTGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAg       >  1:305853/1‑149 (MQ=255)
                                                                                                                                                 aaCGACACCCCACTCTTTTTTTCTCCCTCGCCCCTCAGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACATTACGAACCAGaat    >  2:119099/1‑148 (MQ=255)
                                                                                                                                                   cGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAAtt  >  1:880562/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
TGCAGGGGCTGGGTAAAGAGTTGCCTGCGGCACTGGCTCCCATGGCGAAAGGGAATGCGGCATAGTCGTTTGCGCCCCCTCACCCTAACCCTCTCCCTCAGGGAGAGGGGACCGTTCGGCGCTGTATGTACTCCCTCACTCTGAAACGACACCGCACTCTTTTTTTCTCCCTCGCCCCTCCGGGGAGAGGGCCGGGGTGAGGGGAAAAGGCCGCACTGCTCACCAATTTTCTGCTATCCTGCCCCCGCACTAAGGGGGCAGTATGCTAAACATCGTACTTTACGAACCAGAAATT  >  NC_000913/3780953‑3781247

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