Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,904,595 T→G 17.1% I234S (ATT→AGT)  manZ → mannose‑specific enzyme IID component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,904,5950TG17.1% 253.7 / 45.0 104I234S (ATT→AGT) manZmannose‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base T (44/42);  new base G (7/11);  total (51/53)
Fisher's exact test for biased strand distribution p-value = 4.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 106 total aligned reads displayed.

GGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTG  >  NC_000913/1904496‑1904693
                                                                                                    |                                                                                                  
gggCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTAt                                                                                                    >  1:843108/1‑100 (MQ=255)
 ggCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTAtt                                                                                                   >  1:2047392/1‑100 (MQ=255)
   cATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCt                                                                                                 >  1:579786/1‑100 (MQ=255)
       ggTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGAc                                                                                             <  2:860777/100‑1 (MQ=255)
         ttAACGAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcag                                                                                           <  1:448882/100‑3 (MQ=255)
         ttAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCa                                                                                           >  2:2196377/1‑100 (MQ=255)
         ttAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcag                                                                                           <  2:362724/100‑3 (MQ=255)
         ttAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcag                                                                                           <  1:1840384/100‑3 (MQ=255)
         ttAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcag                                                                                           <  2:33355/100‑3 (MQ=255)
         ttAACAAGTGGAAACAGGTCAACATCCCGCTGGTAGTCTCTCGCATTACTGACCAGACCGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcag                                                                                           <  1:2284732/100‑3 (MQ=255)
         tgaaCAAGTGGACACATGTCAACAGCCCGCTAGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAAGACGTTACTACTTGCCAGAATAGTCTGGAcag                                                                                           <  1:233356/98‑3 (MQ=255)
          tAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAg                                                                                          <  2:1975267/100‑1 (MQ=255)
            aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGtt                                                                                        <  1:153075/100‑1 (MQ=255)
            aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGtt                                                                                        >  1:735986/1‑100 (MQ=255)
            aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcagtag                                                                                        <  2:614586/100‑6 (MQ=255)
            aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcagtag                                                                                        <  1:1423377/100‑6 (MQ=255)
             cAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTa                                                                                       <  1:2333986/100‑1 (MQ=255)
             cAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTa                                                                                       >  1:1443442/1‑100 (MQ=255)
             cAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATGCTGGACCAGTTa                                                                                       >  1:381725/1‑100 (MQ=255)
               aGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAAt                                                                                     <  2:371353/100‑1 (MQ=255)
                gTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATg                                                                                    <  2:2385329/100‑1 (MQ=255)
                gTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATg                                                                                    <  1:731732/100‑1 (MQ=255)
                atGGACATATGGCAGCATCCCTCTGTTTCTCTCGCGCATTACTGACCATACGGTCAAAGAACACGTTACTACTGTCCAGAATAGTCTGGAcagtagta                                                                                      <  1:2093492/97‑8 (MQ=19)
                 tGGACACCTGTCACCATCCCGCTGGTTGTCTCTCGCATTACTGACCATACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGc                                                                                   <  1:51310/100‑1 (MQ=255)
                  ggACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGGTCCAGAATAGTCTGGAcagtagtaacg                                                                                   <  1:2285878/100‑11 (MQ=25)
                  ggACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGcc                                                                                  >  1:444398/1‑100 (MQ=255)
                  ggACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGcc                                                                                  >  2:178347/1‑100 (MQ=255)
                   gACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCa                                                                                 >  1:606477/1‑100 (MQ=255)
                    acacATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAg                                                                                >  2:1914731/1‑100 (MQ=255)
                    acacATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAg                                                                                <  1:178346/100‑1 (MQ=255)
                            cAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTa                                                                        <  2:2002576/100‑1 (MQ=255)
                            cAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTa                                                                        >  1:2141230/1‑100 (MQ=255)
                            cAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTa                                                                        >  2:2589931/1‑100 (MQ=255)
                                   ccGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgc                                                                 >  2:732505/1‑100 (MQ=255)
                                    cGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgc                                                                 <  1:1611199/99‑1 (MQ=255)
                                      cTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTGTTCTGGGCCAGTTAATGCCAGGCCTGGTACCActgctgc                                                              <  1:1464716/100‑1 (MQ=255)
                                      cTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctgc                                                              <  1:1777731/100‑1 (MQ=255)
                                       tGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctgct                                                             >  1:2212839/1‑100 (MQ=255)
                                        ggTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctgctg                                                            <  2:1443450/100‑1 (MQ=255)
                                         gTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGATCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGa                                                           <  2:896013/100‑1 (MQ=255)
                                         gTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGa                                                           >  2:1588320/1‑100 (MQ=255)
                                           tGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAcc                                                         <  2:1023423/100‑1 (MQ=255)
                                           tGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAcc                                                         <  1:95100/100‑1 (MQ=255)
                                            gTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCt                                                        <  2:1483373/100‑1 (MQ=255)
                                              ctctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCttt                                                      >  2:2494961/1‑100 (MQ=255)
                                               tctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTg                                                     <  2:681288/100‑1 (MQ=255)
                                                ctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGc                                                    >  2:56444/1‑100 (MQ=255)
                                                cacGCATTACTTACCCGACGTGCAAAGAACTCCTTACTACTGTCTAGTCAATTCTGGACCAGGTAATGTTAGGCCTGGTACCACTGCTGCTGACCTTTGc                                                    >  1:214927/3‑100 (MQ=255)
                                                 tcGCATTACTGACCAGACGGGCAAAGAACACTTTACTACTGACCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCt                                                   <  2:511912/100‑1 (MQ=255)
                                                  cGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTTGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCtt                                                  <  1:1439263/100‑1 (MQ=255)
                                                  cGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCtt                                                  >  1:1142038/1‑100 (MQ=255)
                                                   gCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTg                                                 >  1:1374872/1‑100 (MQ=255)
                                                     aTTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTa                                               <  2:595883/100‑1 (MQ=255)
                                                       tACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTAtg                                             <  1:837474/100‑1 (MQ=255)
                                                       tACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTAtg                                             >  2:1443985/1‑100 (MQ=255)
                                                       gatctaCCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTAtg                                             <  1:349716/95‑1 (MQ=255)
                                                         cTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTAtgtg                                           <  1:209310/100‑1 (MQ=255)
                                                           gACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGc                                         >  2:45858/1‑100 (MQ=255)
                                                            aCCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTAtgagatc                                        >  2:349718/1‑95 (MQ=255)
                                                            aCCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGc                                         <  1:2196248/99‑1 (MQ=255)
                                                              cAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTAc                                      <  2:766376/100‑1 (MQ=255)
                                                                gACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTg                                    >  1:1019301/1‑100 (MQ=255)
                                                                  cGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTgcg                                  >  1:1054673/1‑100 (MQ=255)
                                                                   gggCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTgcgc                                 >  2:60937/1‑100 (MQ=255)
                                                                     gCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCaa                               >  2:1398453/1‑100 (MQ=255)
                                                                      cAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCaaa                              >  2:142297/1‑100 (MQ=255)
                                                                      cAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCaaa                              >  1:1193731/1‑100 (MQ=255)
                                                                      cAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCGTTGCTTGTATGTGGCTACTGCGCaa                               >  2:1288169/1‑99 (MQ=255)
                                                                        aaGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCaaaaa                            >  2:2023481/1‑100 (MQ=255)
                                                                          gAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAg                          <  2:1675553/100‑1 (MQ=255)
                                                                            acacGTTACTAGTGTCCAGACTATGCTGGACCAGTTAATGCTAGGCCTGGTACCACTGCTGCTTACCTTTTCTTGTATGTGGCTACTGCGCAAAAAAg                          <  1:1932692/98‑1 (MQ=255)
                                                                             cacGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTa                       >  2:1340134/1‑100 (MQ=255)
                                                                                 ttACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAccc                   <  2:1243699/100‑1 (MQ=255)
                                                                                 ttACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAccc                   >  2:86831/1‑100 (MQ=255)
                                                                                  attaactggTCCAGAATAGTCTGGACAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAccc                   >  1:243767/9‑99 (MQ=37)
                                                                                   actactGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGc                 <  1:1463520/100‑1 (MQ=255)
                                                                                    ctactGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCt                >  1:2188539/1‑100 (MQ=255)
                                                                                    ctactGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCt                >  1:899629/1‑100 (MQ=255)
                                                                                     aactggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCt                >  1:2102369/6‑100 (MQ=255)
                                                                                      actggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtg               >  2:2204356/5‑100 (MQ=255)
                                                                                      actGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtgt              <  1:56443/100‑1 (MQ=255)
                                                                                       ctggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtgt              >  2:2003800/4‑100 (MQ=255)
                                                                                       ctGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtgtg             >  2:236552/1‑100 (MQ=255)
                                                                                       ctGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtgtg             >  2:2193857/1‑100 (MQ=255)
                                                                                       ctGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCGTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtgtg             <  1:529875/100‑1 (MQ=255)
                                                                                        tggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtgtg             >  2:225892/3‑100 (MQ=255)
                                                                                        tggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCtgtg             >  1:2231371/3‑100 (MQ=255)
                                                                                         ggTCCAGAATAGTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTgg            >  2:1166866/2‑100 (MQ=255)
                                                                                          gTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGa           >  1:648357/1‑100 (MQ=255)
                                                                                          gTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGa           >  2:1061911/1‑100 (MQ=255)
                                                                                          gTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGCCTACTGCGCAAAAAAGTTAACCCGCTGTGGa           <  2:1740599/100‑1 (MQ=255)
                                                                                             cAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGatca        <  2:1461568/100‑1 (MQ=255)
                                                                                              aGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGatcat       <  1:1525384/100‑1 (MQ=255)
                                                                                              aGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGatcat       >  2:1310949/1‑100 (MQ=255)
                                                                                              aGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGatcat       <  2:1067295/100‑1 (MQ=255)
                                                                                              aGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGatcat       >  2:2242776/1‑100 (MQ=255)
                                                                                                 cTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGt    <  2:2400654/100‑1 (MQ=255)
                                                                                                   aTTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTg  <  1:1251782/100‑1 (MQ=255)
                                                                                                   aTTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTg  <  1:2483155/100‑1 (MQ=255)
                                                                                                   aTTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATCTGGCTGCTGCGCAAAAAAGTTAACCCGCTGTGGatcatctt    >  1:719151/1‑96 (MQ=255)
                                                                                                    |                                                                                                  
GGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTG  >  NC_000913/1904496‑1904693

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