Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,681,932 T→C 100% P33P (CCA→CCG ydgC ← GlpM family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,681,9320TC97.5% 389.2 / ‑1.1 117P33P (CCA→CCGydgCGlpM family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (58/55);  total (60/57)
Fisher's exact test for biased strand distribution p-value = 6.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.36e-01

Only 100 of 117 total aligned reads displayed.

TCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGCCCCA  >  NC_000913/1681835‑1682029
                                                                                                 |                                                                                                 
tCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCgg                                                                                                 <  1:2514311/100‑1 (MQ=255)
  gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCggta                                                                                               <  1:2829417/100‑3 (MQ=255)
  gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                               <  1:3112675/100‑1 (MQ=255)
  gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                               >  2:866443/1‑100 (MQ=255)
  gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                               <  1:1639743/100‑1 (MQ=255)
       caTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGc                                                                                          <  1:1476409/100‑1 (MQ=255)
       caTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGc                                                                                          >  1:2909884/1‑100 (MQ=255)
        aTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGcc                                                                                         <  1:1776320/100‑1 (MQ=255)
         tACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCACAGGTCGGGAAAAGCGGAATCAGcc                                                                                         >  2:4157308/1‑99 (MQ=255)
          aCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCg                                                                                       >  1:494671/1‑100 (MQ=255)
           cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCgg                                                                                      >  1:610449/1‑100 (MQ=255)
           cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCgg                                                                                      >  1:2220824/1‑100 (MQ=255)
           cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCgg                                                                                      >  2:3176505/1‑100 (MQ=255)
           cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCgg                                                                                      >  2:1192268/1‑100 (MQ=255)
           cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGAGGAATCAGCCCgg                                                                                      >  1:2422602/1‑100 (MQ=255)
            tAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGc                                                                                     >  1:229965/1‑100 (MQ=255)
            tAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGc                                                                                     >  1:230589/1‑100 (MQ=255)
            tAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGc                                                                                     >  2:3955906/1‑100 (MQ=255)
             aaaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGc                                                                                     <  1:527145/99‑1 (MQ=255)
                aaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGata                                                                                 <  2:3755476/100‑1 (MQ=255)
                  tgatgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGAtata                                                                               <  1:3243112/100‑1 (MQ=255)
                  tgatgGTAGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCACAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGAtata                                                                               <  2:2231662/100‑1 (MQ=255)
                      ggTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAAGCAGCCCGGCGAGAtaataa                                                                           >  1:3321448/1‑100 (MQ=255)
                         tGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttt                                                                        <  1:3556505/100‑1 (MQ=255)
                         tGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCACGGCGATATAATAAttt                                                                        >  2:328689/1‑100 (MQ=255)
                          gcgcGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAtttt                                                                       >  1:1435339/1‑100 (MQ=255)
                          gcgcGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAACCGGAATCAGCCCGGCGATATAATAAtttt                                                                       >  2:1071235/1‑100 (MQ=255)
                           cgcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttttt                                                                      >  1:2117600/1‑100 (MQ=255)
                           cgcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttttt                                                                      >  2:122693/1‑100 (MQ=255)
                            gcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTc                                                                     <  2:4003941/100‑1 (MQ=255)
                            gcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCAGAATCAGCCCGGCGATATAATAATTTTTc                                                                     >  2:1611935/1‑100 (MQ=255)
                             cgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCg                                                                    >  2:868259/1‑100 (MQ=255)
                               tAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGtt                                                                  >  2:3443345/1‑100 (MQ=255)
                                 aGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGtttt                                                                <  2:3896095/100‑1 (MQ=255)
                                  ggatccAATGCCCCATTCGCCGGAAACTATATAATGCACCATGCGCGCAAAAGTCCGGAAACGCCGAATCAGTCCGGCGATATAATAATTTTTCg                                                                    <  2:3236238/90‑1 (MQ=255)
                                  ggCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGttttt                                                               >  2:905304/1‑100 (MQ=255)
                                  ggCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGttttt                                                               >  1:4140411/1‑100 (MQ=255)
                                   gCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg                                                              <  1:2134487/100‑1 (MQ=255)
                                   gCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg                                                              <  2:683286/100‑1 (MQ=255)
                                    cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                             <  1:4029921/100‑1 (MQ=255)
                                    cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                             <  2:4035649/100‑1 (MQ=255)
                                     ttCAATGCCGCGTTCGCTGGCAACAATATAATGCTCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCt                                                            <  1:2493386/100‑1 (MQ=255)
                                     ttCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCt                                                            <  2:1539113/100‑1 (MQ=255)
                                      tCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTa                                                           >  1:3488155/1‑100 (MQ=255)
                                         aTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaaa                                                        >  1:1581685/1‑100 (MQ=255)
                                          tGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAac                                                       <  2:3376850/100‑1 (MQ=255)
                                          tGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAac                                                       <  1:25167/100‑1 (MQ=255)
                                          tGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAac                                                       <  1:1512325/100‑1 (MQ=255)
                                           gCCGCGTGCGCAGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCGGGCGATATAATAATTTTTCGTTTTTGCTAAAaca                                                      >  2:2189624/1‑100 (MQ=255)
                                             cgcgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACAcc                                                    <  1:3063264/100‑1 (MQ=255)
                                               cgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCaa                                                  >  1:2528388/1‑100 (MQ=255)
                                                gTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAAt                                                 <  2:4140377/100‑1 (MQ=255)
                                                 ttCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATc                                                <  2:2066155/100‑1 (MQ=255)
                                                 ttCGCTGGCAACAATATAATGCGCGATAAGCACAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATc                                                <  1:40000/100‑1 (MQ=255)
                                                  tCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCa                                               <  2:580238/100‑1 (MQ=255)
                                                  tCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCa                                               <  2:580601/100‑1 (MQ=255)
                                                  tCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCa                                               <  2:2704874/100‑1 (MQ=255)
                                                  tCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCa                                               <  1:1884669/100‑1 (MQ=255)
                                                   cGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCaa                                              >  2:2710048/1‑100 (MQ=255)
                                                   cGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCaa                                              >  1:3774492/1‑100 (MQ=255)
                                                    gCTGGCAACAATATAATGCGCGATACGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAAc                                             <  1:4173095/100‑1 (MQ=255)
                                                    gCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAAc                                             >  1:1718381/1‑100 (MQ=255)
                                                    gCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAAc                                             >  2:2165974/1‑100 (MQ=255)
                                                     ccggCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAACAGAATCAGCCCGGCGCTAAAATAACTTTTCGTTTTTGCTAAAACACCAATCa                                               >  2:1850104/3‑97 (MQ=255)
                                                      tGGCAACAATATAATGCGCGATAAGCGCAAAGGTCTGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCaa                                              <  2:850247/97‑1 (MQ=255)
                                                      tGGAAACAATATAGTGCGCGATAAGGGCAAAGGTCGGGAAAAGCGGAAGCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACa                                            <  2:1740625/99‑1 (MQ=255)
                                                       ggCAACAAGATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACa                                            <  2:2829262/98‑1 (MQ=255)
                                                         caacaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCAc                                        >  1:3293968/1‑100 (MQ=255)
                                                          aaaaaTATAATGGGCGATAAGCGCAAAGGGCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCAc                                        <  2:542435/96‑1 (MQ=255)
                                                           acaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTTTTTTTTTTCCAATCAACAGCACTa                                      >  2:3715003/1‑100 (MQ=255)
                                                              atatAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTGTGCTAAAACACCAATCAACAGCACTAcaa                                   >  1:3962520/1‑98 (MQ=255)
                                                                ataATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAgc                                 <  2:3785680/100‑1 (MQ=255)
                                                                 taATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAgcg                                >  1:1011538/1‑100 (MQ=255)
                                                                   aTGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCa                              >  1:3105995/1‑100 (MQ=255)
                                                                   aTGCGCGAGAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGGTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCa                              <  2:1477794/100‑1 (MQ=255)
                                                                    tGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGTGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCAc                             <  1:624509/100‑1 (MQ=255)
                                                                    tGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCAc                             <  2:1049230/100‑1 (MQ=255)
                                                                     gcgcGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCAcc                            >  2:2055436/1‑100 (MQ=255)
                                                                      cgcgATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCg                           <  2:532768/100‑1 (MQ=255)
                                                                      cgcgATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCg                           >  1:2090280/1‑100 (MQ=255)
                                                                      cgcgATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCg                           >  1:3850375/1‑100 (MQ=255)
                                                                          aTAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAgcg                       <  2:3056125/100‑1 (MQ=255)
                                                                             aGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCgg                    >  2:3034780/1‑100 (MQ=255)
                                                                              gcgcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTTCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCgg                    >  2:3654621/1‑99 (MQ=255)
                                                                                gcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCtt                 >  2:1759676/1‑100 (MQ=255)
                                                                                 cAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCttt                >  1:1256828/1‑100 (MQ=255)
                                                                                  aaaGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTg               <  1:328692/100‑1 (MQ=255)
                                                                                   aaGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGa              <  1:1849004/100‑1 (MQ=255)
                                                                                   aaGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGa              <  2:3002250/100‑1 (MQ=255)
                                                                                   aaGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGa              <  2:1991321/100‑1 (MQ=255)
                                                                                      gTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACAACAATCAACACAACAACCAGAGACCCGAGCGCGGcaatga              >  1:3901527/1‑92 (MQ=255)
                                                                                      gTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTa           >  1:810485/1‑100 (MQ=255)
                                                                                        cGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTAcc         <  2:3774463/100‑1 (MQ=255)
                                                                                          ggAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAg       >  2:698845/1‑100 (MQ=255)
                                                                                          ggAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTGTTGCTAAAACACCAATCAACAACAATACCAGCGCTCCGAGCGCGGCTTTGAttggtgt       >  2:2982642/1‑95 (MQ=255)
                                                                                           gAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGc      >  1:2608882/1‑100 (MQ=255)
                                                                                           gAAAAGCGGAATCAGCACGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCCCCGAGCGCGGCTTTGATTACCAg       >  2:3651408/1‑99 (MQ=255)
                                                                                            aaaaGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGcc     <  2:822629/100‑1 (MQ=255)
                                                                                             aaaGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGccc    >  1:1898634/1‑100 (MQ=255)
                                                                                               aGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGCCCCa  <  1:3361353/100‑1 (MQ=255)
                                                                                                 |                                                                                                 
TCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGCCCCA  >  NC_000913/1681835‑1682029

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