Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,352,163 G→A R158C (CGT→TGT)  sapF ← antimicrobial peptide ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,352,1630GA98.8% 319.2 / NA 86R158C (CGT→TGT) sapFantimicrobial peptide ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (50/35);  total (50/36)
Fisher's exact test for biased strand distribution p-value = 4.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATAAATATACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATT  >  NC_000913/1352017‑1352308
                                                                                                                                                  |                                                                                                                                                 
aTAAATATACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGt                                                                                                                                               >  1:1094611/1‑151 (MQ=255)
aTAAATATACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGt                                                                                                                                               >  1:332964/1‑151 (MQ=255)
  aaaTATACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGttt                                                                                                                                             >  1:332787/1‑151 (MQ=255)
   aaTATACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGtttt                                                                                                                                            <  1:821108/151‑1 (MQ=255)
   aaTATACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGtttt                                                                                                                                            >  1:93141/1‑151 (MQ=255)
    atatACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATAGCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGttttt                                                                                                                                           <  2:940210/151‑1 (MQ=255)
     tataCGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCTCCAGACCCAGACACTGttttt                                                                                                                                           <  2:844179/150‑1 (MQ=255)
      ataCGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTg                                                                                                                                          <  2:117921/150‑1 (MQ=255)
       taCGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTc                                                                                                                                        >  1:780888/1‑151 (MQ=255)
         cGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCg                                                                                                                                      >  2:18041/1‑151 (MQ=255)
          gAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCgg                                                                                                                                     <  1:1133803/151‑1 (MQ=255)
            aaTGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGt                                                                                                                                   >  2:112759/1‑151 (MQ=255)
             aTGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTg                                                                                                                                  <  1:483597/151‑1 (MQ=255)
               gCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCGTCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGcc                                                                                                                                >  1:149936/1‑151 (MQ=255)
               gCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGc                                                                                                                                 <  2:1112419/150‑1 (MQ=255)
                gccTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCTGTGAAGCGAGGACTTCATTGGCAATAATAACTTATGAGCCCAATATCAGCGCACGCGTCCGACCCACACACTGTTTTTACCAGGGTg                                                                                                                                  >  2:494471/2‑148 (MQ=255)
                 ccTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCaa                                                                                                                              <  2:436685/151‑1 (MQ=255)
                   tGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAAc                                                                                                                             >  2:1137495/1‑150 (MQ=255)
                    gTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACa                                                                                                                            <  1:350309/150‑1 (MQ=255)
                      ttttCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACatat                                                                                                                         >  2:749366/1‑151 (MQ=255)
                       tttCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATg                                                                                                                        >  2:980384/1‑151 (MQ=255)
                       tttCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATg                                                                                                                        >  1:244678/1‑151 (MQ=255)
                          ccTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCgg                                                                                                                     >  1:770291/1‑151 (MQ=255)
                           cTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCgg                                                                                                                     >  2:1067558/1‑150 (MQ=255)
                           cTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCgg                                                                                                                     >  2:1410144/1‑150 (MQ=255)
                             gTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCTCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATa                                                                                                                  <  2:1196858/151‑1 (MQ=255)
                             gTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGa                                                                                                                    <  1:1221043/149‑1 (MQ=255)
                             gTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATa                                                                                                                  <  1:140917/151‑1 (MQ=255)
                               aaTTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGt                                                                                                                >  1:668550/1‑151 (MQ=255)
                                aTTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTa                                                                                                               <  1:928559/151‑1 (MQ=255)
                                  tCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAAc                                                                                                             >  2:500908/1‑151 (MQ=255)
                                   ccAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAAc                                                                                                             >  2:182238/1‑150 (MQ=255)
                                    cAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTg                                                                                                           >  1:1398971/1‑151 (MQ=255)
                                        aTCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGt                                                                                                       >  2:742428/1‑151 (MQ=255)
                                        aTCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGt                                                                                                       >  1:343197/1‑151 (MQ=255)
                                          cAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTg                                                                                                      <  2:1292746/150‑1 (MQ=255)
                                          cAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTg                                                                                                      >  2:842433/1‑150 (MQ=255)
                                             aTTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATcc                                                                                                  >  1:73292/1‑151 (MQ=255)
                                               tAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCgg                                                                                                >  1:452724/1‑151 (MQ=255)
                                                 aTCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGt                                                                                               >  2:697717/1‑150 (MQ=255)
                                                 aTCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGt                                                                                               >  1:109419/1‑150 (MQ=255)
                                                   cAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTagc                                                                                            <  1:74327/151‑1 (MQ=255)
                                                     aCTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTagc                                                                                            >  2:655141/1‑149 (MQ=255)
                                                       tGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGcc                                                                                        >  2:808206/1‑151 (MQ=255)
                                                        gCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTagcagc                                                                                         <  2:1390854/149‑1 (MQ=255)
                                                        gCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGccc                                                                                       <  1:902677/151‑1 (MQ=255)
                                                         cGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCa                                                                                      <  1:588682/151‑1 (MQ=255)
                                                         cGAACGCATCGACATATCCAGTGAAGCGAGGGCTGCATCGGCAATAATGACTTTTGGGCGCACTATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGCGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGcc                                                                                        >  2:837282/1‑149 (MQ=255)
                                                           aaCGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAAc                                                                                    >  1:1380164/1‑151 (MQ=255)
                                                             cGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAAcc                                                                                   >  1:814625/1‑150 (MQ=255)
                                                               cATCGCCATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCTGTAGCAGCCCAACCAt                                                                                 >  2:89896/1‑150 (MQ=255)
                                                                atagaaatataaaGTGAAGCGAGGGATTCATAGGAAATAATGACTTTTGGGCGCAATATAAGAGTGTGAGCCAGAGTCAGAGACTGTTTTTGTAAGGGTGAAAACATATGAGGATAGTAAATGACGTGATAAGGTAGCAGCGCAACAATAc                                                                               >  2:441398/13‑151 (MQ=255)
                                                                   gACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCa                                                                            <  2:543372/151‑1 (MQ=255)
                                                                    aCATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCAt                                                                           <  1:325620/151‑1 (MQ=255)
                                                                    aCATATCCACGGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCCAGCGCCATACCCAGACACTGTTTTTGTCCCGGTGCCAACATATGCGGATAGTACCTGACGTGATCCGGTAGCAGACCAACCATACGCAt                                                                           <  2:931654/151‑1 (MQ=255)
                                                                     cATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATc                                                                          <  1:259879/151‑1 (MQ=255)
                                                                          ccAGTGAAGCGAGGGCCTCAGCGGCAAGAATGATGTTTGGGCGCAATATCAGCGCGCGTGCTAGAGCCAGACGCTGTTTTTGTCCGGGTGGCAACAGATGCGGATTGTAACTGACGTGATCCGGTAGCAGCCCAACCATTCGGGTCGtt                                                                       <  2:1060278/149‑1 (MQ=255)
                                                                            aGTGAAGCCAGGGCTTCATCGGCAATAATGACTTTTGGGGGCCAAATCAGGGCGCGCGCCAGACCCAGAAACTGTTTTTGTCCGGGTGCCAACATATGCGGAAAGTAACAGACGCGATTCGGTAGGAGCCCAACCATACGCacgggttca                                                                    >  2:426984/1‑142 (MQ=255)
                                                                             gTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCaat                                                                  <  1:112780/151‑1 (MQ=255)
                                                                                    gAGGGCTTCATCGGCAATACTGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATGTGCGGATAGTAACTGATGTGATCCGGTAGCAGCACAACCATACGCATCGTTTCAATAATCTg                                                            >  2:957115/1‑150 (MQ=255)
                                                                                    gAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGt                                                           >  1:1181418/1‑151 (MQ=255)
                                                                                      gggCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGttt                                                         <  1:71500/151‑1 (MQ=255)
                                                                                          ttCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCACTAATCTGTTTGCg                                                      >  2:1378674/1‑150 (MQ=255)
                                                                                            cATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACg                                                   >  1:301976/1‑151 (MQ=255)
                                                                                            cATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACg                                                   <  2:458357/151‑1 (MQ=255)
                                                                                               cGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCt                                                 <  1:875578/150‑1 (MQ=255)
                                                                                                 gCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGtt                                              <  2:332882/151‑1 (MQ=255)
                                                                                                      aatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGt                                         >  2:284632/1‑151 (MQ=255)
                                                                                                         gACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTCAATAATCTGTTTGCGACGCTGTTCCGGTTcc                                      >  1:1268661/1‑150 (MQ=255)
                                                                                                             tttGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAgg                                   <  1:1232038/150‑1 (MQ=255)
                                                                                                             tttGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGCGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAgg                                   >  2:1404984/1‑150 (MQ=255)
                                                                                                               tGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTACCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTc                                 <  2:1038627/150‑1 (MQ=255)
                                                                                                                  gcgcAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACACTGTTCCGGTTCCAGGTCGgtg                             >  2:541448/1‑151 (MQ=255)
                                                                                                                       atatCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAg                        <  2:243100/151‑1 (MQ=255)
                                                                                                                              cgcgcgcgcCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGt                 >  1:51700/1‑151 (MQ=255)
                                                                                                                                cgcgcgcCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTgg               >  2:942771/1‑151 (MQ=255)
                                                                                                                                 gcgcgcCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTgg               >  1:343945/1‑150 (MQ=255)
                                                                                                                                      ccagacccagacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATc         <  1:547724/151‑1 (MQ=255)
                                                                                                                                      ccagacccagacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATc         >  1:1168654/1‑151 (MQ=255)
                                                                                                                                       cagacccagacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATcc        <  1:1410564/151‑1 (MQ=255)
                                                                                                                                        agacccagacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATcc        >  2:1224984/1‑150 (MQ=255)
                                                                                                                                        agacccagacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATca        <  2:1262550/150‑2 (MQ=255)
                                                                                                                                          acccagacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATcc        >  1:612204/1‑148 (MQ=255)
                                                                                                                                          acccagacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATcc        >  1:301954/1‑148 (MQ=255)
                                                                                                                                             cagacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGTTCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAg      <  2:218629/147‑1 (MQ=255)
                                                                                                                                             cagacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAAtt  <  1:873545/151‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
ATAAATATACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATT  >  NC_000913/1352017‑1352308

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