Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,602,080 T→C K3E (AAG→GAG)  ftsX ← putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,602,0800TC97.3% 266.9 / ‑1.4 75K3E (AAG→GAG) ftsXputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (45/28);  total (47/28)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

AAACGTTGGTTTTGCGATTTACCGGTTTTGGCGAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTTATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCGGTTAAACTCTTCAAACAGACGTAAAATGCCTT  >  NC_000913/3601933‑3602227
                                                                                                                                                   |                                                                                                                                                   
aaaCGTTGGTTTTGCGATTTACCGGTTTTGGCGAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCAtt                                                                                                                                                  >  2:1188935/1‑151 (MQ=255)
  aCGTTGGTTTTGCGATTTACCGGTTTTGGCGAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTc                                                                                                                                                 >  1:1058153/1‑150 (MQ=255)
   cGTTGGTTTTGCGATTTACCGGTTTTGGCGAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCa                                                                                                                                                >  1:419637/1‑150 (MQ=255)
    gTTGGTTTTGCGATTTACCGGTTTTGGCGAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATg                                                                                                                                              <  2:883943/151‑1 (MQ=255)
     ttGTTTTTGCGATTTACCGGTTTTGGCGAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCtcattcat                                                                                                                                               <  2:138149/149‑1 (MQ=255)
        gTTTTGCGATTTACCGGTTTTGGCGAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGccc                                                                                                                                          <  2:889793/151‑1 (MQ=255)
         ttttGCGATTTACCGGTTTTGGCGAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGccc                                                                                                                                          >  1:1055943/1‑150 (MQ=255)
         ttttGCGATTTACCGGTTTTGGCGAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGccc                                                                                                                                          <  1:939404/150‑1 (MQ=255)
         ttttGCGATTTACCGGTTTTGGCGAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGccc                                                                                                                                          >  1:13115/1‑150 (MQ=255)
           ttGCGATTTACCGGTTTTGGCGAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGACCACg                                                                                                                                       >  2:1130716/1‑151 (MQ=255)
             gCGATTTACCGGTTTTGGCGAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACg                                                                                                                                       <  1:1049499/149‑1 (MQ=255)
              cGATTTACCGGTTTTGGCGAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCt                                                                                                                                    >  2:421366/1‑151 (MQ=255)
                 tttACCGGTTATGGCGAGGATTACGCGCGTTTTGGAGCGTTACGACAGCCGTAGCCTGAGCCGCCGACCGATTTACGGAAGCGAACAAAACGCACGCCAAACTACCGAATCTGATTCATTGCATCGCGCTCATTCATGGCCCACGCcccc                                                                                                                                  >  1:699641/1‑150 (MQ=255)
                       ggTTTTGGCGAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAg                                                                                                                           <  2:508412/151‑1 (MQ=255)
                          tttGGCGAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTg                                                                                                                         >  1:858808/1‑150 (MQ=255)
                          tttGGCGAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTg                                                                                                                         >  2:544341/1‑150 (MQ=255)
                            tGGCGAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGAcc                                                                                                                      >  2:916033/1‑151 (MQ=255)
                                gAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCg                                                                                                                  >  2:381872/1‑151 (MQ=255)
                                  ggATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCt                                                                                                                >  1:1176008/1‑151 (MQ=255)
                                       tCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAggg                                                                                                           >  1:431565/1‑151 (MQ=255)
                                       tCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAgg                                                                                                            <  1:621417/150‑1 (MQ=255)
                                            cgTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGc                                                                                                      <  2:26532/151‑1 (MQ=255)
                                            cgTATTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGc                                                                                                      <  2:293327/151‑1 (MQ=255)
                                             gTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGc                                                                                                      >  2:1265683/1‑150 (MQ=255)
                                               tttGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCAt                                                                                                    <  2:1074499/150‑1 (MQ=255)
                                                 tGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGc                                                                                                  <  1:525938/150‑1 (MQ=255)
                                                 tGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGc                                                                                                  <  2:43349/150‑1 (MQ=255)
                                                  ggTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGa                                                                                                <  2:1096175/151‑1 (MQ=255)
                                                     gCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGAt                                                                                               >  2:507125/1‑149 (MQ=255)
                                                     gCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGAt                                                                                               >  1:585463/1‑149 (MQ=255)
                                                      cGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAgg                                                                                            <  2:504971/151‑1 (MQ=255)
                                                         tACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAAc                                                                                         >  1:802372/1‑151 (MQ=255)
                                                          aCGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAAcg                                                                                        >  2:861683/1‑151 (MQ=255)
                                                             aCCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAAcgccg                                                                                     >  2:581844/1‑151 (MQ=255)
                                                               cgccgTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAAcgccg                                                                                     <  1:444130/149‑1 (MQ=255)
                                                               cgccgTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCcgcg                                                                                   <  1:456193/151‑1 (MQ=255)
                                                                gccgTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCga                                                                                  <  1:21645/151‑1 (MQ=255)
                                                                     cGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATc                                                                              >  2:678337/1‑150 (MQ=255)
                                                                      gCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATc                                                                              >  2:1223306/1‑149 (MQ=255)
                                                                        cTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGt                                                                          >  1:955957/1‑151 (MQ=255)
                                                                        cTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTAACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCa                                                                             >  1:71446/1‑148 (MQ=255)
                                                                         tGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGtt                                                                         >  1:472550/1‑151 (MQ=255)
                                                                          gAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCGCGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTg                                                                        <  1:1030638/151‑1 (MQ=255)
                                                                          gAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGt                                                                          <  1:1188901/149‑1 (MQ=255)
                                                                            gccgccGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGa                                                                       >  1:483051/1‑150 (MQ=255)
                                                                                ccgaccgaTTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCg                                                                  >  1:459394/1‑151 (MQ=255)
                                                                                        tttACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTTATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGc                                                          >  2:883029/1‑151 (MQ=255)
                                                                                          tACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTTATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGcc                                                         >  2:1210098/1‑150 (MQ=255)
                                                                                           aCGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCAt                                                       >  2:279569/1‑151 (MQ=255)
                                                                                                aGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAAt                                                   >  2:972035/1‑150 (MQ=255)
                                                                                                aGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATa                                                  <  2:1101325/151‑1 (MQ=255)
                                                                                                 gCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATAc                                                 <  2:128045/151‑1 (MQ=255)
                                                                                                   gATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACgg                                               <  1:153119/151‑1 (MQ=255)
                                                                                                   gATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACgg                                               >  1:993576/1‑151 (MQ=255)
                                                                                                     tCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGtt                                             >  1:1292532/1‑151 (MQ=255)
                                                                                                     tCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGtt                                             >  1:310597/1‑151 (MQ=255)
                                                                                                            gcccgccAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCa                                       >  1:1267917/1‑150 (MQ=255)
                                                                                                            gcccgccAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCa                                       >  1:775531/1‑150 (MQ=255)
                                                                                                             cccgccAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAAc                                     <  2:127291/151‑1 (MQ=255)
                                                                                                                 ccAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCgg                                 <  2:676376/151‑1 (MQ=255)
                                                                                                                  cAAACTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCGGt                                >  2:1068387/1‑151 (MQ=255)
                                                                                                                    aaCTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCGGTTa                              <  1:544314/151‑1 (MQ=255)
                                                                                                                      cTGCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCGGtt                               >  2:730603/1‑148 (MQ=255)
                                                                                                                        gCCGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCGGTTAAAct                          >  1:423425/1‑151 (MQ=255)
                                                                                                                         ccGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCGGTTAAActc                         >  2:1064952/1‑151 (MQ=255)
                                                                                                                          cGAATATGATTGATTGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCGGTTAAActct                        >  2:622773/1‑151 (MQ=255)
                                                                                                                                 gattgattGCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCGGTTAAACTCTTCAAACa                 <  2:429106/151‑1 (MQ=255)
                                                                                                                                    tgattgCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCGGTTAAACTCTTCAAACAGa               <  2:946694/150‑1 (MQ=255)
                                                                                                                                     gattgCATCGCGCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCGGTTAAACTCTTCAAACAGACg             <  2:943850/151‑1 (MQ=255)
                                                                                                                                            tCGCCCTCATTCATTGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCGGTTAAACTCTTCAAACAGACGTAAAAt       >  1:987089/1‑150 (MQ=255)
                                                                                                                                             cgcgCTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCGGTTAAACTCTTCAAACAGACGTAAAATGc     >  1:6022/1‑151 (MQ=255)
                                                                                                                                              gcgcTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCGGTTAAACTCTTCAAACAGACGTAAAATg      >  1:450780/1‑149 (MQ=255)
                                                                                                                                              gcgcTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCGGTTAAACTCTTCAAACAGACGTAAAATGcc    >  1:377590/1‑151 (MQ=255)
                                                                                                                                               cgcTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCGGTTAAACTCTTCAAACAGACGTAAAATGCCt   >  2:556695/1‑151 (MQ=255)
                                                                                                                                               cgcTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCGGTTAAACTCTTCAAACAGACGTAAAATGCCt   >  1:815330/1‑151 (MQ=255)
                                                                                                                                                gcTCATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCGGTTAAACTCTTCAAACAGACGTAAAATGCCtt  <  1:1149581/151‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AAACGTTGGTTTTGCGATTTACCGGTTTTGGCGAGGATTTCGCGCGTTTTGGTGCGTTACGACCGCCGTCGCCTGAGCCGCCGACCGATTTACGGAAGCGATCAAGACGCCCGCCAAACTGCCGAATATGATTGATTGCATCGCGCTTATTCATGGCCCACGCCTCCATGCAAGTGACCATCGCTCAGGGTGAGCATGCGATAGGAACGCCGCGAGATCAGGTTGATGTCGTGCGTTGCCATCAATACGGTTACCCCAACGCGGTTAAACTCTTCAAACAGACGTAAAATGCCTT  >  NC_000913/3601933‑3602227

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