Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 520,993 T→G 100% V193G (GTC→GGC)  ybbP → putative ABC transporter permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913520,9930TG98.2% 572.6 / ‑1.4 174V193G (GTC→GGC) ybbPputative ABC transporter permease
Reads supporting (aligned to +/- strand):  ref base T (1/4);  new base G (95/73);  total (96/78)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 177 total aligned reads displayed.

GGATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGTCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACGG  >  NC_000913/520898‑521089
                                                                                               |                                                                                                
ggATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCAcc                                                                                              >  1:2773451/1‑100 (MQ=255)
   ttCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTgg                                                                                           >  2:2314570/1‑100 (MQ=255)
     ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGTGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                         >  2:155636/1‑100 (MQ=255)
     ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                         >  2:1082377/1‑100 (MQ=255)
     ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                         >  1:1630109/1‑100 (MQ=255)
     ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                         >  1:333696/1‑100 (MQ=255)
     ccGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCg                                                                                         >  2:2872453/1‑100 (MQ=255)
      cGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGt                                                                                        >  2:2219421/1‑100 (MQ=255)
       ggTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGtt                                                                                       >  1:789264/1‑100 (MQ=255)
        gTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTta                                                                                      <  2:1089749/100‑1 (MQ=255)
          tttAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCCCCTGGCGTtata                                                                                    >  1:3945550/1‑100 (MQ=255)
          tttAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTtata                                                                                    <  1:4052774/100‑1 (MQ=255)
          tttAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTtata                                                                                    <  2:3621485/100‑1 (MQ=255)
           ttAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATaa                                                                                   >  1:325062/1‑100 (MQ=255)
           ttAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATaa                                                                                   >  2:1634329/1‑100 (MQ=255)
            tAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATaaa                                                                                  >  2:1522460/1‑100 (MQ=255)
             aaCCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAAt                                                                                 <  1:874026/100‑1 (MQ=255)
                cccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcg                                                                              >  1:2711355/1‑100 (MQ=255)
                cccTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcg                                                                              <  1:2872469/100‑1 (MQ=255)
                  cTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggc                                                                            >  1:3275711/1‑100 (MQ=255)
                   ttCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGTCACCTGGCGTTATAAATTcggcg                                                                           <  1:4047146/100‑1 (MQ=255)
                   ttCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggcg                                                                           <  1:1706113/100‑1 (MQ=255)
                    tCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggcgg                                                                          <  2:149711/100‑1 (MQ=255)
                    tCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTcggcgg                                                                          >  1:2811731/1‑100 (MQ=255)
                     ccAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCggcggc                                                                         <  1:3837487/100‑1 (MQ=255)
                         aTGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACg                                                                     >  1:3590352/1‑100 (MQ=255)
                          tGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGCCGGGGCACCTGGCGTTATAAATTCGGCGGCAACga                                                                    >  2:3986769/1‑100 (MQ=255)
                            gCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGACGGGGCACATGGCGTTATAAATTCGGCGGCAACgaga                                                                  >  2:2765436/1‑100 (MQ=255)
                             ctgcgcgTCTGATGATGAAGCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGaa                                                                 <  2:3297791/97‑1 (MQ=255)
                             cTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGaa                                                                 >  1:2761878/1‑100 (MQ=255)
                              tCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAAc                                                                <  2:3733191/100‑1 (MQ=255)
                                cgcgTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCa                                                              >  1:2706139/1‑100 (MQ=255)
                                cgcgTCTGATGATGAAGCTGGCGGATGTCGATAAACTCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACgag                                                                   <  2:1696523/95‑1 (MQ=255)
                                  cgTCTGATGATGAATCTGGCTGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGc                                                            <  2:3873506/100‑1 (MQ=255)
                                  cgTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGc                                                            <  1:518195/100‑1 (MQ=255)
                                  cgTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGc                                                            <  2:2738190/100‑1 (MQ=255)
                                     cTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCGAGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCg                                                         >  2:2288718/1‑100 (MQ=255)
                                      tgatgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGa                                                        <  1:3786809/100‑1 (MQ=255)
                                       gatgatgaCTCTGTCGTAGGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCg                                                         <  2:2353509/98‑1 (MQ=255)
                                        atgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACg                                                      >  1:169950/1‑100 (MQ=255)
                                        atgatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACg                                                      >  1:1484870/1‑100 (MQ=255)
                                          gatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGTCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGc                                                    >  1:458177/1‑100 (MQ=255)
                                          gatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACg                                                      >  2:2626868/1‑98 (MQ=255)
                                          gatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGc                                                    >  1:2627614/1‑100 (MQ=255)
                                          gatgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGc                                                    <  1:3407931/100‑1 (MQ=255)
                                           atgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCt                                                   >  2:1836000/1‑100 (MQ=255)
                                           atgaATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGAGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGc                                                    >  1:899156/1‑99 (MQ=255)
                                             gaagatGGCGGATGTCGATAGACCTGGAGCGGTGCAGCCCGGGAGTCGGGGCACCAGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTAt                                                 <  2:3682643/95‑1 (MQ=255)
                                              aaTCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCCGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATg                                                <  2:2241927/100‑1 (MQ=255)
                                               aTCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATga                                               >  2:45620/1‑100 (MQ=255)
                                                  tGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGTTCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATgag                                              <  2:1011844/98‑1 (MQ=255)
                                                  tGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGaa                                            >  1:629794/1‑100 (MQ=255)
                                                   ggCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGaaa                                           >  1:1060640/1‑100 (MQ=255)
                                                     cGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATg                                         <  1:329617/100‑1 (MQ=255)
                                                     cGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATg                                         >  1:3204433/1‑100 (MQ=255)
                                                     cGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATg                                         <  1:2498268/100‑1 (MQ=255)
                                                          gTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGAGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttgt                                    >  2:56396/1‑100 (MQ=255)
                                                           tCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGgttgtt                                   <  1:1780567/100‑1 (MQ=255)
                                                            cGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTa                                  >  2:533091/1‑100 (MQ=255)
                                                            cGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTa                                  >  1:2609004/1‑100 (MQ=255)
                                                            cGATAAAACCGGAGCCGTGCAGCCGGGGAGGCGGGGCACCTGGCCTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTa                                  >  1:1087519/1‑100 (MQ=255)
                                                             gATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTAc                                 <  1:2650114/100‑1 (MQ=255)
                                                              aTAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTAcc                                <  1:1169274/100‑1 (MQ=255)
                                                              aTAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGGCAACTGGCGATATAAATTCGGCGGCGACGAGAACCAGCGCGACGGCTAAGAGAAATgggggctgc                                 >  2:4115220/1‑92 (MQ=255)
                                                              aTAAAACCGGAGACGGGCAGCCGGGGAGTCGGGCCACCTGGCGTCATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTAcc                                <  2:2447492/100‑1 (MQ=255)
                                                                aaaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCt                               >  2:1335002/1‑99 (MQ=255)
                                                                 aaaCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCa                             >  1:2432180/1‑100 (MQ=255)
                                                                  naCCGGAGCCGTGCAGCCGGGGAGTCGGGGCACCCGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTAcccaa                             >  2:2183167/2‑96 (MQ=255)
                                                                    ccGGAGCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCt                          <  1:240151/100‑1 (MQ=255)
                                                                      ggAGCCGTGCAGCCGGGGAGTCGGGGTACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGc                           <  2:585740/97‑1 (MQ=255)
                                                                         gCCGTGCAGCCGGGGAGTCGGGGCACCTTGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAc                     <  2:3484115/100‑1 (MQ=255)
                                                                         gCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAc                     <  2:303570/100‑1 (MQ=255)
                                                                         gCCGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAc                     <  1:3125821/100‑1 (MQ=255)
                                                                         gCCGTGCAGCCGGGGAGGCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCATAAAc                     >  2:3020986/1‑100 (MQ=255)
                                                                          ccGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAcc                    <  2:928324/100‑1 (MQ=255)
                                                                          ccGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAcc                    <  1:45620/100‑1 (MQ=255)
                                                                           cGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAccc                   <  2:882943/100‑1 (MQ=255)
                                                                           cGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAccc                   >  1:2086429/1‑100 (MQ=255)
                                                                           cGTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAccc                   >  1:1960571/1‑100 (MQ=255)
                                                                            gTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCg                  >  1:1578656/1‑100 (MQ=255)
                                                                            gTGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAAGGGTTGTTGCCTCAGCTTAAACCCg                  >  2:610230/1‑100 (MQ=255)
                                                                             tGCAGCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGGACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAAACCGa                 >  2:419816/1‑100 (MQ=255)
                                                                              gCAGCCGGGGAGTCGGGGCACCTGGCGTAATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaa                >  1:705570/1‑100 (MQ=255)
                                                                               cggCCGGGGGGTCGGGTCACCTGGCGTTATAAATGCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaac               <  2:757652/98‑1 (MQ=255)
                                                                                 gCCGGGGAGTCGGGTCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaac               <  2:2418300/98‑1 (MQ=255)
                                                                                 gCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaacaa             >  1:2421353/1‑100 (MQ=255)
                                                                                 gCCGGGGAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGaacaa             >  2:3956576/1‑100 (MQ=255)
                                                                                    ggggAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGc          >  2:3098662/1‑100 (MQ=255)
                                                                                    ggggAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGc          >  1:3078491/1‑100 (MQ=255)
                                                                                     gggAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCt         <  1:1665923/100‑1 (MQ=255)
                                                                                     gggAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCt         <  1:2295320/100‑1 (MQ=255)
                                                                                       gAGTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTgg       <  2:2617202/100‑1 (MQ=255)
                                                                                        aGTCGGGGCACCTGGCGGTATAAATTCGGCGGCAACGAGCACCCGCTCGACGGCCATGAGAACGGACTGTTACCTCTGCTTGAACCCGACCAACg           >  2:2179880/1‑95 (MQ=255)
                                                                                         gTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCCGCCCGACGGCTATGAGAGATGGTCATTACCTCAGCCTAAACCCGTCCAACGc          >  2:3909168/1‑95 (MQ=255)
                                                                                         gTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTa     <  1:468581/100‑1 (MQ=255)
                                                                                         gTCGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTa     <  2:2086416/100‑1 (MQ=255)
                                                                                           cGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACTGCTATGATAAATGGCAGTTACCACAGCTTCAACCCGAACAACGCTGGCAc    >  2:1628920/1‑99 (MQ=255)
                                                                                           cGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACg   >  1:1400089/1‑100 (MQ=255)
                                                                                           cGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACg   >  1:4094833/1‑100 (MQ=255)
                                                                                           cGGGGCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGCTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACg   >  2:1383327/1‑100 (MQ=255)
                                                                                               |                                                                                                
GGATTCCGGTTTTAACCCCTTCCAGATGGCTCCGCGTCTGATGATGAATCTGGCGGATGTCGATAAAACCGGAGCCGTGCAGCCGGGGAGTCGGGTCACCTGGCGTTATAAATTCGGCGGCAACGAGAACCAGCTCGACGGCTATGAGAAATGGTTGTTACCTCAGCTTAAACCCGAACAACGCTGGTACGG  >  NC_000913/520898‑521089

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