Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,901,029 Δ1 bp 6.7% coding (557/1557 nt) yoaE ← putative membrane protein/conserved protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,901,0290A. 6.7% 361.4 / 18.1 105coding (557/1557 nt)yoaEputative membrane protein/conserved protein
Reads supporting (aligned to +/- strand):  ref base A (48/50);  new base . (3/4);  total (51/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.99e-01

Only 100 of 108 total aligned reads displayed.

TTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGAC  >  NC_000913/1900885‑1901174
                                                                                                                                                |                                                                                                                                                  
ttGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGcaccaccac                                                                                                                                                <  2:862825/149‑1 (MQ=255)
    gACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGt                                                                                                                                             <  4:1030495/148‑1 (MQ=255)
     aCGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTg                                                                                                                                            <  4:922546/148‑1 (MQ=255)
       gcgcAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGCAGCTCAGCCAGAGCACAAGCACCGTgg‑gg                                                                                                                                        <  2:963546/149‑1 (MQ=255)
       gcgcAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgg‑gg                                                                                                                                        >  2:659843/1‑149 (MQ=255)
        cgcAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgg‑gg                                                                                                                                        >  3:932898/1‑148 (MQ=255)
        cgcAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgg‑gg                                                                                                                                        <  4:719367/148‑1 (MQ=255)
        cgcAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGt                                                                                                                                       >  2:511827/1‑149 (MQ=255)
           aaTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGc                                                                                                                                     <  3:13277/148‑1 (MQ=255)
           aaTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCt                                                                                                                                    <  1:266797/149‑1 (MQ=255)
            aTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCt                                                                                                                                    <  3:24704/148‑1 (MQ=255)
              cTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGt                                                                                                                                 <  1:837327/149‑1 (MQ=255)
               tGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAAGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGCCAGAGCACCACCACCGTGG‑GGTGCTGGtt                                                                                                                                <  2:438345/149‑1 (MQ=255)
                 gTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTa                                                                                                                               <  3:430155/148‑1 (MQ=255)
                  ttAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAAc                                                                                                                             >  2:421474/1‑149 (MQ=255)
                    aaaCACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGa                                                                                                                           >  2:667858/1‑149 (MQ=255)
                    aaaCACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGCtct                                                                                                                                  >  2:683804/1‑140 (MQ=255)
                    aaaCACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGCtct                                                                                                                                  <  1:683804/142‑3 (MQ=255)
                      acacTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCGCCAGACTCAGATCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAAt                                                                                                                         >  2:767314/1‑149 (MQ=255)
                      acacTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAAt                                                                                                                         >  2:625931/1‑149 (MQ=255)
                       cacTTCGATGATGATCGAGAAGCCAATCGCGGTATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAAt                                                                                                                         <  3:337453/148‑1 (MQ=255)
                       cacTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATc                                                                                                                        <  2:445214/149‑1 (MQ=255)
                         cTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAATCTCAGACAGAGCACCACCACCGTGG‑CGTGCTGGTTATCGAATcg                                                                                                                       <  4:579533/148‑1 (MQ=255)
                           tCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATcgcg                                                                                                                     <  4:136622/148‑1 (MQ=255)
                             gatgatgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTc                                                                                                                   <  3:342553/148‑1 (MQ=255)
                                 atgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCgg                                                                                                              <  1:55913/149‑1 (MQ=255)
                                    atCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGttt                                                                                                           <  1:872564/149‑1 (MQ=255)
                                     tCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGttt                                                                                                           <  3:926064/148‑1 (MQ=255)
                                      cGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTg                                                                                                          >  4:414128/1‑148 (MQ=255)
                                       gagaAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTgg                                                                                                         <  4:133801/148‑1 (MQ=255)
                                       gagaAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGa                                                                                                        >  1:1037330/1‑149 (MQ=255)
                                            gCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                   <  1:247321/149‑1 (MQ=255)
                                             ccAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                   >  4:1013112/1‑148 (MQ=255)
                                             ccAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCa                                                                                                   >  3:367142/1‑148 (MQ=255)
                                             ccAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAg                                                                                                  <  1:659843/149‑1 (MQ=255)
                                                         tACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGc                                                                                      <  1:511827/149‑1 (MQ=255)
                                                               tAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCt                                                                                                                                    >  4:661092/1‑97 (MQ=255)
                                                               tAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCt                                                                                                                                    <  3:965493/97‑1 (MQ=255)
                                                               tAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATcg                                                                                                                       >  3:776899/1‑110 (MQ=255)
                                                               tAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATcg                                                                                                                       <  4:776899/110‑1 (MQ=255)
                                                               tAACCTTTCGGAATGTGGAAACCGAAACCTTCTGACACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCt                                                                                                                                    <  3:661092/97‑1 (MQ=255)
                                                                aaCCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGc                                                                                >  3:1043552/1‑148 (MQ=255)
                                                                aaCCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCa                                                                               >  2:384922/1‑149 (MQ=255)
                                                                 aCCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCa                                                                               <  4:275463/148‑1 (MQ=255)
                                                                  ccTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCaa                                                                              <  3:965943/148‑1 (MQ=255)
                                                                   cTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGAAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGCCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAAt                                                                             >  4:891301/1‑148 (MQ=255)
                                                                    tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAGCAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATc                                                                            <  3:799459/148‑1 (MQ=255)
                                                                    tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATg                                                                                                      <  3:1022451/122‑1 (MQ=255)
                                                                    tttCGGAATGTGGAAAACGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATc                                                                            <  3:1013112/148‑1 (MQ=255)
                                                                     ttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGCtc                                                                                                                                   >  3:607699/1‑91 (MQ=37)
                                                                      tCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACt                                                                         >  1:205506/1‑149 (MQ=255)
                                                                        ggAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTa                                                                        >  3:51647/1‑148 (MQ=255)
                                                                          aaTGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTAcc                                                                      >  3:185047/1‑148 (MQ=255)
                                                                           aTGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTAccgc                                                                    >  1:453303/1‑149 (MQ=255)
                                                                              tgGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTAccgccg                                                                  <  4:932898/148‑1 (MQ=255)
                                                                              tgGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GCTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTAccgccg                                                                  >  4:635796/1‑148 (MQ=255)
                                                                                gaaaccgaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCgccgcc                                                                >  3:165577/1‑148 (MQ=255)
                                                                                   accgaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCcatca                                                            <  1:936178/149‑1 (MQ=255)
                                                                                     cgaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCatcatc                                                          >  2:803380/1‑149 (MQ=255)
                                                                                            ttCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGc                                                                                                                                     >  3:409059/1‑67 (MQ=37)
                                                                                            ttCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGc                                                                                                                                     <  4:409059/67‑1 (MQ=37)
                                                                                             tCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGCTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCggta                                                   >  4:908349/1‑146 (MQ=255)
                                                                                                  caccaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATgg                                              <  4:51647/148‑1 (MQ=255)
                                                                                                     cacaCTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTa                                           <  3:251544/145‑1 (MQ=255)
                                                                                                      aGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTaa                                          >  4:820757/1‑148 (MQ=255)
                                                                                                      aGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAAc                                         <  1:384922/149‑1 (MQ=255)
                                                                                                        aCTCAGACGAATCATTAACAGGAAGCGCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCAATACAGCCGCCATCATCACCGGCAGATGGTTGACCa                                       >  2:602844/1‑149 (MQ=255)
                                                                                                          tCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGa                                                                                                        >  2:599879/1‑82 (MQ=255)
                                                                                                           cAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATcc                                    <  2:205506/149‑1 (MQ=255)
                                                                                                            aGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAAt                                                                             >  1:233026/1‑107 (MQ=255)
                                                                                                             gACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATcc                                    >  2:655002/1‑147 (MQ=255)
                                                                                                             gACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATcc                                    <  1:655002/147‑1 (MQ=255)
                                                                                                              aCCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCt                                  <  3:414128/148‑1 (MQ=255)
                                                                                                              aCCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTa                                 >  1:392837/1‑149 (MQ=255)
                                                                                                                  aTCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGc                             >  2:90786/1‑149 (MQ=255)
                                                                                                                      ttAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGt                          >  3:414657/1‑148 (MQ=255)
                                                                                                                       tAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTa                         <  3:247435/148‑1 (MQ=255)
                                                                                                                         aCAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAtt                      >  2:1028844/1‑149 (MQ=255)
                                                                                                                         aCAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAt                       >  4:1022391/1‑148 (MQ=255)
                                                                                                                          cAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAtt                      >  4:625950/1‑148 (MQ=255)
                                                                                                                            ggAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTAc                    >  4:745238/1‑148 (MQ=255)
                                                                                                                             gAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTAcc                   >  4:898137/1‑148 (MQ=255)
                                                                                                                                gCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCAcc                                                       <  1:219542/109‑1 (MQ=255)
                                                                                                                                gCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCAcc                                                       >  2:219542/1‑109 (MQ=255)
                                                                                                                                  tCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATcc             >  1:998541/1‑149 (MQ=255)
                                                                                                                                   cagacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATcc             >  4:10250/1‑148 (MQ=255)
                                                                                                                                   cagacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATcc             >  4:371414/1‑148 (MQ=255)
                                                                                                                                   cagacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATcc             >  3:855761/1‑148 (MQ=255)
                                                                                                                                   cagacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATcc             >  4:14541/1‑148 (MQ=255)
                                                                                                                                     gacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAAc          <  2:1061555/149‑1 (MQ=255)
                                                                                                                                     gacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAAc          <  2:392837/149‑1 (MQ=255)
                                                                                                                                     gacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAAc          <  2:754773/149‑1 (MQ=255)
                                                                                                                                       cagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACg         <  3:745238/148‑1 (MQ=255)
                                                                                                                                         gagCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACgag       >  3:717607/1‑148 (MQ=255)
                                                                                                                                          agCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACgaga      <  4:414657/148‑1 (MQ=255)
                                                                                                                                          agCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACgaga      <  4:1043552/148‑1 (MQ=255)
                                                                                                                                          agCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGagaa     <  1:1028844/149‑1 (MQ=255)
                                                                                                                                          agCACCACCACCGTGG‑GGTGCTAGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTGCCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGagaa     <  1:964103/149‑1 (MQ=255)
                                                                                                                                            caccaccacCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGagaaga   <  2:82457/149‑1 (MQ=255)
                                                                                                                                             accaccaccGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGAc  >  1:23804/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
TTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGAC  >  NC_000913/1900885‑1901174

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑