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,9320TC99.3% 501.8 / ‑5.1 143P33P (CCA→CCGydgCGlpM family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (63/79);  total (63/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 143 total aligned reads displayed.

ATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGCCCC  >  NC_000913/1681834‑1682028
                                                                                                  |                                                                                                
 tCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCgg                                                                                                <  2:1660589/100‑1 (MQ=255)
   gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                              <  2:1900821/100‑1 (MQ=255)
   gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                              >  1:2038628/1‑100 (MQ=255)
     ccACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATc                                                                                            <  2:2655319/100‑1 (MQ=255)
       acaTACTAAAACTGATGGTTGCGCGTCAGGCCTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATACGCGCAAAGGTCGGGAAAAGCGGAACCa                                                                                           >  1:2064767/1‑99 (MQ=255)
       acaTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAg                                                                                          >  1:4165624/1‑100 (MQ=255)
        caTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAAGGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGc                                                                                         <  1:86673/100‑1 (MQ=255)
          tACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGccc                                                                                       >  2:4174629/1‑100 (MQ=255)
          tACTAAAAATGATGGTTGCGCGTAAGGCGTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGccc                                                                                       <  1:2916677/100‑1 (MQ=255)
           aCTAAAAATGATGTTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGc                                                                                         <  1:1407218/97‑1 (MQ=255)
            cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCgg                                                                                     >  2:513364/1‑100 (MQ=255)
            cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCgg                                                                                     <  1:575267/100‑1 (MQ=255)
            cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCgg                                                                                     >  2:425389/1‑100 (MQ=255)
            cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGCTCGCTGGCCACAATATAACCCGCGATAAGCGCAAAGGTCCGCACAAGCGGACTctcccc                                                                                       >  1:3459376/1‑93 (MQ=255)
             tAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGc                                                                                    >  1:2099055/1‑100 (MQ=255)
             tAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGc                                                                                    >  2:3925620/1‑100 (MQ=255)
             tAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGc                                                                                    <  2:2265401/100‑1 (MQ=255)
             tAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGc                                                                                    <  1:1712430/100‑1 (MQ=255)
             tAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGc                                                                                    >  2:787253/1‑100 (MQ=255)
              aaaatTGACGGTTGAGCGAAAGGCTCCAATGCCGCGTTCGCAGGCAACAATATAATGCGCGATAAGCGAAAAGGTGGTGAACAGCGTAATCAGCCCGGc                                                                                    <  1:1419463/99‑1 (MQ=255)
                aaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTGCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGat                                                                                 <  2:2204965/100‑1 (MQ=255)
                 aaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGata                                                                                <  2:4253622/100‑1 (MQ=255)
                 aaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGata                                                                                <  2:491802/100‑1 (MQ=255)
                  atgatgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGatat                                                                               <  2:1026876/100‑1 (MQ=255)
                  atgatgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGatat                                                                               <  1:1413507/100‑1 (MQ=255)
                     atgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATataat                                                                            <  2:131959/100‑1 (MQ=255)
                       ggTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATAtaataa                                                                          <  2:1050812/100‑1 (MQ=255)
                        gTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATaataat                                                                         >  2:392791/1‑100 (MQ=255)
                         tggcgcGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATaataat                                                                         <  1:2010696/97‑1 (MQ=255)
                          tGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttt                                                                       <  1:3168652/100‑1 (MQ=255)
                          tGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttt                                                                       >  2:1714085/1‑100 (MQ=255)
                              cgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCg                                                                   <  2:151079/100‑1 (MQ=255)
                              cgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCg                                                                   >  2:3410179/1‑100 (MQ=255)
                              cgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCg                                                                   >  1:599016/1‑100 (MQ=255)
                              cgTAAGGCTTCAAGGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCg                                                                   <  1:2790179/100‑1 (MQ=255)
                                tAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGtt                                                                 >  1:1986767/1‑100 (MQ=255)
                                 aaGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTTTTTAAAAGCGTTATCAGCCCGGCGATATAATAATTTTTCGttt                                                                <  1:2707076/100‑1 (MQ=255)
                                 aaGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGAGATATAATAATTTTTCGttt                                                                >  1:3420312/1‑100 (MQ=255)
                                    gCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg                                                             >  2:530521/1‑100 (MQ=255)
                                    gCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg                                                             <  2:1166069/100‑1 (MQ=255)
                                    gCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                            >  2:4030364/1‑100 (MQ=255)
                                     cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg                                                             >  2:24187/1‑99 (MQ=255)
                                     cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                            <  1:673448/100‑1 (MQ=255)
                                     cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                            <  2:334922/100‑1 (MQ=255)
                                     cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                            <  2:653020/100‑1 (MQ=255)
                                      ttCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCt                                                           <  1:186561/100‑1 (MQ=255)
                                       tCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTa                                                          >  2:3634184/1‑100 (MQ=255)
                                       tCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTa                                                          <  1:1672462/100‑1 (MQ=255)
                                        cAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaa                                                         <  2:341655/100‑1 (MQ=255)
                                         aaTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaa                                                        <  2:1752030/100‑1 (MQ=255)
                                         aaTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaa                                                        <  1:1039365/100‑1 (MQ=255)
                                         aaTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaa                                                        >  1:679670/1‑100 (MQ=255)
                                          aTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGGCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaa                                                         <  1:1033873/98‑1 (MQ=255)
                                           tGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAac                                                      <  2:3851965/100‑1 (MQ=255)
                                           tGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAac                                                      >  2:3718983/1‑100 (MQ=255)
                                            gCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaaaga                                                     >  1:2586856/1‑98 (MQ=255)
                                            gCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAaca                                                     >  1:1332604/1‑100 (MQ=255)
                                            gCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAaca                                                     >  2:170394/1‑100 (MQ=255)
                                            gCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAaca                                                     >  2:2737482/1‑100 (MQ=255)
                                            gCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAac                                                      <  1:4087728/99‑1 (MQ=255)
                                             ccGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAacac                                                    <  2:1226499/100‑1 (MQ=255)
                                              cgcgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACAcc                                                   <  2:2396461/100‑1 (MQ=255)
                                              cccgTGCGTTGGCAAGAATATACTGCGCCATAAGCACAAAGGTCGGGAAAAGCGGAATCAGACCGGCGATATAATAATTTTTCGTTTTTGCTAAAac                                                      <  1:3620602/95‑1 (MQ=255)
                                                    cGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCTGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCaa                                             >  1:2172373/1‑100 (MQ=255)
                                                    cGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCaa                                             >  2:3898741/1‑100 (MQ=255)
                                                     gCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAAc                                            >  2:3720993/1‑100 (MQ=255)
                                                     gCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAAc                                            <  1:418530/100‑1 (MQ=255)
                                                     gCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAAc                                            <  2:4233080/100‑1 (MQ=255)
                                                      cTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACa                                           >  2:2439834/1‑100 (MQ=255)
                                                      cTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACa                                           <  1:3770477/100‑1 (MQ=255)
                                                        ggCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGc                                         <  2:3324776/100‑1 (MQ=255)
                                                         gCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCa                                        <  1:3134994/100‑1 (MQ=255)
                                                          caacaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTTTTTTTTTTTAAAACACCAATCAACAGCAc                                       >  1:3858738/1‑100 (MQ=255)
                                                          caacaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCAc                                       >  2:1564959/1‑100 (MQ=255)
                                                          caacaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCAc                                       >  2:3821956/1‑100 (MQ=255)
                                                              aaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTAcc                                   <  2:131970/100‑1 (MQ=255)
                                                              aaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTAcc                                   <  2:3168141/100‑1 (MQ=255)
                                                                 ataATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAgc                                >  1:3337060/1‑100 (MQ=255)
                                                                 ataATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAgc                                >  2:127213/1‑100 (MQ=255)
                                                                         cgATAAGCGCAAAGGTCGGAAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTTTTTAAACACCAATCAACAGCACTACCAGCGCACCGa                         >  1:2458794/1‑99 (MQ=255)
                                                                              aGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCgg                   <  2:2826320/100‑1 (MQ=255)
                                                                              aGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCgg                   >  2:1467007/1‑100 (MQ=255)
                                                                              aGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCgg                   <  1:3698284/100‑1 (MQ=255)
                                                                               gcgcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCgg                   >  1:2978092/1‑99 (MQ=255)
                                                                               gcgcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGc                  >  1:2507044/1‑100 (MQ=255)
                                                                               gcgcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGc                  >  2:2410892/1‑100 (MQ=255)
                                                                                  cAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCttt               <  2:3493205/100‑1 (MQ=255)
                                                                                   aaaGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTg              >  2:191728/1‑100 (MQ=255)
                                                                                   aaaGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTg              <  2:531629/100‑1 (MQ=255)
                                                                                     aGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGAt            <  1:917028/100‑1 (MQ=255)
                                                                                     aGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGAt            <  1:1136834/100‑1 (MQ=255)
                                                                                       gTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTa          >  1:2884144/1‑100 (MQ=255)
                                                                                       gTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTa          <  2:3671568/100‑1 (MQ=255)
                                                                                        tCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTAc         >  1:2839837/1‑100 (MQ=255)
                                                                                         cGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTAcc        <  2:2839848/100‑1 (MQ=255)
                                                                                          gggAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCa       <  2:3735568/100‑1 (MQ=255)
                                                                                          gggAAAAGCGGAACCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCa       <  2:3868072/100‑1 (MQ=255)
                                                                                           ggAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAg      >  2:2433518/1‑100 (MQ=255)
                                                                                               aaGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGcccc  <  1:2686687/100‑1 (MQ=255)
                                                                                               aaGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGcccc  >  1:915098/1‑100 (MQ=255)
                                                                                                  |                                                                                                
ATCGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGCCCC  >  NC_000913/1681834‑1682028

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: