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,1630GA96.6% 311.6 / ‑1.3 88R158C (CGT→TGT) sapFantimicrobial peptide ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (2/1);  new base A (50/35);  total (52/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

ACATAAATATACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTT  >  NC_000913/1352015‑1352309
                                                                                                                                                    |                                                                                                                                                  
aCATAAATATACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCt                                                                                                                                                  >  1:78183/1‑151 (MQ=255)
   tAAATATACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGtt                                                                                                                                               >  1:44098/1‑151 (MQ=255)
    aaaTATACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGttt                                                                                                                                              <  2:891081/151‑1 (MQ=255)
     aaTATACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGtttt                                                                                                                                             <  1:568069/151‑1 (MQ=255)
     aaTATACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGtttt                                                                                                                                             >  2:811638/1‑151 (MQ=255)
     aaTATACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGtttt                                                                                                                                             >  1:1074806/1‑151 (MQ=255)
       tataCGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTg                                                                                                                                           >  1:1054500/1‑151 (MQ=255)
        ataCGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGt                                                                                                                                          >  2:916222/1‑151 (MQ=255)
        ataCGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGt                                                                                                                                          >  1:1037451/1‑151 (MQ=255)
         taCGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTc                                                                                                                                         >  2:574101/1‑151 (MQ=255)
          aCGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTcc                                                                                                                                        >  1:196077/1‑151 (MQ=255)
            gAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCgg                                                                                                                                      >  2:825522/1‑151 (MQ=255)
            gAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCgg                                                                                                                                      >  2:772263/1‑151 (MQ=255)
            gAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCgg                                                                                                                                      >  2:426627/1‑151 (MQ=255)
                      gTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAAc                                                                                                                              <  2:189647/149‑1 (MQ=255)
                      gTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACat                                                                                                                            <  1:31728/151‑1 (MQ=255)
                      gTTTTTCATGTAAATCCAGCATCTGATTAATCAACTTCGCACGCATCGTCATATCCAGTGAAGCGAGGGCTTAATCGGCAATAATGACCCATGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCa                                                                                                                                >  1:393617/1‑147 (MQ=255)
                        ttttCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACata                                                                                                                           <  2:1106555/150‑1 (MQ=255)
                            ccTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCgg                                                                                                                      >  1:821719/1‑151 (MQ=255)
                            ccTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCgg                                                                                                                      <  2:887968/151‑1 (MQ=255)
                               gTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGAt                                                                                                                    <  2:135811/150‑1 (MQ=255)
                               gTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATa                                                                                                                   <  2:249345/151‑1 (MQ=255)
                                  aTTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGt                                                                                                                 >  1:441714/1‑150 (MQ=255)
                                  aTTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGt                                                                                                                 >  2:956027/1‑150 (MQ=255)
                                  aTTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTa                                                                                                                <  1:604573/151‑1 (MQ=255)
                                      cAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTg                                                                                                            >  1:420667/1‑151 (MQ=255)
                                               aTTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATcc                                                                                                   >  2:775645/1‑151 (MQ=255)
                                                 tAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCg                                                                                                  <  2:1025173/150‑1 (MQ=255)
                                                 tAATCAACTACGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCgg                                                                                                 >  2:617497/1‑151 (MQ=255)
                                                  aaTCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCgg                                                                                                 <  2:128374/150‑1 (MQ=255)
                                                  aaTCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCgg                                                                                                 <  2:1105157/150‑1 (MQ=255)
                                                  aaTCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGt                                                                                                >  2:1015065/1‑151 (MQ=255)
                                                   aTCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTa                                                                                               <  2:927304/151‑1 (MQ=255)
                                                     cAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTa                                                                                               <  1:572251/149‑1 (MQ=255)
                                                     cAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTag                                                                                              <  2:989290/150‑1 (MQ=255)
                                                        cTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTagcagc                                                                                          >  2:461593/1‑151 (MQ=255)
                                                         tGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGcc                                                                                         >  1:1023749/1‑151 (MQ=255)
                                                             aaCGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAAc                                                                                     >  2:624299/1‑151 (MQ=255)
                                                               cGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAAcc                                                                                    >  2:350719/1‑150 (MQ=255)
                                                                 cATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCAt                                                                                  >  1:718602/1‑150 (MQ=255)
                                                                  aTCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATAc                                                                                >  1:101738/1‑151 (MQ=255)
                                                                     gACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGc                                                                              >  2:813197/1‑150 (MQ=255)
                                                                       cATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCAt                                                                            <  1:87876/150‑1 (MQ=255)
                                                                        atatCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCg                                                                          >  2:63626/1‑151 (MQ=255)
                                                                        atatCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCg                                                                          <  2:142370/151‑1 (MQ=255)
                                                                          atCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGtt                                                                        >  1:762214/1‑151 (MQ=255)
                                                                                  aaGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCaataat                                                                <  2:312158/151‑1 (MQ=255)
                                                                                        gggCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGttt                                                          <  2:1054635/151‑1 (MQ=255)
                                                                                            ttCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCg                                                       >  2:824958/1‑150 (MQ=255)
                                                                                              cATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGAc                                                     <  2:122486/150‑1 (MQ=255)
                                                                                               aTCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGc                                                   >  2:315833/1‑151 (MQ=255)
                                                                                                 cGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTg                                                 >  2:159961/1‑151 (MQ=255)
                                                                                                 cGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTg                                                 >  2:770658/1‑151 (MQ=255)
                                                                                                   gCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGtt                                               >  2:484941/1‑151 (MQ=255)
                                                                                                    cAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTc                                              >  1:248078/1‑151 (MQ=255)
                                                                                                      ataatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCg                                            >  2:826005/1‑151 (MQ=255)
                                                                                                       taatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCgg                                           >  1:417482/1‑151 (MQ=255)
                                                                                                       taatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCgg                                           >  1:934412/1‑151 (MQ=255)
                                                                                                        aatGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCgg                                           <  2:544532/150‑1 (MQ=255)
                                                                                                         atGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGtt                                         >  2:64037/1‑151 (MQ=255)
                                                                                                           gACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTc                                        <  2:111006/150‑1 (MQ=255)
                                                                                                            aCTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCa                                      <  2:5841/151‑1 (MQ=255)
                                                                                                             cTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTcc                                       >  1:663258/1‑149 (MQ=255)
                                                                                                             cTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTcc                                       >  1:84173/1‑149 (MQ=255)
                                                                                                              ttttGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAgg                                    >  2:880648/1‑151 (MQ=255)
                                                                                                                 tGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTc                                  <  2:854019/150‑1 (MQ=255)
                                                                                                                   ggCGCAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCgg                                >  1:538189/1‑150 (MQ=255)
                                                                                                                    gcgcAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGgtg                              >  2:370014/1‑151 (MQ=255)
                                                                                                                      gcAATATCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGtt                            <  2:315721/151‑1 (MQ=255)
                                                                                                                         atatCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAg                         <  1:1128603/151‑1 (MQ=255)
                                                                                                                           atCAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAgg                        <  1:577713/150‑1 (MQ=255)
                                                                                                                             cAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGgcgc                     <  2:685695/151‑1 (MQ=255)
                                                                                                                             cAGCGCGCGCGCCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGgcg                      <  1:881553/150‑1 (MQ=255)
                                                                                                                                cgcgcgcgcCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGt                  <  2:647151/151‑1 (MQ=255)
                                                                                                                                cgcgcgcgcCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGt                  <  2:968831/151‑1 (MQ=255)
                                                                                                                                 gcgcgcgcCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTg                 <  1:402609/151‑1 (MQ=255)
                                                                                                                                 gcgcgcgcCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTg                 >  1:465460/1‑151 (MQ=255)
                                                                                                                                  cgcgcgcCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTg                 <  2:288818/150‑1 (MQ=255)
                                                                                                                                  cgcgcgcCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTg                 <  1:1025415/150‑1 (MQ=255)
                                                                                                                                    cgcgcCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGaa              <  2:825519/151‑1 (MQ=255)
                                                                                                                                       gcCAGACCCAGACACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAAt           >  1:29134/1‑151 (MQ=255)
                                                                                                                                         cagacccagacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATcc         >  2:764623/1‑151 (MQ=255)
                                                                                                                                           gacccagacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATcc         >  2:61671/1‑149 (MQ=255)
                                                                                                                                            acccagacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATcc         >  2:21341/1‑148 (MQ=255)
                                                                                                                                             cccagacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTGACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGa      >  1:119651/1‑150 (MQ=255)
                                                                                                                                             cccagacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGaa     <  2:651472/151‑1 (MQ=255)
                                                                                                                                              ccagacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAAt    >  1:14278/1‑151 (MQ=255)
                                                                                                                                                agacACTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAAttt  >  1:1061476/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACATAAATATACGAAATGCCCTGTTTTTCCTGTAATTCCAGCATCAGATTAATCAACTGCGAACGCATCGACATATCCAGTGAAGCGAGGGCTTCATCGGCAATAATGACTTTTGGGCGCAATATCAGCGCGCGCGCCAGACCCAGACGCTGTTTTTGTCCGGGTGCCAACATATGCGGATAGTAACTGACGTGATCCGGTAGCAGCCCAACCATACGCATCGTTTCAATAATCTGTTTGCGACGCTGTTCCGGTTCCAGGTCGGTGTTCAGGCGCAGTGGAAAATCCAGAATTT  >  NC_000913/1352015‑1352309

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