Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,311,285 A→C 6.4% E66A (GAG→GCG)  tonB → membrane spanning protein in TonB‑ExbB‑ExbD transport complex

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,311,2850AC 6.4% 299.7 / 4.9 99E66A (GAG→GCG) tonBmembrane spanning protein in TonB‑ExbB‑ExbD transport complex
Reads supporting (aligned to +/- strand):  ref base A (45/47);  new base C (4/3);  total (49/50)
Fisher's exact test for biased strand distribution p-value = 7.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTC  >  NC_000913/1311139‑1311430
                                                                                                                                                  |                                                                                                                                                 
cTGCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATTTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGcc                                                                                                                                                <  2:594208/150‑1 (MQ=255)
 tGCATTCATGGTGCTGTTGTGGCGGGTCTGCGCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATATCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCg                                                                                                                                               <  1:970320/150‑1 (MQ=255)
 tGCAGTCATGGTGCTGTTGTGGCGGGTCTGCTCTCTACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCCCGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCCGCCGCCCCCGGAGCCg                                                                                                                                               >  2:262854/1‑150 (MQ=255)
  gCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGTTTTCTGTCACGATGGTTACGCCTGCTGATCTCGACCCGCCCCAAGCCGTTCAGCCGCCACCGGCGCCgg                                                                                                                                              >  2:311519/1‑150 (MQ=255)
   cATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACCTCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCggg                                                                                                                                             >  1:364260/1‑149 (MQ=255)
   cATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCACGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCggg                                                                                                                                             >  2:145792/1‑149 (MQ=255)
    aTTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCggtg                                                                                                                                            <  2:1397484/150‑1 (MQ=255)
    aTTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCggtg                                                                                                                                            <  1:1364321/150‑1 (MQ=255)
     ttCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCggtgg                                                                                                                                           <  2:1498484/150‑1 (MQ=255)
     ttCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCggtgg                                                                                                                                           <  2:1558252/150‑1 (MQ=255)
     atCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGAGGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCggtgg                                                                                                                                           <  2:1373824/149‑1 (MQ=255)
       cATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTa                                                                                                                                         >  2:719592/1‑150 (MQ=255)
       cATGGTGCTGTTGTGGCGGGGCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTa                                                                                                                                         <  1:1187252/150‑1 (MQ=255)
        aTGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTag                                                                                                                                        <  2:473029/150‑1 (MQ=255)
         tGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTaga                                                                                                                                       >  1:802150/1‑150 (MQ=255)
         gggTGCTGTTGTGGCGGGGCTGCTCTCTACCTCGGTACATCAGGTTATTGAACTCCCTGCGCCGGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGCTCCggtggcggc                                                                                                                                       <  2:1245311/149‑5 (MQ=255)
         gggTGCTGTGGTGGCGGGTCTGCTCTATACCTCGGTCCATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTGCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGCTCCggtggcggc                                                                                                                                       <  1:1685251/149‑5 (MQ=255)
          ggTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTagag                                                                                                                                      <  1:640737/150‑1 (MQ=255)
            tGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGCTCCggtggcggctga                                                                                                                                    <  2:1706582/150‑8 (MQ=255)
            tGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGcc                                                                                                                                    >  1:1430435/1‑150 (MQ=255)
             gCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCa                                                                                                                                   <  2:419460/150‑1 (MQ=255)
               tgttgtGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTCATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCCGTGGTAGAGCCAGa                                                                                                                                 >  2:1590369/1‑150 (MQ=255)
                gttgtGGCGGGTCTGCTCTATACCTAGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCGGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACACGCCGGTCAGCCGCCACCGGAGCCGGTGGGAGAGCAAGaa                                                                                                                                >  2:503036/1‑150 (MQ=255)
                   gtgGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCg                                                                                                                             <  2:771943/150‑1 (MQ=255)
                       ggggTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACTGGAGCCGGTGGTAGAGCCAGAACCGGAAc                                                                                                                         <  1:1862878/149‑1 (MQ=255)
                       cGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAAc                                                                                                                         <  1:727785/150‑1 (MQ=255)
                        gggTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAAcc                                                                                                                        <  2:905969/150‑1 (MQ=255)
                        gggTCTGCTCTATACCTCGGTACATCAGATTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGaatc                                                                                                                        <  1:1830013/150‑3 (MQ=255)
                          gTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTg                                                                                                                      >  2:1207473/1‑150 (MQ=255)
                              gCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCGTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGcc                                                                                                                  <  2:490445/150‑1 (MQ=255)
                              gCTCTATACCTCGGTACATCAGGGTCTTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGATACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGcc                                                                                                                  <  2:134172/150‑1 (MQ=255)
                                 ctATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGAt                                                                                                               >  1:700016/1‑150 (MQ=255)
                                    taCCTCTGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATccc                                                                                                            <  2:558350/150‑1 (MQ=255)
                                    taCCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGCGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATccc                                                                                                            >  1:1188687/1‑150 (MQ=255)
                                      ccTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCg                                                                                                          >  1:1730148/1‑150 (MQ=255)
                                      ccTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCg                                                                                                          >  1:1574238/1‑150 (MQ=255)
                                                tCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCaaa                                                                                                >  2:390753/1‑150 (MQ=255)
                                                tCAGGTTATTGAACTACCTGCGCCTACGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCaaa                                                                                                >  1:1418923/1‑150 (MQ=255)
                                                 cAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCaaaa                                                                                               >  1:511517/1‑150 (MQ=255)
                                                  aGGTTATTGAACTACCTGAGCCTGCGCAGCCGATATCTGTCAAGATGGTTAAGCCTGCTGATCTCGAACAGCCTCAAGCCGTTCAGCCTCCACCGGAACCGCTTGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCaaaat                                                                                              >  2:288807/1‑149 (MQ=255)
                                                   ggTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCATTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAagga                                                                                             >  1:1513996/1‑150 (MQ=255)
                                                     ttATTGAACTACCTGCGCCTGCGCAGCCGATGTCTGTCCCGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTGCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAaagaag                                                                                           <  2:950887/150‑1 (MQ=255)
                                                      tATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGCTCCCCGAACCGCCAAAAGAAGc                                                                                          >  1:143804/1‑150 (MQ=255)
                                                      tATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTCAGCCGATCCCCGACCCGCCAAAAGAAGc                                                                                          >  1:680587/1‑150 (MQ=255)
                                                        ttGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCTAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCAc                                                                                        <  2:1259539/150‑1 (MQ=255)
                                                            aCTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCggt                                                                                    <  1:368965/150‑1 (MQ=255)
                                                               aCCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGCACCTGCGCCGCTCCCCGAACCGCCAAAAGAAGCACCggtggt                                                                                 >  2:372238/1‑150 (MQ=255)
                                                               aCCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCAGCggtggg                                                                                 >  2:884384/1‑149 (MQ=255)
                                                               aCCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCggtggt                                                                                 >  1:807387/1‑150 (MQ=255)
                                                                cctgcgcctgcgCAGCCGATTTCTGTCACGATGGTTACGACTCCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGCGCCAGAACCGGAACCTGAGCCGATCCCCGAACCCCCAAAAGAAGCACCGGTGGTc                                                                                >  1:973496/1‑150 (MQ=255)
                                                                   gcgcctgcgcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCAtt                                                                             <  2:1411021/150‑1 (MQ=255)
                                                                    cgcctgcgcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTg                                                                            <  2:1063358/150‑1 (MQ=255)
                                                                      cctgcgcAGCCGATTTCTGGCACGATGGCTACGCCTGCTGATCTCGACCCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGGGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGaa                                                                          <  2:364260/150‑1 (MQ=255)
                                                                           gcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGTCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaagcc                                                                     <  1:1267041/150‑1 (MQ=255)
                                                                           gcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAatgcc                                                                     >  2:386174/1‑150 (MQ=255)
                                                                           gcAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaagcc                                                                     >  1:1705663/1‑150 (MQ=255)
                                                                              gCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGCACCGGAACCTGCGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaagccgaa                                                                  >  2:1005294/1‑150 (MQ=255)
                                                                               ccGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAaagccgaag                                                                 <  2:1050520/150‑1 (MQ=255)
                                                                                     tCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGCGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCAGaaa                                                           >  1:1427248/1‑150 (MQ=255)
                                                                                     gcTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGACCCGCCAAAAGCAGCACCGGTGGTCAGTGAAAAGCCGCAGCCGcaa                                                           >  1:527079/2‑150 (MQ=255)
                                                                                      cTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGTCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAAc                                                          >  2:731268/1‑150 (MQ=255)
                                                                                       tGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAAcc                                                         <  2:35278/150‑1 (MQ=255)
                                                                                         tCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCTGAACCGGTACCTTAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTa                                                       <  1:503036/150‑1 (MQ=255)
                                                                                            cGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGGTCAGGCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAGACCTAAGc                                                    <  1:1504454/150‑1 (MQ=255)
                                                                                                ggTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCaaa                                                >  1:376154/1‑150 (MQ=255)
                                                                                                 gTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGACCCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGACGCCGACACCTAAGCCacaa                                               >  1:130399/1‑150 (MQ=255)
                                                                                                    aCGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCg                                            <  1:1750361/150‑1 (MQ=255)
                                                                                                     cGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGa                                           <  2:210545/150‑1 (MQ=255)
                                                                                                       ccTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGACCCGGTGGTAGAGCCAGAACCGGATCCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCAATCAAAAGCCGAAGCCGAAACCTAAGCCAACACCGAAg                                         >  1:541687/1‑150 (MQ=255)
                                                                                                         tgctgATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGcc                                       <  1:807415/150‑1 (MQ=255)
                                                                                                           ctgATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCgg                                     <  1:1379302/150‑1 (MQ=255)
                                                                                                           ctgATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCgg                                     >  2:1028689/1‑150 (MQ=255)
                                                                                                               tctcGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAAGCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGACATTGAAAAGCCGAAGCCGAACCCTTCGCCAAAACCGAAGCAggagga                                 >  2:1119845/1‑146 (MQ=255)
                                                                                                                ctcGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGaaa                                >  2:1246948/1‑150 (MQ=255)
                                                                                                                 tcGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGACGCACCGGTGGTCATTGAAAAGCCGAGGCCGAAACCTAAGCCATACCCGCAGCCGGTGaaaa                               >  2:1189887/1‑150 (MQ=255)
                                                                                                                 tcGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGaaaa                               >  1:1458669/1‑150 (MQ=255)
                                                                                                                   gAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGACCCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAg                             >  2:784063/1‑150 (MQ=255)
                                                                                                                    aaCCGCCACAAGCCGGTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAgg                            <  1:1730804/150‑1 (MQ=255)
                                                                                                                     aCCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAGGGt                           <  2:1228740/150‑1 (MQ=255)
                                                                                                                     aCCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGt                           <  2:631156/150‑1 (MQ=255)
                                                                                                                     aCCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAAACTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGt                           <  1:1018362/150‑1 (MQ=255)
                                                                                                                              aGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGCGCCAGAACCGGAACCTGCGCCGCTCCCCGCACCGCCAAAAGAAGCCCCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAACACCGAGGCCGGTGAAACAGGTCCAGGaggc                  >  2:762926/1‑148 (MQ=255)
                                                                                                                               gCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGGGAAAAAGGTACAGGagcag                 <  2:1135700/150‑1 (MQ=255)
                                                                                                                               cccGGTCAGCCGCCACCGGAGCCGGGGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGagcag                 <  1:1682280/149‑1 (MQ=255)
                                                                                                                                    tCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCaaa            >  1:1558108/1‑150 (MQ=255)
                                                                                                                                        ctaccACCGGCTCCGGTCGTAGAGCCAGAGCCGGAGCCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGATAAGCCGCAGCCGAACCCAAAGCCAAACCCGAAGCCGGTG‑CAAAGGTACAGGAGCAGCCAAAAcccg       >  2:1608541/4‑150 (MQ=255)
                                                                                                                                        ccgccACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAAcgc        >  1:308529/1‑150 (MQ=255)
                                                                                                                                        ccgccACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAAcgc        >  1:736539/1‑150 (MQ=255)
                                                                                                                                        ccgccACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAAcgc        <  2:1740164/150‑1 (MQ=255)
                                                                                                                                         taccACCGGCTCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAACCCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAAcgcg       >  1:1314532/3‑150 (MQ=255)
                                                                                                                                         cgccACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAAcgcg       >  1:181395/1‑150 (MQ=255)
                                                                                                                                           ccACCGGAGCCGGTGGTAGAGCCAGAACCTGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTTGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGAt     <  2:1405967/150‑1 (MQ=255)
                                                                                                                                           ccACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGAt     <  2:1574238/150‑1 (MQ=255)
                                                                                                                                            cACCGGAGCCGGTGGTAAAGCCAGAACCGGAACCTGAGCCGATCCCCGACCCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAACCCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAACACGCGATg    >  2:1344858/1‑150 (MQ=255)
                                                                                                                                             aCCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGt   <  2:1730148/150‑1 (MQ=255)
                                                                                                                                              cgggAGCCGGTGGTAGAGCCAGAACCGTAACCTGAGCCGATCCCCGAAACGCCAAAATAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACATAAGCTAAAAGCGATGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTc  <  1:145792/148‑1 (MQ=255)
                                                                                                                                              ccGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTc  <  1:748355/150‑1 (MQ=255)
                                                                                                                                              ccGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTAATTGAAAAGCCTAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTc  <  1:1590369/150‑1 (MQ=255)
                                                                                                                                              ccGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGAACCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGTGCAGCCAAAACGCGATGTc  >  1:366739/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
CTGCATTCATGGTGCTGTTGTGGCGGGTCTGCTCTATACCTCGGTACATCAGGTTATTGAACTACCTGCGCCTGCGCAGCCGATTTCTGTCACGATGGTTACGCCTGCTGATCTCGAACCGCCACAAGCCGTTCAGCCGCCACCGGAGCCGGTGGTAGAGCCAGAACCGGAACCTGAGCCGATCCCCGAACCGCCAAAAGAAGCACCGGTGGTCATTGAAAAGCCGAAGCCGAAACCTAAGCCAAAACCGAAGCCGGTGAAAAAGGTACAGGAGCAGCCAAAACGCGATGTC  >  NC_000913/1311139‑1311430

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