Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,434,094 G→A 12.6% R216R (CGC→CGT)  dedA ← DedA family protein DedA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,434,0940GA12.6% 252.9 / 23.2 99R216R (CGC→CGT) dedADedA family protein DedA
Reads supporting (aligned to +/- strand):  ref base G (46/38);  new base A (6/6);  total (53/46)
Fisher's exact test for biased strand distribution p-value = 7.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.27e-01

CTGGTTCGATGCCAGTTTTATCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTAAAAAGCAGTACC  >  NC_000913/2433948‑2434241
                                                                                                                                                  |                                                                                                                                                   
cTGGTTCGATGCCAGTTTTATCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGt                                                                                                                                                  >  1:181603/1‑150 (MQ=255)
     tCGATGCCAGTTTTATCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGc                                                                                                                                             >  1:1233156/1‑150 (MQ=255)
      cGATGCCAGTTTTATCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCa                                                                                                                                            >  1:1465893/1‑150 (MQ=255)
      cGATGCCAGTTTTATCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCa                                                                                                                                            >  2:343635/1‑150 (MQ=255)
       gATGCCAGTTTTATCTTTGGGGACACATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAg                                                                                                                                           >  1:1574579/1‑150 (MQ=255)
         tGCCAGTTTTATCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCACGCGCCGCaaaa                                                                                                                                         <  1:1029785/150‑4 (MQ=255)
         tGCCAGTTTTATCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCACGCGCCGCaaaa                                                                                                                                         <  1:374020/150‑4 (MQ=255)
          gCCAGTTTTATCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATACCAGCCCGAAAGTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCACGCGCCGCaaaat                                                                                                                                        <  2:456957/150‑5 (MQ=255)
          gCCAGTTTTATCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCCCGCGCCGCaaaat                                                                                                                                        <  1:1531491/150‑5 (MQ=255)
          gCCAGTTTTATCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCACGCGCCGCaaaat                                                                                                                                        <  1:497143/150‑5 (MQ=255)
          gCCAGTTTTATCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCACGCGCCGCaaaat                                                                                                                                        <  1:996278/150‑5 (MQ=255)
           ccAGTTTTATCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGcgcg                                                                                                                                       >  2:1116934/1‑150 (MQ=255)
           ccAGTTTTATCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCACGCGCCGCaaaata                                                                                                                                       <  1:665448/150‑6 (MQ=255)
              gTTTTATCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGttt                                                                                                                                    <  2:913594/150‑1 (MQ=255)
               ttttATCTTTTGGTAAGCATAATGCCATTTTTGCCCCCAAAAGACCATGAATGTTGCACATTAAAACATAACATCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCTCGTGCCGCATCGCGTTtg                                                                                                                                   <  1:1536659/150‑1 (MQ=255)
               ttttATCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTtg                                                                                                                                   <  2:982654/150‑1 (MQ=255)
               ttttATCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTtg                                                                                                                                   <  1:1542419/150‑1 (MQ=255)
                 ttATCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTtgtg                                                                                                                                 <  1:569932/150‑1 (MQ=255)
                   aTCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGAc                                                                                                                               >  1:108329/1‑150 (MQ=255)
                    tCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACg                                                                                                                              >  1:206928/1‑150 (MQ=255)
                     cTTTGGGGACGCATAATGCCATTTTTGGCCCGAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAATAGTGGTCGTACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACgg                                                                                                                             <  2:156047/150‑1 (MQ=255)
                     cTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACgg                                                                                                                             >  2:856419/1‑150 (MQ=255)
                     cTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACgg                                                                                                                             >  1:1271573/1‑150 (MQ=255)
                         ggggACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGatta                                                                                                                         >  2:1523669/1‑150 (MQ=255)
                            gACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGGCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTAttt                                                                                                                      <  1:1350848/150‑1 (MQ=255)
                                cATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTTGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGAt                                                                                                                  >  1:1477535/1‑150 (MQ=255)
                                     tGCCATTTTTGCCCCCACCAGACCATGCATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAACAACGTGGTCGCACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGGGACGGAGTATTTCGATGCCGc                                                                                                             >  2:671332/1‑150 (MQ=255)
                                     tGCCATTTTTGCCCCCAACCGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACGTTGGATAAAAAAGGGGTCGACCCGCGGAGGTACTTTTTATTTTGCGGCGCGTGCCGCCGCGCGGTTGTGCCGGATTATTTCGATGCCGc                                                                                                             >  2:666025/1‑150 (MQ=255)
                                      gCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAATTGGTCGAACCCCGAATTTAATTTTAATTTTGCGCCGCGTCCCCCACCCCTTTTGTGCCGTATTATTACGATGACGcc                                                                                                            >  2:310326/1‑150 (MQ=255)
                                         aTTTTTGCCCCCAACAGACCATGAATGTTGCACATTTACAGATAACAGCCCGAAACTTTGTATAAAAAAGTGGTCGAACCGCGGAGTTTCTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCgg                                                                                                         <  1:1098741/150‑1 (MQ=255)
                                         aTTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCgg                                                                                                         <  1:320439/150‑1 (MQ=255)
                                          tttttGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGc                                                                                                        <  2:489342/150‑1 (MQ=255)
                                              tGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCaaaa                                                                                                    <  1:1195907/150‑1 (MQ=255)
                                               gCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAAt                                                                                                   <  1:780075/150‑1 (MQ=255)
                                                cccccAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATa                                                                                                  >  1:168917/1‑150 (MQ=255)
                                                  cccAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGa                                                                                                >  2:562533/1‑150 (MQ=255)
                                                      aCAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAAcc                                                                                            <  1:1523669/150‑1 (MQ=255)
                                                         gACCATGAATGTTGCACATTAAAACATAACAGCCCGAAGCTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACcaca                                                                                         >  1:992246/1‑150 (MQ=255)
                                                                aaTGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATcc                                                                                  >  2:815331/1‑150 (MQ=255)
                                                                 aTGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATccc                                                                                 <  1:1129275/150‑1 (MQ=255)
                                                                   gTTGCACATTAAAACATAACAGACAGAAACTTTGGATAAAAAAGTGGTCGAACAGAGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCcga                                                                               >  2:485123/1‑150 (MQ=255)
                                                                    ttGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCcgac                                                                              <  1:343635/150‑1 (MQ=255)
                                                                       cacaTTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGAt                                                                           >  1:239970/1‑150 (MQ=255)
                                                                        acaTTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATc                                                                          <  2:1271573/150‑1 (MQ=255)
                                                                            tAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATcagca                                                                      <  2:1521842/150‑1 (MQ=255)
                                                                             aaaaCATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATcagcag                                                                     >  1:149809/1‑150 (MQ=255)
                                                                                 cATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTa                                                                 >  2:1353847/1‑150 (MQ=255)
                                                                                      cAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGtt                                                            >  1:1471063/1‑150 (MQ=255)
                                                                                        gCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTAt                                                          >  1:1167104/1‑150 (MQ=255)
                                                                                         cccGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGtatta                                                         >  1:1604629/1‑146 (MQ=255)
                                                                                          ccGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATcc                                                        >  1:1499225/1‑150 (MQ=255)
                                                                                            gAAACTTTTGATAAAAAAGTGGTCTAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCTTACGATCAGTAGCTTAAGGTTATCCTg                                                      <  2:693876/150‑1 (MQ=255)
                                                                                            gAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTg                                                      >  1:1247906/1‑150 (MQ=255)
                                                                                               aCTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAAc                                                   <  1:598244/150‑1 (MQ=255)
                                                                                                cTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAAcc                                                  >  2:217820/1‑150 (MQ=255)
                                                                                                cTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAAcc                                                  >  1:150759/1‑150 (MQ=255)
                                                                                                   tGGCTAAAAAAGGGGGCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCTTTTGTTACGGATTATTTCGATGACGCCCGGCAAAATAGAAAACAGAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATc                                               <  2:1465893/150‑1 (MQ=255)
                                                                                                         aaaaaGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGAtt                                         >  2:885351/1‑150 (MQ=255)
                                                                                                         aaaaaGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCAAGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGAtt                                         >  2:417856/1‑150 (MQ=255)
                                                                                                           aaaGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGt                                       <  2:1233156/150‑1 (MQ=255)
                                                                                                            aaGTGGTCGAACCGCGGAGTTACTTTTAATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGCGATTGTa                                      >  2:1138621/1‑150 (MQ=255)
                                                                                                              gTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTAcc                                    <  2:1076315/150‑1 (MQ=255)
                                                                                                              gTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTAcc                                    >  1:1599469/1‑150 (MQ=255)
                                                                                                                ggTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCTCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGCTCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCga                                  >  2:1530522/1‑150 (MQ=255)
                                                                                                                ggTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCga                                  >  2:19152/1‑150 (MQ=255)
                                                                                                                 gTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCgaa                                 <  1:519702/150‑1 (MQ=255)
                                                                                                                 gTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCgaa                                 >  1:313656/1‑150 (MQ=255)
                                                                                                                  tCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCgaag                                <  1:598959/150‑1 (MQ=255)
                                                                                                                  tCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCgaag                                <  2:181603/150‑1 (MQ=255)
                                                                                                                  tCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCgaag                                <  2:1599469/150‑1 (MQ=255)
                                                                                                                      aCCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAAt                            <  1:1284443/150‑1 (MQ=255)
                                                                                                                       ccGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATa                           >  2:1565632/1‑150 (MQ=255)
                                                                                                                         gcgGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATTGGGATTGTACCGAAGAAATAGc                         >  2:1120096/1‑150 (MQ=255)
                                                                                                                         gcgGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGc                         >  2:1384818/1‑150 (MQ=255)
                                                                                                                         gcgGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGc                         >  2:1531719/1‑150 (MQ=255)
                                                                                                                         gcgGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGc                         <  1:249232/150‑1 (MQ=255)
                                                                                                                            gAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTg                      <  2:1494157/150‑1 (MQ=255)
                                                                                                                                tACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAgg                                                              <  1:25604/106‑1 (MQ=255)
                                                                                                                                tACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAgg                                                              >  2:25604/1‑106 (MQ=255)
                                                                                                                                tACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTa                  <  1:1557835/150‑1 (MQ=255)
                                                                                                                                tACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTa                  <  1:1353847/150‑1 (MQ=255)
                                                                                                                                tACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTa                  >  1:1345756/1‑150 (MQ=255)
                                                                                                                                   tttttATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGt               <  1:310326/150‑1 (MQ=255)
                                                                                                                                    ttttATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTa              <  1:815331/150‑1 (MQ=255)
                                                                                                                                    ttttATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTa              <  1:1116934/150‑1 (MQ=255)
                                                                                                                                       tATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTaaaa           >  1:24997/1‑150 (MQ=255)
                                                                                                                                       cgcgctGCGGCACGCGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTaaaa           >  2:194605/6‑150 (MQ=255)
                                                                                                                                        gcgctGCGGCACGCGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTaaaaa          >  1:642625/5‑150 (MQ=255)
                                                                                                                                        aTTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTaaaaa          >  2:608157/1‑150 (MQ=255)
                                                                                                                                        aTTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTaaaaa          >  1:207055/1‑150 (MQ=255)
                                                                                                                                        aTTTTGCGGCCCGTGCCGAACCCCTTTTGTGACGGATTATTTAGATGACGCCCGGCAAAATAGAAACCACAATAATCCAGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTaaaaa          >  2:540693/1‑150 (MQ=255)
                                                                                                                                         ttttGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTAAAAAg         <  2:25155/150‑1 (MQ=255)
                                                                                                                                         cgctGCGGCACGCGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTaaaacg         >  2:1474805/4‑148 (MQ=255)
                                                                                                                                         cgctGCGGCACGCGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTAAAAAg         >  2:207287/4‑150 (MQ=255)
                                                                                                                                         cgctGCGGCACGCGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTAAAAAg         >  2:1074463/4‑150 (MQ=255)
                                                                                                                                         cgctGCGGCACGCGCCGCAGCGCGTTTGTGACGGACTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTAAAAAg         >  1:1467346/4‑150 (MQ=255)
                                                                                                                                            tGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCTATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTAAAAAGCAg      <  1:666025/150‑1 (MQ=255)
                                                                                                                                            tGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTCAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTAAAAAGCAg      <  1:671332/150‑1 (MQ=255)
                                                                                                                                            tGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTAAAAAGCAg      <  1:19152/150‑1 (MQ=255)
                                                                                                                                                gcgcGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTAAAAAGCAGTAcc  <  2:108329/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTGGTTCGATGCCAGTTTTATCTTTGGGGACGCATAATGCCATTTTTGCCCCCAACAGACCATGAATGTTGCACATTAAAACATAACAGCCCGAAACTTTGGATAAAAAAGTGGTCGAACCGCGGAGTTACTTTTTATTTTGCGGCGCGTGCCGCAGCGCGTTTGTGACGGATTATTTCGATGACGCCCGGCAAAATAGAAACCACAATAATCCCGACGATCAGCAGCTTAAGGTTATCCTGAACCATCGGGATTGTACCGAAGAAATAGCCTGCGTAGGTAAAAAGCAGTACC  >  NC_000913/2433948‑2434241

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