Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,015,912 T→A 16.5% intergenic (+30/+40) rmf → / ← fabA ribosome modulation factor/beta‑hydroxydecanoyl thioester dehydrase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,015,9120TA16.5% 451.0 / 73.1 181intergenic (+30/+40)rmf/fabAribosome modulation factor/beta‑hydroxydecanoyl thioester dehydrase
Reads supporting (aligned to +/- strand):  ref base T (70/81);  new base A (17/13);  total (87/94)
Fisher's exact test for biased strand distribution p-value = 3.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 183 total aligned reads displayed.

CTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAG  >  NC_000913/1015814‑1016011
                                                                                                  |                                                                                                   
    aTCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATTGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCgg                                                                                                >  1:3790593/1‑100 (MQ=255)
    aTCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCgg                                                                                                >  1:1535097/1‑100 (MQ=255)
    aTCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCgg                                                                                                >  1:2206840/1‑100 (MQ=255)
    aTCAAAGGTCACAATGGCTGGGAGGCGGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCgg                                                                                                <  2:1296091/100‑1 (MQ=255)
      cAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAg                                                                                              >  1:2139142/1‑100 (MQ=255)
       aaaGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAgg                                                                                             <  1:3429108/100‑1 (MQ=255)
         aGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGtt                                                                                           >  2:1119601/1‑100 (MQ=255)
         aGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGtt                                                                                           <  2:727680/100‑1 (MQ=255)
         aGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGtt                                                                                           >  2:2712475/1‑100 (MQ=255)
         aGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGGGGtt                                                                                           <  1:3406966/100‑1 (MQ=255)
          ggTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGttt                                                                                          <  2:3178404/100‑1 (MQ=255)
          ggTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGttt                                                                                          >  1:2792901/1‑100 (MQ=255)
           gTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAGGTTTc                                                                                         <  1:1812137/100‑1 (MQ=255)
            tCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAGGTTTct                                                                                        <  1:3192852/100‑2 (MQ=255)
              acaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGcc                                                                                      <  1:1733427/100‑1 (MQ=255)
              acaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGcc                                                                                      <  2:3207165/100‑1 (MQ=255)
              acaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAGGTTTcttt                                                                                      <  1:2477196/100‑4 (MQ=255)
              acaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCCGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAGGTTTcttt                                                                                      <  1:1453436/100‑4 (MQ=255)
              acaATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGATAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGc                                                                                       <  2:3633668/99‑1 (MQ=255)
               caATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCt                                                                                     >  1:224140/1‑100 (MQ=255)
                 aTGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCttt                                                                                   >  2:613559/1‑100 (MQ=255)
                   ggCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTg                                                                                 <  1:1006049/100‑1 (MQ=255)
                   ggCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTg                                                                                 >  1:1736273/1‑100 (MQ=255)
                   ggCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTg                                                                                 <  1:1625756/100‑1 (MQ=255)
                   ggCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCAATGCGGAGGTTTctttttaaa                                                                                 <  1:937540/100‑9 (MQ=38)
                    gCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGa                                                                                >  1:1672415/1‑100 (MQ=255)
                     cTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGAt                                                                               >  2:555298/1‑100 (MQ=255)
                         gAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATActc                                                                           >  1:1841574/1‑100 (MQ=255)
                           ggCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACtctg                                                                         >  2:2978331/1‑100 (MQ=255)
                            gCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACtctgt                                                                        >  2:3873082/1‑100 (MQ=255)
                             cTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACtctgtc                                                                       <  1:2646261/100‑1 (MQ=255)
                              tggcttctcgCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACtctgtct                                                                      >  1:1111177/10‑100 (MQ=255)
                              tGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGagatcggaaga                                                                      >  2:751369/1‑90 (MQ=255)
                                  gagaAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGAtt                                                                  >  1:3172599/1‑100 (MQ=255)
                                   agaAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATta                                                                 <  2:1168546/100‑1 (MQ=255)
                                      aGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATaa                                                              <  1:2813115/100‑1 (MQ=255)
                                       gCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAAt                                                             >  1:242316/1‑100 (MQ=255)
                                            ggCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGaa                                                        <  2:1353970/100‑1 (MQ=255)
                                             gCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAg                                                       >  1:815991/1‑100 (MQ=255)
                                             gCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAg                                                       >  1:2356675/1‑100 (MQ=255)
                                               ggACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGc                                                     >  2:301931/1‑100 (MQ=255)
                                                gACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCa                                                    <  2:2201223/100‑1 (MQ=255)
                                                   aGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGacagagtat                                                 <  2:264257/100‑9 (MQ=255)
                                                   aGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGAc                                                 <  2:2551364/100‑1 (MQ=255)
                                                    gggTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACg                                                >  1:3354523/1‑100 (MQ=255)
                                                     ggTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTATGTCTGATTATAATCAGAAGGCAGACGt                                               >  2:3928038/1‑100 (MQ=255)
                                                      gtagtaATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTa                                              >  2:2087885/1‑100 (MQ=255)
                                                        agtaATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATc                                            >  2:661340/1‑100 (MQ=255)
                                                          taATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGTCGTATCCt                                          >  2:344169/1‑100 (MQ=255)
                                                           aaTGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTg                                         <  2:1551637/100‑1 (MQ=255)
                                                           aaTGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTg                                         >  2:1584546/1‑100 (MQ=255)
                                                           aaTGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTg                                         <  2:3387811/100‑1 (MQ=255)
                                                             tGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGa                                       <  2:1496876/100‑1 (MQ=255)
                                                               gCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAAc                                     >  1:1130677/1‑100 (MQ=255)
                                                                 cTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAg                                   <  2:1002628/100‑1 (MQ=255)
                                                                        tGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACt                            >  2:3103847/1‑100 (MQ=255)
                                                                         ntctcTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACtt                           >  2:3654246/2‑100 (MQ=255)
                                                                            tctTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCa                        <  1:3743216/100‑1 (MQ=255)
                                                                            tctTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCa                        >  1:2260691/1‑100 (MQ=255)
                                                                            tctTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCa                        >  2:1773353/1‑100 (MQ=255)
                                                                             ctTTAAAAAGAAACCTCGGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAg                       >  2:1383679/1‑100 (MQ=255)
                                                                             ctTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAg                       <  2:1841565/100‑1 (MQ=255)
                                                                              tttAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAgg                      <  2:2925593/100‑1 (MQ=255)
                                                                              tttAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAgg                      >  1:102867/1‑100 (MQ=255)
                                                                              tttAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAgg                      >  1:359200/1‑100 (MQ=255)
                                                                                tAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTATCTGATTATAATCAGAAGTCAGACGTATCCTGGAACAGACCGACTTTCa                        >  2:1484175/1‑96 (MQ=255)
                                                                                tAAAAAGAAACCTCCGCATTGAGGAGGTTTCGCCCTTTGATACTCTGTCTGATTATAATCATAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGt                     <  2:640368/99‑1 (MQ=255)
                                                                                 aaaaaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCg                   >  1:2528639/1‑100 (MQ=255)
                                                                                 aaaaaGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCg                   >  1:584798/1‑100 (MQ=255)
                                                                                  aggcgAAACCTCCGCAATGCGGAGGTTTCTCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGc                  >  2:2135182/5‑100 (MQ=38)
                                                                                  aaaaGAAACCTCCGCATTGCGGAGGTTTCGTCTTTTGATTCTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACTGACTTTCAGGTCg                   <  2:466856/99‑1 (MQ=255)
                                                                                   ggcgAAACCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCt                 >  1:3271223/4‑100 (MQ=255)
                                                                                   ggcgAAACCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCt                 >  2:3367936/4‑100 (MQ=255)
                                                                                    gcgAAACCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGCCTTTCAGGTCGCt                 >  2:240937/3‑99 (MQ=255)
                                                                                    gcgAAACCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTg                >  1:2048298/3‑100 (MQ=255)
                                                                                    gcgAAACCTCCGCAATGCGGAGGTTTCGCCTTTTGATACCCTGTCTGATTATCATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGACGCTg                >  2:2241362/3‑100 (MQ=38)
                                                                                     cgAAACCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGTCGTATCCTGGAACAGACCGACTTTCAGGTCGCt                 >  2:776679/2‑98 (MQ=255)
                                                                                     cgAAACCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTgg               >  1:2838605/2‑100 (MQ=255)
                                                                                     aGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTgg               <  1:1377145/100‑1 (MQ=255)
                                                                                      gAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGc              >  1:1178463/1‑100 (MQ=255)
                                                                                        agccTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCgg            <  2:2889773/98‑1 (MQ=255)
                                                                                        aaCCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCgg            >  1:1903659/1‑100 (MQ=255)
                                                                                        aaCCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCgg            <  2:1096040/100‑1 (MQ=255)
                                                                                        aaCCTCCGCAATGCGGAGGTTTCTCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCgg            <  1:2135193/100‑1 (MQ=255)
                                                                                        aaCCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGCCTGATTATAATCAGAAGGCGGGCGTATCATGAAACAGACCTAGTTTCa                        >  2:3094960/1‑88 (MQ=37)
                                                                                         aCCTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGt           >  1:2609577/1‑100 (MQ=255)
                                                                                          ccTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCgtctgatcta          >  1:2914296/1‑91 (MQ=255)
                                                                                          ccTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGta          >  1:3932910/1‑100 (MQ=255)
                                                                                          ccTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGt           <  1:3522897/99‑1 (MQ=255)
                                                                                          ccTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCATACATATCATAGAACAGACCGACTTTCAGGTCGCTGACGGta          >  1:321935/1‑100 (MQ=255)
                                                                                          ccTCCGCAATGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGta          <  2:2048287/100‑1 (MQ=255)
                                                                                           cTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGtat         >  1:1119428/1‑100 (MQ=255)
                                                                                            tCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGtata        <  1:842972/100‑1 (MQ=255)
                                                                                              cGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGa      >  2:361088/1‑100 (MQ=255)
                                                                                              cGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGa      <  2:359200/100‑1 (MQ=255)
                                                                                               gCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGAt     <  2:2055538/100‑1 (MQ=255)
                                                                                                cATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATc    <  1:361088/100‑1 (MQ=255)
                                                                                                cATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATc    <  2:1417515/100‑1 (MQ=255)
                                                                                                cATTGCCGAGGTTTCGAATTTTGATACCCTGTCTGATACTAATCAGAAGGCAGACGTATCTTGGATCATACTGACTTTTAGGTCGCTGGCTGTATAGATc    <  2:325123/100‑1 (MQ=255)
                                                                                                 aTTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCa   <  1:2858645/100‑1 (MQ=255)
                                                                                                  |                                                                                                   
CTGAATCAAAGGTCACAATGGCTGGGAGGCTGGCGAGAAGCCATGGCGGACAGGGTAGTAATGGCCTGATTCTGTCTCTTTAAAAAGAAACCTCCGCATTGCGGAGGTTTCGCCTTTTGATACTCTGTCTGATTATAATCAGAAGGCAGACGTATCCTGGAACAGACCGACTTTCAGGTCGCTGGCGGTATAGATCAG  >  NC_000913/1015814‑1016011

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