Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,901,034 C→G 6.6% T184T (ACG→ACC)  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,0340CG 6.6% 368.6 / 15.1 106T184T (ACG→ACC) yoaEputative membrane protein/conserved protein
Reads supporting (aligned to +/- strand):  ref base C (48/51);  new base G (3/4);  total (51/55)
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.04e-01

Only 100 of 107 total aligned reads displayed.

GACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGAC  >  NC_000913/1900889‑1901174
                                                                                                                                                 |                                                                                                                                             
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                                                                                                                                        <  4:719367/148‑1 (MQ=255)
    cgcAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgg‑gg                                                                                                                                        >  3:932898/1‑148 (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)
            ggTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGCtctgt                                                                                                                                <  3:614518/148‑5 (MQ=37)
             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                                                                                                                                  <  1:683804/142‑3 (MQ=255)
                aaaCACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGCtct                                                                                                                                  >  2:683804/1‑140 (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)
                                                     tACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGc                                                                                      <  1:511827/149‑1 (MQ=255)
                                                           tAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCt                                                                                                                                    <  3:965493/97‑1 (MQ=255)
                                                           tAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCt                                                                                                                                    >  4:661092/1‑97 (MQ=255)
                                                           tAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATcg                                                                                                                       <  4:776899/110‑1 (MQ=255)
                                                           tAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATcg                                                                                                                       >  3:776899/1‑110 (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)
                                                               cTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGAAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGCCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAAt                                                                             >  4:891301/1‑148 (MQ=255)
                                                                tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAGCAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATc                                                                            <  3:799459/148‑1 (MQ=255)
                                                                tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATg                                                                                                      <  3:1022451/122‑1 (MQ=255)
                                                                tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATg                                                                                                      >  4:1022451/1‑122 (MQ=255)
                                                                tttCGGAATGTGGAAAACGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATc                                                                            <  3:1013112/148‑1 (MQ=255)
                                                                 ttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGCtc                                                                                                                                   <  4:607699/92‑2 (MQ=37)
                                                                 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)
                                                                                               accaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGtt                                            >  2:732869/1‑149 (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)
                                                                                                      tCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGa                                                                                                        <  1:599879/82‑1 (MQ=255)
                                                                                                       cAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATcc                                    <  2:205506/149‑1 (MQ=255)
                                                                                                        aGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAAt                                                                             <  2:233026/107‑1 (MQ=255)
                                                                                                        aGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAAt                                                                             >  1:233026/1‑107 (MQ=255)
                                                                                                         gACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATcc                                    <  1:655002/147‑1 (MQ=255)
                                                                                                         gACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATcc                                    >  2:655002/1‑147 (MQ=255)
                                                                                                          aCCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTa                                 >  1:392837/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)
                                                                                                                        ggAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTAc                    >  3:175825/1‑148 (MQ=255)
                                                                                                                         gAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTAcc                   >  4:898137/1‑148 (MQ=255)
                                                                                                                            gCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCAcc                                                       >  2:219542/1‑109 (MQ=255)
                                                                                                                            gCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCAcc                                                       <  1:219542/109‑1 (MQ=255)
                                                                                                                              tCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATcc             >  1:998541/1‑149 (MQ=255)
                                                                                                                               cagacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATcc             >  4:14541/1‑148 (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)
                                                                                                                                 gacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAAc          <  2:754773/149‑1 (MQ=255)
                                                                                                                                 gacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAAc          <  2:1061555/149‑1 (MQ=255)
                                                                                                                                 gacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAAc          <  2:392837/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)
                                                                                                                                                 |                                                                                                                                             
GACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGAC  >  NC_000913/1900889‑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: ‑