Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,354,687 C→A V262F (GTT→TTT)  sapB ← antimicrobial peptide transport ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,354,6870CA98.7% 285.3 / NA 75V262F (GTT→TTT) sapBantimicrobial peptide transport ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (35/39);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATA  >  NC_000913/1354540‑1354831
                                                                                                                                                   |                                                                                                                                                
ggTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACAt                                                                                                                                               <  1:641390/151‑1 (MQ=255)
    aGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCg                                                                                                                                           >  1:723782/1‑151 (MQ=255)
     gCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCgg                                                                                                                                          <  1:1029269/151‑1 (MQ=255)
     gCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCgg                                                                                                                                          <  1:159089/151‑1 (MQ=255)
       cATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTg                                                                                                                                        <  1:714700/151‑1 (MQ=255)
        aTGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTg                                                                                                                                        >  1:37353/1‑150 (MQ=255)
              cccAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCg                                                                                                                                 >  1:1070102/1‑151 (MQ=255)
                 aaaaTATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGcc                                                                                                                               >  2:57564/1‑150 (MQ=255)
                  aaaTATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGcc                                                                                                                               >  2:711254/1‑149 (MQ=255)
                   aaTATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCGAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCatc                                                                                                                            <  2:1275508/151‑3 (MQ=255)
                    atatCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCg                                                                                                                           >  1:191058/1‑151 (MQ=255)
                    atatCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCg                                                                                                                           >  2:1345732/1‑151 (MQ=255)
                           aaaTCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCAt                                                                                                                    >  2:306470/1‑151 (MQ=255)
                            aaTCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCAtg                                                                                                                   <  1:623107/151‑1 (MQ=255)
                             aTCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGa                                                                                                                                       >  1:1203226/1‑130 (MQ=255)
                             aTCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGa                                                                                                                                       <  2:1203201/130‑1 (MQ=255)
                             aTCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCAtgg                                                                                                                  >  1:749344/1‑151 (MQ=255)
                                   ttAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGaa                                                                                                             >  2:1343836/1‑150 (MQ=255)
                                    tAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGaaa                                                                                                            >  1:488097/1‑150 (MQ=255)
                                     aaCAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGGCAGCGTCAGCATGGTGGAAAAc                                                                                                          <  2:972905/151‑1 (MQ=255)
                                      aCAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACt                                                                                                         >  1:429555/1‑151 (MQ=255)
                                        aataatGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGt                                                                                                       >  1:52070/1‑151 (MQ=255)
                                        aataatGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGGGGAAAACTg                                                                                                        >  2:1070888/1‑150 (MQ=255)
                                         ataatGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGTTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTa                                                                                                      <  2:801583/151‑1 (MQ=255)
                                         ataatGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTa                                                                                                      <  2:500689/151‑1 (MQ=255)
                                           aatGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTaa                                                                                                     <  1:402760/150‑1 (MQ=255)
                                           aatGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAg                                                                                                    <  2:398259/151‑1 (MQ=255)
                                              gACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGccc                                                                                                 >  2:310846/1‑151 (MQ=255)
                                               aCCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGccc                                                                                                 >  2:1268084/1‑150 (MQ=255)
                                                ccAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAg                                                                                               <  1:1001929/151‑1 (MQ=255)
                                                  aGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGc                                                                                             >  1:629286/1‑151 (MQ=255)
                                                     gAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAg                                                                                          <  2:148154/151‑1 (MQ=255)
                                                         cacacaCCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAAt                                                                                      <  2:996807/151‑1 (MQ=255)
                                                             cacCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGgaagaaag                                                                                  <  1:142449/151‑3 (MQ=255)
                                                             cacCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAAc                                                                                   <  2:1016908/150‑1 (MQ=255)
                                                                cATCCCTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACcggc                                                                               <  2:1070095/151‑1 (MQ=255)
                                                                 aTCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACcggcg                                                                              >  1:70527/1‑151 (MQ=255)
                                                                        cggcggAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAAcgc                                                                       <  1:25831/151‑1 (MQ=255)
                                                                          gcggAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGt                                                                     <  1:637906/151‑1 (MQ=255)
                                                                           cggAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGtt                                                                    <  2:843750/151‑1 (MQ=255)
                                                                            ggAAATGGCTGCGTGCTCCTGCCGTCGGCTGGCGTTAAGTTACCAGGGCCCTAATTCCGGACAGCTAAAAACCATTACTGGGAGCCTCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAAcc                                                                                  <  2:1093748/136‑1 (MQ=255)
                                                                            ggAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGACAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGAGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGtt                                                                    >  1:138147/1‑150 (MQ=255)
                                                                              aaaTGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATg                                                                 >  2:590507/1‑151 (MQ=255)
                                                                              aaaTGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATg                                                                 >  2:233357/1‑151 (MQ=255)
                                                                                   gCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGaa                                                            <  1:1020728/151‑1 (MQ=255)
                                                                                     tGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACg                                                          >  1:513171/1‑151 (MQ=255)
                                                                                      gCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAAc                                                           <  1:1187291/149‑1 (MQ=255)
                                                                                      gCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAAc                                                           <  1:389488/149‑1 (MQ=255)
                                                                                        ggaaTCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTg                                                        <  2:480311/148‑1 (MQ=255)
                                                                                         tAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTgg                                                       >  2:1273213/1‑150 (MQ=255)
                                                                                          aaTCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATGACCGGCGGTAACGCGTTATGCAGAACGTGGc                                                      <  1:108000/150‑1 (MQ=255)
                                                                                          aaTCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGc                                                      <  1:202683/150‑1 (MQ=255)
                                                                                           aTCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGa                                                    <  2:1323997/151‑1 (MQ=255)
                                                                                             ccTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGAc                                                   >  2:9689/1‑150 (MQ=255)
                                                                                              ctgctgGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGc                                                 >  1:476565/1‑151 (MQ=255)
                                                                                                gctgGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCa                                                <  1:357521/150‑1 (MQ=255)
                                                                                                   ggCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCaaaa                                             <  1:1040478/150‑1 (MQ=255)
                                                                                                   ggCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAAt                                            <  1:235819/151‑1 (MQ=255)
                                                                                                    gCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATa                                           <  2:459354/151‑1 (MQ=255)
                                                                                                     cGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAAt                                            >  2:363160/1‑149 (MQ=255)
                                                                                                       gATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTa                                        <  1:1147881/151‑1 (MQ=255)
                                                                                                          ggCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTaaa                                      <  1:1049409/150‑1 (MQ=255)
                                                                                                             gTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGCGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAgcgc                                  >  1:879648/1‑151 (MQ=255)
                                                                                                                        aGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACcgcg                       >  1:1257552/1‑151 (MQ=255)
                                                                                                                        aGCGCCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACcgcg                       >  1:78596/1‑151 (MQ=255)
                                                                                                                         gcgcCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCgg                      <  1:957799/151‑1 (MQ=255)
                                                                                                                           gcCCTAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTGAAGCGCGACAAACCGCGGGt                    >  1:758952/1‑151 (MQ=255)
                                                                                                                               tAAACCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTcgcc                >  2:425426/1‑151 (MQ=255)
                                                                                                                                aaaCCCGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTcgcc                <  2:988320/150‑1 (MQ=255)
                                                                                                                                   cccGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTcgccgccg            <  2:41552/151‑1 (MQ=255)
                                                                                                                                   accGTCCCGCTAAGAAACATTTCGGGGATCATCGCCAGCGTCAGCATGGCGTAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGGAACTCGTGATGCATAACGTGGCTACGCAAAATAGGAAGGCGCGACAAGCcgcg                       >  2:681372/2‑140 (MQ=255)
                                                                                                                                    ccGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTcgccgccgc           <  1:590516/151‑1 (MQ=255)
                                                                                                                                     cGGCCAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCt          >  1:483663/1‑151 (MQ=255)
                                                                                                                                        ccAGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCtttt       <  2:1071799/151‑1 (MQ=255)
                                                                                                                                          aGCTAAAAAACATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTa      >  2:294279/1‑150 (MQ=255)
                                                                                                                                              aaaaaaCATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACAt   <  2:270054/149‑1 (MQ=255)
                                                                                                                                              aaaaaaCATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACAt   >  1:178201/1‑149 (MQ=255)
                                                                                                                                              aaaaaaCATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATa  <  1:1003535/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GGTTAGCCATGGCACCCAAAATATCAGAAATCACGTTAACAATAATGACCAGTGAGCCACACACCATCACTCCGGCGGAAATGGCTGCGTAATCCTGCTGGCGGATGGCGTTAATTAACCAGCGCCCTAAACCCGGCCAGCTAAAAACCATTTCGGTGATCATCGCCAGCGTCAGCATGGTGGAAAACTGTAAGCCCAGGCGAGGAATAACCGGCGGTAACGCGTTATGCAGAACGTGGCGACGCAAAATAGTAAAGCGCGACAAACCGCGGGTCGCCGCCGCTTTTACATA  >  NC_000913/1354540‑1354831

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