Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,219 G→A 6.0% T131T (ACG→ACA yhcG → DUF1016 family protein in the PD‑(D/E)XK nuclease superfamily

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,368,2190GA 6.0% 561.6 / 15.5 167T131T (ACG→ACAyhcGDUF1016 family protein in the PD‑(D/E)XK nuclease superfamily
Reads supporting (aligned to +/- strand):  ref base G (69/87);  new base A (6/4);  total (76/91)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 170 total aligned reads displayed.

TTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAG  >  NC_000913/3368070‑3368366
                                                                                                                                                     |                                                                                                                                                   
   tctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTg                                                                                                                                                  >  1:3422095/1‑150 (MQ=255)
   tctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTg                                                                                                                                                  >  1:1784363/1‑150 (MQ=255)
      gCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCt                                                                                                                                               <  2:2281938/150‑1 (MQ=255)
        aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGc                                                                                                                                             <  1:2836934/150‑1 (MQ=255)
         aaaaaCTTCCGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGcc                                                                                                                                            <  2:455896/150‑1 (MQ=255)
          aaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCa                                                                                                                                           <  2:304164/150‑1 (MQ=255)
           aaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCaa                                                                                                                                          >  1:733003/1‑150 (MQ=255)
            aaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACTTGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGATTTCCCATCGGCAAAGATTTGGCAGACAGCGTCTGCCaaa                                                                                                                                         <  1:1722047/150‑1 (MQ=255)
            aaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCaaa                                                                                                                                         <  1:973839/150‑1 (MQ=255)
                   gCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCtttac                                                                                                                                  >  1:3689000/1‑148 (MQ=255)
                   gCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTcc                                                                                                                                  >  1:52264/1‑150 (MQ=255)
                   gCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTcc                                                                                                                                  >  1:2406933/1‑150 (MQ=255)
                    ccAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGTCGAATTAACACCATTGGGAATTCCGCAGACACCGTATGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTccc                                                                                                                                 <  2:2614468/150‑1 (MQ=255)
                    ccAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGACTTCCCATCGGCAAAGATTTGGCAGACGCGGTCTGCCAAAACTTTccc                                                                                                                                 >  1:2070888/1‑150 (MQ=255)
                    ccAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTccc                                                                                                                                 <  2:239168/150‑1 (MQ=255)
                     cAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTcccc                                                                                                                                >  2:743320/1‑150 (MQ=255)
                     cAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATATGGCAGACGCGGTCTGCCAAATCGTTcccc                                                                                                                                >  2:1390310/1‑150 (MQ=255)
                             gCTTTTTTCCCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGATTTAACCCCATGGGGAAGTCCGCAGACACCGTCGGCAGAATTCCCATCGGCAAAGATTGGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAg                                                                                                                        <  1:3923402/150‑1 (MQ=255)
                             gCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAg                                                                                                                        <  2:2110519/150‑1 (MQ=255)
                             gCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAg                                                                                                                        <  1:2043720/150‑1 (MQ=255)
                             gCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAg                                                                                                                        >  2:2860029/1‑150 (MQ=255)
                             gCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGAAGCTGTCTGCCCAATCTTTCCCCTTGCCCAg                                                                                                                        >  1:893274/1‑150 (MQ=255)
                             gCTTTTTAACCTCTTTTTTCAACACGTCGAAATTCACCAGCCATTGTCGGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAg                                                                                                                        <  1:3143371/150‑1 (MQ=255)
                                 ttttACCTCTTTTTTCAACCTGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCGGTCTGCCAAATCTTTCCCCTTGCCCAGGTcc                                                                                                                    <  2:2039612/150‑1 (MQ=255)
                                 ttttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTcc                                                                                                                    <  1:2594554/150‑1 (MQ=255)
                                 ttttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCGGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTcc                                                                                                                    <  1:844533/150‑1 (MQ=255)
                                  tttACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCa                                                                                                                   <  1:2272325/150‑1 (MQ=255)
                                        tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACg                                                                                                             >  2:976539/1‑150 (MQ=255)
                                         ctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGt                                                                                                            >  1:3131465/1‑150 (MQ=255)
                                           tttttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTcc                                                                                                          <  1:3547470/150‑1 (MQ=255)
                                              ttCAACATGTTGATATTCACCAGACAATGTCTGGCGATTTAACACCATTGTGAATTACGCAGACACCGTCTGCAGATTTCCCATGGGCAAAGATTTGGCTGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCAGtt                                                                                                       <  1:1207552/150‑1 (MQ=255)
                                              gacAACATGCTGAAATTCACCAGACACTGGCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCACAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGtt                                                                                                       <  2:2570331/148‑1 (MQ=255)
                                                  aCATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACt                                                                                                   >  1:1081375/1‑150 (MQ=255)
                                                  aCATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGGTTACt                                                                                                   >  1:1081320/1‑150 (MQ=255)
                                                  aCATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCGGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACt                                                                                                   <  2:374791/150‑1 (MQ=255)
                                                     tGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTc                                                                                                <  1:2952914/150‑1 (MQ=255)
                                                       ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCg                                                                                              >  1:3481114/1‑150 (MQ=255)
                                                       ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGATTACTTTCCg                                                                                              <  2:1777791/150‑1 (MQ=255)
                                                       ttGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCATCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATATTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCg                                                                                              >  1:1406038/1‑150 (MQ=255)
                                                        ggAAATTCACCAGACAAGGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGt                                                                                             <  2:848892/149‑1 (MQ=255)
                                                         gAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCCTCGGCAATGATTTGGCAGACGCGGTCTTCCAAATCTGTCCCCTGGGCCCGGTCCACCAACGGCCGGTTAATTTGCggt                                                                                            >  1:479811/1‑148 (MQ=255)
                                                         gAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGtt                                                                                            >  1:1971415/1‑150 (MQ=255)
                                                         gAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGtt                                                                                            >  1:479835/1‑150 (MQ=255)
                                                          aaaTTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTa                                                                                           <  2:2511043/150‑1 (MQ=255)
                                                             ttCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaaa                                                                                        <  1:3516914/150‑1 (MQ=255)
                                                             ttCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaaa                                                                                        <  1:599552/150‑1 (MQ=255)
                                                             ttCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTGTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaaa                                                                                        >  2:1752636/1‑150 (MQ=255)
                                                             ttCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTAGGCAGACGCTGTCTGCCAAATCTTTCACCTAGCCCAGGTCCACCTACGTCCGATTACTTTCCGTTaaaa                                                                                        >  2:2307929/1‑150 (MQ=255)
                                                              tCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaaaaa                                                                                       <  1:2963932/150‑1 (MQ=255)
                                                                 ccAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGc                                                                                    <  2:60812/150‑1 (MQ=255)
                                                                  cAGACAATGTCTGGCTAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCt                                                                                   <  1:1268966/150‑1 (MQ=255)
                                                                  cAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCt                                                                                   <  2:3911707/150‑1 (MQ=255)
                                                                  cAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCt                                                                                   <  2:722937/150‑1 (MQ=255)
                                                                       aaTGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGTCGc                                                                              >  1:1491044/1‑150 (MQ=255)
                                                                               gCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATATTTCCCCTTGCCCAGGTCCACCAAAGTCCGGTTACGTTGCGTTAAAAAAGCTGACGCCCGCAGCt                                                                      >  1:3003332/1‑150 (MQ=255)
                                                                                 gAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCATAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCttt                                                                    <  1:1891306/150‑1 (MQ=255)
                                                                                  aaTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGGCGCCCGCAGCtttt                                                                   >  2:2722527/1‑150 (MQ=255)
                                                                                   aTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCATGTCCACCTAGGTCCTTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                                                  <  1:3518097/150‑1 (MQ=255)
                                                                                   aTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                                                  <  1:1238950/150‑1 (MQ=255)
                                                                                   aTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                                                  <  1:2591329/150‑1 (MQ=255)
                                                                                   aTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                                                  <  2:2535553/150‑1 (MQ=255)
                                                                                   aTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTa                                                                  <  1:3266800/150‑1 (MQ=255)
                                                                                         acCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaaa                                                            <  2:448182/150‑1 (MQ=255)
                                                                                          ccATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaaaa                                                           <  2:878072/150‑1 (MQ=255)
                                                                                             ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAATGATTTGGCAGACGCTGTCTGCCAAATCTTTCGCGTTGCGCAGGTGCACCTAGGTCGGTTTACTTTCCGTTATTATCGCTGACGCGCGCAGCTTTTATGAATAtgag                                                        >  2:3583223/1‑150 (MQ=255)
                                                                                             ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAagag                                                        >  2:3814692/1‑150 (MQ=255)
                                                                                             ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAagag                                                        >  2:199096/1‑150 (MQ=255)
                                                                                              tGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAgaga                                                       >  2:1264974/1‑150 (MQ=255)
                                                                                              tGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAgaga                                                       >  2:1264981/1‑150 (MQ=255)
                                                                                                  aaTTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCt                                                   >  2:563444/1‑150 (MQ=255)
                                                                                                         cAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCtgt                                            <  1:2001315/150‑1 (MQ=255)
                                                                                                         cAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCtgt                                            <  2:3608867/150‑1 (MQ=255)
                                                                                                         cAGACACCGTCTGCAGAATTCCCATCGGAAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCtgt                                            <  2:1785204/150‑1 (MQ=255)
                                                                                                           cacacCGTCTGCAGACTCCCCCTCGGCAAACATTTGGCAGACGCGGTCTGCCAAATCTTTCCCCTGGCCCAGGTCCACCTACGTCCTTTTACTTTCCGTTAAAAACGCTGCCGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTgg                                          <  1:380415/149‑1 (MQ=255)
                                                                                                            acacCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGc                                         >  2:2641808/1‑150 (MQ=255)
                                                                                                                    tGCAGAATTCCCATCGGCAATGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAAAGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGt                                 >  2:1488591/1‑150 (MQ=255)
                                                                                                                      cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTc                               <  2:3012197/150‑1 (MQ=255)
                                                                                                                      cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTc                               <  1:3539248/150‑1 (MQ=255)
                                                                                                                      cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTc                               <  2:1653875/150‑1 (MQ=255)
                                                                                                                      cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTc                               <  1:3317989/150‑1 (MQ=255)
                                                                                                                       aGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCg                              <  1:1520308/150‑1 (MQ=255)
                                                                                                                          aTTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCa                           <  2:422787/150‑1 (MQ=255)
                                                                                                                               caaggggaaaGATTTGGCAGACCGCGTCTGCCAAATCTTTGCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTAGGAAACAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTag                      >  2:953379/8‑150 (MQ=38)
                                                                                                                               cATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTag                      <  1:824481/150‑1 (MQ=255)
                                                                                                                                 tCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTagag                    >  2:3393953/1‑150 (MQ=255)
                                                                                                                                   ggCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCg                  <  2:3220098/150‑1 (MQ=255)
                                                                                                                                     cAAAGATTTGGCAGACGCTGTCTGCTAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCTGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGc                <  2:1456204/150‑1 (MQ=255)
                                                                                                                                      aaaGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTAATTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCa               >  1:2058517/1‑150 (MQ=255)
                                                                                                                                      aaaGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTAACTTTCCGTAAAAAACGCTGACGCCCGCAGCTGTTATGAAAAAGAGACGCTCCGCTCTGGCTGGTCTGTTCGCCAGCTAGAGCGGCa               >  2:2813526/1‑150 (MQ=255)
                                                                                                                                      aaaGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCATGCCCAGGTCCACCTACGTCCGTTTACTTTACGTTAAAAAAGCAGACGCCCGCAGCTTTTATGAAAAAGATACGCTCCGCGGTGGCTGGTCTGTTCGCCACCTAGAGCGGCa               >  1:2675480/1‑150 (MQ=255)
                                                                                                                                       aaGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCaa              >  2:2260332/1‑150 (MQ=255)
                                                                                                                                        aGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCGTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCaaa             <  2:1834959/150‑1 (MQ=255)
                                                                                                                                         gATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAAt            >  2:1940782/1‑150 (MQ=255)
                                                                                                                                          aTTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTAATTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAAtt           >  1:603966/1‑150 (MQ=255)
                                                                                                                                           tttGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTg          >  1:3246731/1‑150 (MQ=255)
                                                                                                                                           tttGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTg          >  1:225549/1‑150 (MQ=255)
                                                                                                                                            ttGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGc         <  2:3422095/150‑1 (MQ=255)
                                                                                                                                            ttGGCAGACAGCGTCTGCCAAATCTTTGCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGc         >  2:2037874/1‑150 (MQ=255)
                                                                                                                                            ttGGCAGACAGCGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGc         >  2:3190008/1‑150 (MQ=255)
                                                                                                                                            ttGGCAGACAGCGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGc         >  2:1931954/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
TTCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCAG  >  NC_000913/3368070‑3368366

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