Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,368,219 G→A 5.1% 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 5.1% 257.9 / 3.0 96T131T (ACG→ACA) yhcGDUF1016 family protein in the PD‑(D/E)XK nuclease superfamily
Reads supporting (aligned to +/- strand):  ref base G (48/43);  new base A (2/3);  total (50/46)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

TCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCA  >  NC_000913/3368071‑3368365
                                                                                                                                                    |                                                                                                                                                  
tctcTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACg                                                                                                                                                    >  2:1033183/1‑149 (MQ=255)
 ctctGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGCCACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGc                                                                                                                                                   >  2:562445/1‑149 (MQ=255)
    tGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGt                                                                                                                                                >  1:1137850/1‑149 (MQ=255)
    tGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACAACGTCTGCAGAATACCCATCGGCAAATATTTTGCAGACGCTGt                                                                                                                                                >  2:1125559/1‑149 (MQ=255)
     tcAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTc                                                                                                                                               <  1:679092/148‑1 (MQ=255)
      cAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCt                                                                                                                                              >  1:1162049/1‑149 (MQ=255)
       aaaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTACCATCGGCAAAGATTTGGCAGACGCTGTCTg                                                                                                                                             >  2:358795/1‑149 (MQ=255)
        aaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGc                                                                                                                                            <  2:388874/149‑1 (MQ=255)
        aaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCACATCGGCAAAGATTTGGCAGACGCTGTCTGc                                                                                                                                            >  1:209015/1‑149 (MQ=255)
        aaaaaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGATGAAATTCACCAGACAATGTCTGTCGAATTAACACCATTGGGAATTACGCAGACACCGTCTGCATAAATCCCATCGGCAAAGATTTGTCAGACGCTGTATga                                                                                                                                            >  1:709554/1‑148 (MQ=255)
           aaCTTACGCCAATTGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGGCTGGAGAATTAACACCATTGGGAATTCCGCAGACACCCTCTGCAGAATTCCCATCGTCAAAGATTTGGCAGACAGCGTCTGCCaa                                                                                                                                         <  1:1173225/149‑1 (MQ=255)
           aaCTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACAGCGTCTGCCaa                                                                                                                                         <  2:302840/149‑1 (MQ=255)
              ttACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGTAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATc                                                                                                                                      >  1:743104/1‑149 (MQ=255)
              ttACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACCCCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATc                                                                                                                                      >  1:109578/1‑149 (MQ=255)
              ttACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATc                                                                                                                                      >  2:727767/1‑149 (MQ=255)
              ttACGCCAAATGAGGCTTTTTTACCTCTTATTTCAACATGTAGAAATTCACCAGACACTGTCTGGCGAATTAACTCCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTCGCAGACAGCGTCTGCCAAATc                                                                                                                                      <  1:653987/149‑1 (MQ=255)
                  gCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTc                                                                                                                                  <  1:183604/149‑1 (MQ=255)
                  gCCAAATGAGGCTTTTGTACCTCTTTTTTCAACATGTTGAAATTCACAAGACAATGTATTGCGAATAAACACAATTGGGAATTACGCAGAAAACGTCTCAAGAATTCACAAAGTCAAAGAATTAGAATACATTGTCTACAAAATcgtaa                                                                                                                                  >  2:768923/1‑145 (MQ=255)
                   ccAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTcc                                                                                                                                 <  1:983166/149‑1 (MQ=255)
                     aaaTGAGGCTTTTTTACATCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTcccc                                                                                                                               <  1:215804/149‑1 (MQ=255)
                          aGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGcc                                                                                                                          >  1:147040/1‑149 (MQ=255)
                          aGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGcc                                                                                                                          >  1:752132/1‑149 (MQ=255)
                            gCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCa                                                                                                                        >  1:170656/1‑149 (MQ=255)
                            gCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCa                                                                                                                        >  1:1163023/1‑149 (MQ=255)
                                  ttACCTCTTTTTTCAACATGTTGAAATTCAACAGACAATGTCTGTCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCTGGTCCa                                                                                                                  <  1:1115957/149‑1 (MQ=255)
                                   tACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCAc                                                                                                                 >  2:415625/1‑149 (MQ=255)
                                       tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCAAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                                             >  2:133816/1‑149 (MQ=255)
                                       tctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAAATCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCAAAATCTTTCCCCTTGCCCAGGTCCACCTAc                                                                                                             >  2:652982/1‑149 (MQ=255)
                                        ctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACg                                                                                                            >  2:628515/1‑149 (MQ=255)
                                        ctTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACg                                                                                                            >  2:44241/1‑149 (MQ=255)
                                          tttttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTc                                                                                                          >  2:864791/1‑149 (MQ=255)
                                            tttCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCg                                                                                                        >  2:867142/1‑149 (MQ=255)
                                            tttCAACAAGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGCATTGCCATCGGCAAAGAATTGGCAGCCGGTGTCTGCCAAATCTTTACCCTTGCCCACGTGCACCTACGTCCg                                                                                                        >  2:1159346/1‑149 (MQ=255)
                                              tCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGtt                                                                                                      >  2:1086767/1‑149 (MQ=255)
                                                         aaaTTCACCAGACAATGTCTGGCGAATTTACACCATTGGGAAATCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGtt                                                                                           <  2:316155/149‑1 (MQ=255)
                                                           aTTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaa                                                                                         <  2:752132/149‑1 (MQ=255)
                                                           aTTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTaa                                                                                         <  1:460161/149‑1 (MQ=255)
                                                               aCCAGTCAATGTCTGGCGAATTATCATCATTGGGTATTCGGCAGACACCGTCTGAAGAATTCCCATCGGCAAAGTTTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAAc                                                                                     <  2:953060/149‑1 (MQ=255)
                                                                ccAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACg                                                                                    <  2:59051/149‑1 (MQ=255)
                                                                 cAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACTCCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCGAATTCTTTCCCCGTGCCCGGGGACGCCTACGGGCGTTTAGTGTCCGGTatgagcgc                                                                                   >  2:178222/1‑142 (MQ=255)
                                                                 cAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGc                                                                                   <  1:405340/149‑1 (MQ=255)
                                                                 cAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGc                                                                                   <  1:823925/149‑1 (MQ=255)
                                                                      aaTGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACg                                                                              >  2:909521/1‑149 (MQ=255)
                                                                      aaTGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACg                                                                              >  2:914282/1‑149 (MQ=255)
                                                                       aTGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTGTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCACAGGTCCACCTACGTCCGTTTACTTTACGTTAAAAACGCTGACGc                                                                             >  1:834137/1‑149 (MQ=255)
                                                                       aTGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGc                                                                             <  1:534957/149‑1 (MQ=255)
                                                                             ggCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGCCGCCCGCAg                                                                       >  1:1011740/1‑149 (MQ=255)
                                                                                 aaTTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCCACGTCCGTTTACTTTCCGTTAAACACGCTGACGCCCGCAGCttt                                                                   >  1:17053/1‑149 (MQ=255)
                                                                                    tAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTAt                                                                <  1:1181795/149‑1 (MQ=255)
                                                                                    tAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTAt                                                                <  2:204060/149‑1 (MQ=255)
                                                                                     aaCACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATg                                                               <  1:636536/149‑1 (MQ=255)
                                                                                      acacCATTGGGAATTCCGCATTCACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGa                                                              <  1:703819/149‑1 (MQ=255)
                                                                                        acCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaa                                                            >  2:3125/1‑149 (MQ=255)
                                                                                        acCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTACCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaa                                                            <  1:845197/149‑1 (MQ=255)
                                                                                          cATTGGGAATTCAGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGaaaaa                                                          <  1:868340/149‑1 (MQ=255)
                                                                                            ttGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGAAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAaga                                                        <  2:650227/149‑1 (MQ=255)
                                                                                              gggAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAgaga                                                      <  2:749891/149‑1 (MQ=255)
                                                                                                  aTTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCt                                                  >  2:992280/1‑149 (MQ=255)
                                                                                                    tCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTcc                                                <  2:678335/149‑1 (MQ=255)
                                                                                                      cGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGCCGCCCGCAGGTTTTATGAAAAAGAGACGCTCCGc                                              >  2:726992/1‑149 (MQ=255)
                                                                                                      cGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGc                                              <  2:147040/149‑1 (MQ=255)
                                                                                                        cAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACGTTCCGTTAAAAACGCTGACGCCCGCAGGTTTTATGAAAAAGAGACGCTCCGCtg                                            >  2:370700/1‑149 (MQ=255)
                                                                                                         aGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCtgt                                           <  2:284201/149‑1 (MQ=255)
                                                                                                           acacCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTgg                                         <  2:728590/149‑1 (MQ=255)
                                                                                                           acacCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTACGTTTACTTTCCGTTAAAAACGCTGACGCGCGCAGCTTTTATGTTAAAGAGACGCTCCCCtgttg                                         >  2:401739/1‑147 (MQ=255)
                                                                                                              ccGTCTGCATAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGACCACCTACGTCCGTTTACTTTCCGTTAAAAACGCAGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTg                                      <  2:105544/149‑1 (MQ=255)
                                                                                                              ccGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGATGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTg                                      <  1:628515/149‑1 (MQ=255)
                                                                                                                gTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGCCGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGt                                    >  2:191555/1‑149 (MQ=255)
                                                                                                                gTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGt                                    <  2:931130/149‑1 (MQ=255)
                                                                                                                 tCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTc                                   >  2:342200/1‑149 (MQ=255)
                                                                                                                   ggCAGAATTCCCATCGGCAAAGATATGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTg                                 <  1:97053/148‑1 (MQ=255)
                                                                                                                     cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGCCGCTCCGCTGTGGCTGGTCTGtt                               >  1:565427/1‑149 (MQ=255)
                                                                                                                     cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGtt                               >  1:1019308/1‑149 (MQ=255)
                                                                                                                     cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTAGGAAAAAGAGCCGCTCCGCTGTGGCTGGTCTGtt                               >  2:864916/1‑149 (MQ=255)
                                                                                                                     cAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCAACCTACGTCCGTTTACTTTCCGTTAAAACCGCTGACGCCCGCAGCTTTTATGAAAAAGAGAAGCTCCGCTGTCGCTGGTCTGtt                               >  1:570006/1‑149 (MQ=255)
                                                                                                                     cAGAATTCCCATCGGCAAAGAGTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTCTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGtt                               <  2:919696/149‑1 (MQ=255)
                                                                                                                      aGAATTCCCGTCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTTCCCAGGTCCCCCTACGTCCGTGTACTGTCCGTTTAAACCGCTGCCGCCCCCAGCGGTTATGAAAAAGAGACGCTCCGCTGTTGCTGGTCTTTTc                              >  1:794773/1‑149 (MQ=255)
                                                                                                                      aGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTc                              <  2:170656/149‑1 (MQ=255)
                                                                                                                      aGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTc                              <  1:1086767/149‑1 (MQ=255)
                                                                                                                      aGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTc                              <  1:329561/149‑1 (MQ=255)
                                                                                                                      aGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACCTCCGGTTTCTTTCCGTTAAAAACGCGGACGCCCGCAGCTTTTAGGAAAAAGAGAAGCTCCGCTGTGGCTGGTCTGTTc                              >  1:1070581/1‑149 (MQ=255)
                                                                                                                       gAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCg                             >  1:757342/1‑149 (MQ=255)
                                                                                                                        aaTTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCGCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGc                            <  1:946981/149‑1 (MQ=255)
                                                                                                                          agcccATCGGCAAAGATTTTGTTGCCGCGGTGTGCGAAATCTTTAAACTTGATAAGGGCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCa                          <  2:1185550/147‑1 (MQ=255)
                                                                                                                            cccATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGc                        <  1:914282/149‑1 (MQ=255)
                                                                                                                            cccATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGc                        <  1:909521/149‑1 (MQ=255)
                                                                                                                             ccATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACGTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCt                       >  1:607962/1‑149 (MQ=255)
                                                                                                                                 ggggaaaGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACTTACGTCCGTTTACTTTCCGTTAAAAACGATGACGCCCGCAGCTTTAATGAATAAGAGACGCTCCGCTGTGGCTGGTATGTTCGACAGCTagag                   >  1:852840/5‑149 (MQ=255)
                                                                                                                                  ggCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGc                  <  1:652982/149‑1 (MQ=255)
                                                                                                                                   gCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCg                 >  1:340901/1‑149 (MQ=255)
                                                                                                                                     aaaGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGCCGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGCCTGGTCTGTTCGCCAGCTAGAGAGGc               >  2:198439/1‑149 (MQ=255)
                                                                                                                                       aGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCACAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCaa             <  2:1163023/149‑1 (MQ=255)
                                                                                                                                        cctgtgtcAGACGCTTGTTGCCAAATCTTTCCCCTTGACCAGGTCCACCTACTTCCGTTTACTTTCCTTTATAAACGCTGACGCTAGCAGCTTTTATGTAAAAGAGACGCTCCGCTGTGGCTGGGCTGTTCGACAGCTAGAGCGGCaaa            <  1:768923/142‑1 (MQ=255)
                                                                                                                                           ttGGCAGACAGCGTCTGCCAAATCTTTGCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCATCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTg         >  2:1101479/1‑149 (MQ=255)
                                                                                                                                              gCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGa      <  2:1162049/149‑1 (MQ=255)
                                                                                                                                              gCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGa      >  1:892686/1‑149 (MQ=255)
                                                                                                                                                  aCGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCa  <  2:1035014/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCTCTGCAAAAAACTTACGCCAAATGAGGCTTTTTTACCTCTTTTTTCAACATGTTGAAATTCACCAGACAATGTCTGGCGAATTAACACCATTGGGAATTCCGCAGACACCGTCTGCAGAATTCCCATCGGCAAAGATTTGGCAGACGCTGTCTGCCAAATCTTTCCCCTTGCCCAGGTCCACCTACGTCCGTTTACTTTCCGTTAAAAACGCTGACGCCCGCAGCTTTTATGAAAAAGAGACGCTCCGCTGTGGCTGGTCTGTTCGCCAGCTAGAGCGGCAAATTGCGACCCA  >  NC_000913/3368071‑3368365

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