Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,901,034 C→G 5.7% 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 5.7% 616.0 / 21.4 175T184T (ACG→ACCyoaEputative membrane protein/conserved protein
Reads supporting (aligned to +/- strand):  ref base C (78/87);  new base G (4/6);  total (82/93)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.38e-01

Only 100 of 184 total aligned reads displayed.

GACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACGGCGTCA  >  NC_000913/1900889‑1901181
                                                                                                                                                 |                                                                                                                                                    
gACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTg                                                                                                                                                   >  4:1217420/1‑149 (MQ=255)
gACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTg                                                                                                                                                   <  3:689751/149‑1 (MQ=255)
 aCGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgg                                                                                                                                                  >  1:1001310/1‑149 (MQ=255)
 aCGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgg                                                                                                                                                  <  2:150220/149‑1 (MQ=255)
  cgcgCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTggg                                                                                                                                                <  1:280475/149‑1 (MQ=255)
  cgcgCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTggg                                                                                                                                                <  1:249674/149‑1 (MQ=255)
  cgcgCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTggg                                                                                                                                                <  1:1537851/149‑1 (MQ=255)
   gcgcAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgggg                                                                                                                                               >  2:964730/1‑149 (MQ=255)
    cgcAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCTTGGGGt                                                                                                                                              <  2:1128088/149‑1 (MQ=255)
     gcAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTg                                                                                                                                             >  2:374672/1‑149 (MQ=255)
       aaTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCt                                                                                                                                           <  1:517919/149‑1 (MQ=255)
       aaTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCt                                                                                                                                           >  3:699537/1‑149 (MQ=255)
        aTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTg                                                                                                                                          >  3:475753/1‑149 (MQ=255)
        aTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGATAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCCCACGGTGGTGGTGCtc                                                                                                                                          <  2:1312889/149‑2 (MQ=255)
        aTCTGGTTAAACACCTCGATGATGATCGAGAGGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCAAAACCTTCTGCGACCAGACTCAGACCAATCTTTAGCAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTg                                                                                                                                          <  3:1313475/149‑1 (MQ=255)
         tCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTgg                                                                                                                                         >  2:591214/1‑149 (MQ=255)
         tCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTgg                                                                                                                                         >  2:1568053/1‑149 (MQ=255)
             gTTAAACCCTTTGTTGAAGATCGAATAGCCAATCGCGGCATGCAGGTAACCTTTAGGATTGTGGAAACCGAAACCTTCTGCCACCAGACACAGACCAATCATTAACAAGAAGCTCAGACAGAGCAACACCTCCGTGGGGTGCTGTTTaa                                                                                                                                     <  4:144220/149‑1 (MQ=37)
             gTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTaa                                                                                                                                     >  4:1565624/1‑149 (MQ=255)
             gTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCCCACGGTGGTGGTGCtctgtct                                                                                                                                     <  1:1295324/149‑7 (MQ=37)
                 aaCACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGaa                                                                                                                                 <  2:1250111/149‑1 (MQ=255)
                    acTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATcg                                                                                                                              >  1:686360/1‑149 (MQ=255)
                     cTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATcgc                                                                                                                             <  2:1265413/149‑1 (MQ=255)
                     cTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATcgc                                                                                                                             <  1:1219948/149‑1 (MQ=255)
                      ttCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATcgcg                                                                                                                            >  4:842456/1‑149 (MQ=255)
                      ttCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCCCACGGTGGTGGTGCtct                                                                                                                                         <  1:394842/136‑3 (MQ=37)
                      ttCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCGCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATcgcg                                                                                                                            <  1:1101291/149‑1 (MQ=255)
                         gatgatgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCa                                                                                                                         >  4:1398812/1‑149 (MQ=255)
                           tgatgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGc                                                                                                                       <  2:250549/149‑1 (MQ=255)
                               gatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAAc                                                                                                                                    >  1:971710/1‑132 (MQ=255)
                               gatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGtt                                                                                                                   <  3:1401009/149‑1 (MQ=255)
                               gatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGtt                                                                                                                   <  3:845751/149‑1 (MQ=255)
                                atCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGttt                                                                                                                  <  2:686360/149‑1 (MQ=255)
                                  cGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAAACATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTgg                                                                                                                <  4:270869/149‑1 (MQ=255)
                                   gagaAGCCAATCGCGGCATACATGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGa                                                                                                               <  1:1366020/149‑1 (MQ=255)
                                    agaAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGAt                                                                                                              <  3:842456/149‑1 (MQ=255)
                                    agaAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCCCACGGTGGTGGTGCtctgtct                                                                                                                                     <  1:863659/126‑7 (MQ=37)
                                           aaTCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCa                                                                                                       >  4:1229748/1‑149 (MQ=255)
                                            aTCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCaa                                                                                                      <  1:647532/149‑1 (MQ=255)
                                                  gCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAAc                                                                                                <  1:1193792/149‑1 (MQ=255)
                                                    aTACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCg                                                                                              >  4:195661/1‑149 (MQ=255)
                                                    aTACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCg                                                                                              >  4:533265/1‑149 (MQ=255)
                                                      aCAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGcc                                                                                            <  4:178554/149‑1 (MQ=255)
                                                         ggTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGTACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATc                                                                                         <  1:1591062/149‑1 (MQ=255)
                                                          gTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCg                                                                                        >  3:1387253/1‑149 (MQ=255)
                                                               cTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATc                                                                                   <  3:151941/149‑1 (MQ=255)
                                                                tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCTCCATCGCaa                                                                                     >  1:1259241/1‑146 (MQ=255)
                                                                tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCa                                                                                  >  4:1061909/1‑149 (MQ=255)
                                                                tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCa                                                                                  >  2:928098/1‑149 (MQ=255)
                                                                tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCa                                                                                  >  1:526184/1‑149 (MQ=255)
                                                                tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCa                                                                                  >  2:1378407/1‑149 (MQ=255)
                                                                    ggAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGTGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTAc                                                                              <  2:837406/149‑1 (MQ=255)
                                                                       aTGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTAccgc                                                                           <  4:908446/149‑1 (MQ=255)
                                                                          tgGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACcgccgc                                                                        <  1:1008594/149‑1 (MQ=255)
                                                                           ggAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCgccgcc                                                                       >  4:368429/1‑149 (MQ=255)
                                                                            gaaaccgaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAAc                                                                                                                                    <  2:824966/87‑1 (MQ=255)
                                                                            gaaaccgaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAAc                                                                                                                                    >  1:824966/1‑87 (MQ=255)
                                                                             aaaccgaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTggg                                                                                                                                                >  2:209106/1‑74 (MQ=255)
                                                                              aaccgaaaccTTCTGCCACCAGACTCAGACCAGTCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCcatc                                                                    >  3:858813/1‑149 (MQ=255)
                                                                              aaccgaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCcatc                                                                    <  1:374672/149‑1 (MQ=255)
                                                                                  gaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCAtcatca                                                                >  4:902671/1‑149 (MQ=255)
                                                                                     accTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCg                                                             <  3:1061909/149‑1 (MQ=255)
                                                                                        ttCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCa                                                                                           <  1:196380/116‑1 (MQ=255)
                                                                                        ttCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCa                                                          >  2:1026344/1‑149 (MQ=255)
                                                                                           tGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCTAGCAACATAACCGCCATCGCAATCACTACCGCGGCCATCATCACCGGCAGAt                                                       <  4:1332779/149‑1 (MQ=255)
                                                                                            gCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGtt                                                                                                                                       >  4:432864/1‑68 (MQ=255)
                                                                                            gCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGtt                                                                                                                                       <  3:432864/68‑1 (MQ=255)
                                                                                             ccaccaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATgg                                                     <  2:17354/149‑1 (MQ=255)
                                                                                               accaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGtt                                                   <  2:345989/149‑1 (MQ=255)
                                                                                                ccaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGCTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTa                                                  >  1:78629/1‑149 (MQ=255)
                                                                                                 caGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTaa                                                 <  3:368429/149‑1 (MQ=255)
                                                                                                   gACTCAGACCAATCATTATCAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAAcc                                               <  2:700087/149‑1 (MQ=255)
                                                                                                      tCAGACCAATCATTAACAGGACGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACTATc                                            >  4:822937/1‑149 (MQ=255)
                                                                                                      tCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGttt                                                                                                                  <  3:411248/79‑1 (MQ=255)
                                                                                                        aGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATccc                                          <  2:1001310/149‑1 (MQ=255)
                                                                                                          aCCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTa                                        <  1:1378407/149‑1 (MQ=255)
                                                                                                           ccAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTAc                                       <  3:1190841/149‑1 (MQ=255)
                                                                                                            cAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCcatcat                                                                  <  1:828936/121‑1 (MQ=255)
                                                                                                             aaTCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTg                                     >  2:876430/1‑149 (MQ=255)
                                                                                                              aTCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCa                                                                                  >  4:1029185/1‑103 (MQ=255)
                                                                                                              aTCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCa                                                                                  <  3:1029185/103‑1 (MQ=255)
                                                                                                                 aTTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGt                                 <  3:122423/149‑1 (MQ=255)
                                                                                                                  ttAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGc                                                                                             <  2:723421/88‑1 (MQ=255)
                                                                                                                  ttAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGc                                                                                             >  1:723421/1‑88 (MQ=255)
                                                                                                                     aCAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAtt                             >  3:628048/1‑149 (MQ=255)
                                                                                                                     aCAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAtt                             >  3:736246/1‑149 (MQ=255)
                                                                                                                        ggAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAAcc                                                                                               >  1:769851/1‑80 (MQ=255)
                                                                                                                        ggAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAAcc                                                                                               <  2:769851/80‑1 (MQ=255)
                                                                                                                        ggAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTAcc                          <  1:1370631/149‑1 (MQ=255)
                                                                                                                              tCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATcc                    >  4:495966/1‑149 (MQ=255)
                                                                                                                                 gacagaGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAAc                 >  4:44047/1‑149 (MQ=255)
                                                                                                                                   cagaGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACga               <  2:523168/149‑1 (MQ=255)
                                                                                                                                    agagCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACgag              >  2:1355332/1‑149 (MQ=255)
                                                                                                                                      agCACCGCCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATAACCGGCAGATGGTTAACCATCCCTGCTGCAGTAATTACCGCAACCAACGagaa            <  4:1508865/149‑1 (MQ=255)
                                                                                                                                         accaccaccGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGAc         >  2:140219/1‑149 (MQ=255)
                                                                                                                                          ccaccaccGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCATCATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACg        >  1:1579259/1‑149 (MQ=255)
                                                                                                                                           caccaccGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACgg       >  3:262456/1‑149 (MQ=255)
                                                                                                                                            accaccGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACGGc      >  1:917842/1‑149 (MQ=255)
                                                                                                                                              taccGTGGGGTGCTGTTTAACGAATCGCGTCAACGGTTTGGTTCTCAGCGACATAGCCGCCGTCGATATCACCGCCGCCGCCATCATCACCGGCAGTTGGCCAGCCGTCCCTACTGCAGTAATTATCGCATCCAACGAGAAGCCgccct    >  4:706198/2‑145 (MQ=255)
                                                                                                                                              caccGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACGGCGt    >  2:268192/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGGGGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACGGCGTCA  >  NC_000913/1900889‑1901181

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