Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,901,029 Δ1 bp 5.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. 5.7% 589.3 / 27.0 174coding (557/1557 nt)yoaEputative membrane protein/conserved protein
Reads supporting (aligned to +/- strand):  ref base A (76/85);  new base . (4/6);  total (82/92)
Fisher's exact test for biased strand distribution p-value = 7.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.07e-01

Only 100 of 185 total aligned reads displayed.

AAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACGGC  >  NC_000913/1900882‑1901177
                                                                                                                                                   |                                                                                                                                                     
aaGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGcaccac                                                                                                                                                      <  4:995089/149‑1 (MQ=255)
  gTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGcaccacca                                                                                                                                                    >  1:150220/1‑149 (MQ=255)
       gACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTg                                                                                                                                               <  3:689751/149‑1 (MQ=255)
       gACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTg                                                                                                                                               >  4:1217420/1‑149 (MQ=255)
        aCGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgg                                                                                                                                              <  2:150220/149‑1 (MQ=255)
        aCGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgg                                                                                                                                              >  1:1001310/1‑149 (MQ=255)
         cgcgCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgg‑g                                                                                                                                            <  3:1217420/149‑1 (MQ=255)
         cgcgCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgg‑g                                                                                                                                            <  3:40966/149‑1 (MQ=255)
          gcgcAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgg‑gg                                                                                                                                           >  2:964730/1‑149 (MQ=255)
           cgcAATCTGGTTAATCACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGt                                                                                                                                          <  4:162727/149‑1 (MQ=255)
           cgcAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCTTGG‑GGt                                                                                                                                          <  2:1128088/149‑1 (MQ=255)
            gcAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTg                                                                                                                                         >  1:947705/1‑149 (MQ=255)
            gcAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTg                                                                                                                                         >  2:374672/1‑149 (MQ=255)
              aaTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCt                                                                                                                                       >  3:699537/1‑149 (MQ=255)
              aaTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCt                                                                                                                                       >  1:250549/1‑149 (MQ=255)
               aTCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTg                                                                                                                                      >  3:475753/1‑149 (MQ=255)
                tCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTgg                                                                                                                                     >  2:591214/1‑149 (MQ=255)
                tCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGCtct                                                                                                                                     >  2:1244992/1‑147 (MQ=38)
                   ggTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGCtctgtc                                                                                                                                  <  2:1275966/149‑6 (MQ=255)
                    gTTAAACCCTTTGTTGAAGATCGAATAGCCAATCGCGGCATGCAGGTAACCTTTAGGATTGTGGAAACCGAAACCTTCTGCCACCAGACACAGACCAATCATTAACAAGAAGCTCAGACAGAGCAACACCTCCGTGG‑GGTGCTGTTTaa                                                                                                                                 <  4:144220/149‑1 (MQ=37)
                    gTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTaa                                                                                                                                 >  4:1565624/1‑149 (MQ=255)
                    gTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGCtctgtct                                                                                                                                 <  1:1295324/149‑7 (MQ=37)
                        aaCACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGaa                                                                                                                             >  1:978173/1‑149 (MQ=255)
                           acTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATcg                                                                                                                          >  1:686360/1‑149 (MQ=255)
                            cTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATcgc                                                                                                                         <  1:1219948/149‑1 (MQ=255)
                             ttCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCGCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATcgcg                                                                                                                        <  1:1101291/149‑1 (MQ=255)
                               cGATGATGATCGAGAAGCTAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAATCTCAGACAGTGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTc                                                                                                                      >  4:122423/1‑149 (MQ=255)
                                gatgatgatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCa                                                                                                                     >  4:1398812/1‑149 (MQ=255)
                                      gatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGtt                                                                                                               <  3:1401009/149‑1 (MQ=255)
                                      gatCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGtt                                                                                                               <  3:845751/149‑1 (MQ=255)
                                       atCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGttt                                                                                                              <  4:83340/149‑1 (MQ=255)
                                       atCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGttt                                                                                                              <  2:1339144/149‑1 (MQ=255)
                                       atCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGttt                                                                                                              <  2:686360/149‑1 (MQ=255)
                                         cGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAAACATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTgg                                                                                                            <  4:270869/149‑1 (MQ=255)
                                          gagaAGCCAATCGCGGCATACATGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGa                                                                                                           <  1:1366020/149‑1 (MQ=255)
                                          gagaAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGa                                                                                                           >  2:369749/1‑149 (MQ=255)
                                           agaAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGCtctgtct                                                                                                                                 <  1:863659/126‑7 (MQ=37)
                                           agaAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACC‑CCACGGTGGTGGTGCtctgtct                                                                                                                                 >  2:863659/1‑120 (MQ=37)
                                              aGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGcc                                                                                                       <  3:781793/149‑1 (MQ=255)
                                                ccAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCCCCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAg                                                                                                     <  4:896183/149‑1 (MQ=255)
                                                   aTCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCaa                                                                                                  <  1:647532/149‑1 (MQ=255)
                                                         gCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAAc                                                                                            <  1:1193792/149‑1 (MQ=255)
                                                           aTACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCg                                                                                          >  4:533265/1‑149 (MQ=255)
                                                           aTACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCg                                                                                          >  4:195661/1‑149 (MQ=255)
                                                            tACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGc                                                                                         >  2:1546233/1‑149 (MQ=255)
                                                               aGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTgg                                                                                                            <  2:1589108/127‑1 (MQ=255)
                                                               aGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTgg                                                                                                            >  1:1589108/1‑127 (MQ=255)
                                                               aGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGa                                                                                                           <  1:907075/128‑1 (MQ=255)
                                                               aGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGa                                                                                                           >  2:907075/1‑128 (MQ=255)
                                                                 gTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCg                                                                                    >  3:1387253/1‑149 (MQ=255)
                                                                      cTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATc                                                                               >  3:296789/1‑149 (MQ=255)
                                                                      cTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATc                                                                               <  3:151941/149‑1 (MQ=255)
                                                                       tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCTCCATCGCaa                                                                                 >  1:1259241/1‑146 (MQ=255)
                                                                       tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCaa                                                                                 <  2:1259241/146‑1 (MQ=255)
                                                                       tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCa                                                                              >  4:1061909/1‑149 (MQ=255)
                                                                       tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCa                                                                              >  1:526184/1‑149 (MQ=255)
                                                                       tttCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCa                                                                              >  2:1008594/1‑149 (MQ=255)
                                                                                 tgGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACcgccgc                                                                    <  1:1008594/149‑1 (MQ=255)
                                                                                  ggAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCgccgcc                                                                   >  4:368429/1‑149 (MQ=255)
                                                                                   gaaaccgaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAAc                                                                                                                                >  1:824966/1‑87 (MQ=255)
                                                                                   gaaaccgaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAAc                                                                                                                                <  2:824966/87‑1 (MQ=255)
                                                                                    aaaccgaaaccTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTgg‑g                                                                                                                                            <  1:209106/74‑1 (MQ=255)
                                                                                     aaccgaaaccTTCTGCCACCAGACTCAGACCAGTCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCcatc                                                                >  3:858813/1‑149 (MQ=255)
                                                                                            accTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCg                                                         <  3:1061909/149‑1 (MQ=255)
                                                                                               ttCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCa                                                                                       <  1:196380/116‑1 (MQ=255)
                                                                                               ttCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCa                                                                                       >  2:196380/1‑116 (MQ=255)
                                                                                                  tGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCTAGCAACATAACCGCCATCGCAATCACTACCGCGGCCATCATCACCGGCAGAt                                                   <  4:1332779/149‑1 (MQ=255)
                                                                                                   gCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGtt                                                                                                                                   >  4:432864/1‑68 (MQ=255)
                                                                                                   gCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGtt                                                                                                                                   <  3:432864/68‑1 (MQ=255)
                                                                                                      accaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGtt                                               <  1:369749/149‑1 (MQ=255)
                                                                                                       ccaGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGCTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTa                                              >  1:78629/1‑149 (MQ=255)
                                                                                                        caGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTaa                                             <  1:1568053/149‑1 (MQ=255)
                                                                                                                 aCCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTa                                    <  1:1378407/149‑1 (MQ=255)
                                                                                                                  ccAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTAc                                   <  3:1190841/149‑1 (MQ=255)
                                                                                                                   cAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTATCGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTAc                                                                          <  2:1366689/109‑1 (MQ=255)
                                                                                                                   cAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCcatcat                                                              <  1:828936/121‑1 (MQ=255)
                                                                                                                   cAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCcatcat                                                              >  2:828936/1‑121 (MQ=255)
                                                                                                                     aTCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCa                                                                              <  3:1029185/103‑1 (MQ=255)
                                                                                                                     aTCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCa                                                                              >  4:1029185/1‑103 (MQ=255)
                                                                                                                        aTTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGCTAACCAGCCCTACTGCGGt                             >  3:1247116/1‑149 (MQ=255)
                                                                                                                         ttAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGc                                                                                         >  1:723421/1‑88 (MQ=255)
                                                                                                                         ttAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGc                                                                                         <  2:723421/88‑1 (MQ=255)
                                                                                                                            aCAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAtt                         <  1:623752/149‑1 (MQ=255)
                                                                                                                            aCAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAtt                         >  3:628048/1‑149 (MQ=255)
                                                                                                                               ggAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAAcc                                                                                           >  1:769851/1‑80 (MQ=255)
                                                                                                                               ggAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAAcc                                                                                           <  2:769851/80‑1 (MQ=255)
                                                                                                                               ggAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTAcc                      <  2:21956/149‑1 (MQ=255)
                                                                                                                               ggAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTAcc                      <  1:1370631/149‑1 (MQ=255)
                                                                                                                                gAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCg                     <  1:51466/149‑1 (MQ=255)
                                                                                                                                  aGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCCCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCa                   <  4:904894/149‑1 (MQ=255)
                                                                                                                                     tCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATcc                >  4:495966/1‑149 (MQ=255)
                                                                                                                                     tCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATcc                >  1:389587/1‑149 (MQ=255)
                                                                                                                                       agacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCaa              >  3:72580/1‑149 (MQ=255)
                                                                                                                                        gacagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAAc             >  4:44047/1‑149 (MQ=255)
                                                                                                                                          cagaGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACga           <  2:523168/149‑1 (MQ=255)
                                                                                                                                             agCACCGCCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATAACCGGCAGATGGTTAACCATCCCTGCTGCAGTAATTACCGCAACCAACGagaa        <  4:1508865/149‑1 (MQ=255)
                                                                                                                                             agCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGagaa        <  4:72580/149‑1 (MQ=255)
                                                                                                                                               caccaccacCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCACCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAAATACCGCATCCAACGagaaga      <  2:323985/149‑1 (MQ=255)
                                                                                                                                                accaccaccGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGAc     >  2:140219/1‑149 (MQ=255)
                                                                                                                                                  caccaccGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACgg   >  3:262456/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AAGTTGCGACGCGCAATCTGGTTAAACACTTCGATGATGATCGAGAAGCCAATCGCGGCATACAGGTAACCTTTCGGAATGTGGAAACCGAAACCTTCTGCCACCAGACTCAGACCAATCATTAACAGGAAGCTCAGACAGAGCACCACCACCGTGG‑GGTGCTGGTTAACGAATCGCGTCAGCGGTTTGGATGCCAGCAACATAACCGCCATCGCAATCACTACCGCCGCCATCATCACCGGCAGATGGTTAACCATCCCTACTGCAGTAATTACCGCATCCAACGAGAAGACGGC  >  NC_000913/1900882‑1901177

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