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,9320TC96.8% 300.9 / ‑0.9 93P33P (CCA→CCGydgCGlpM family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (43/47);  total (45/48)
Fisher's exact test for biased strand distribution p-value = 6.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

CGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGCC  >  NC_000913/1681836‑1682026
                                                                                                |                                                                                              
cGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGa                                                                                             >  1:233161/1‑100 (MQ=255)
 gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                            <  1:490483/100‑1 (MQ=255)
 gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                            <  1:503821/100‑1 (MQ=255)
 gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                            <  2:1335980/100‑1 (MQ=255)
 gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                            >  2:791942/1‑100 (MQ=255)
     acaTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCa                                                                                         <  2:2863362/99‑1 (MQ=255)
      caTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGc                                                                                       <  1:50670/100‑1 (MQ=255)
       aTACTAAAAATGATGGTTGCGCGTAAGTCTTCAATGCTTCGTTCGCTGGTAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGc                                                                                       <  1:1899233/99‑1 (MQ=255)
        tACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGccc                                                                                     >  1:2046736/1‑100 (MQ=255)
          cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGccc                                                                                     >  2:967468/1‑98 (MQ=255)
          cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCgg                                                                                   <  2:2018282/100‑1 (MQ=255)
          cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCgg                                                                                   >  2:998183/1‑100 (MQ=255)
          cTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCgg                                                                                   <  1:504349/100‑1 (MQ=255)
           tAAAAATGATGGTTGCTCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGc                                                                                  <  1:1210240/100‑1 (MQ=255)
            aaaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGGGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCg                                                                                 <  2:1455947/100‑1 (MQ=255)
            aaaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCg                                                                                 >  2:2362146/1‑100 (MQ=255)
             aaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGa                                                                                <  1:1464897/100‑1 (MQ=255)
              aaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGat                                                                               >  1:1921925/1‑100 (MQ=255)
                atgatgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGatat                                                                             >  2:628964/1‑100 (MQ=255)
                atgatgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGatat                                                                             <  1:205879/100‑1 (MQ=255)
                  gatgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATataa                                                                           >  1:2163938/1‑100 (MQ=255)
                   atgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATataat                                                                          <  1:1596680/100‑1 (MQ=255)
                     ggTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATAtaataa                                                                        >  1:1788049/1‑100 (MQ=255)
                     ggTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATAtaataa                                                                        >  1:2368867/1‑100 (MQ=255)
                      gTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATaataat                                                                       >  1:2865068/1‑100 (MQ=255)
                      gTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATaataat                                                                       >  2:3047388/1‑100 (MQ=255)
                      gTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATaataat                                                                       >  1:2225598/1‑100 (MQ=255)
                        tGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCTGGCGATATAATAAtt                                                                      >  2:1960342/1‑99 (MQ=255)
                          cgcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttttt                                                                   <  2:1591702/100‑1 (MQ=255)
                          cgcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCAGGCGATAAAATAAttttt                                                                   >  1:2809188/1‑100 (MQ=255)
                           gcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTc                                                                  >  1:1453816/1‑100 (MQ=255)
                            cgTAAGGCTTTAATGTCTTGTTCGCTGGTAATAATATAATGCGCGATAAGCGCAAAGGGCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttt                                                                     <  1:1123282/96‑1 (MQ=255)
                            cgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCg                                                                 <  2:2951754/100‑1 (MQ=255)
                               aaGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGttt                                                              <  1:355255/100‑1 (MQ=255)
                                aGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGtttt                                                             >  2:806449/1‑100 (MQ=255)
                                 ggCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGttttt                                                            >  2:2005998/1‑100 (MQ=255)
                                 ggCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGttttt                                                            >  2:2818864/1‑100 (MQ=255)
                                  gCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg                                                           >  1:242705/1‑100 (MQ=255)
                                   cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTATTTGc                                                          >  2:14393/1‑100 (MQ=255)
                                      cAATGCCGCGTTCGCTGGCAACAATATAATGCGCGCTAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaa                                                       >  1:637178/1‑100 (MQ=255)
                                      cAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTCCTaa                                                       >  2:2419561/1‑100 (MQ=255)
                                       aaTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaa                                                      <  2:2616776/100‑1 (MQ=255)
                                       aaTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaa                                                      >  2:2010213/1‑100 (MQ=255)
                                          gCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAaca                                                   >  1:383751/1‑100 (MQ=255)
                                            cgcgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACAcc                                                 >  2:1206758/1‑100 (MQ=255)
                                             gcgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCa                                                >  1:1362237/1‑100 (MQ=255)
                                               gTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAAt                                              <  2:821584/100‑1 (MQ=255)
                                                ttCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATc                                             <  1:1842391/100‑1 (MQ=255)
                                                 tCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCa                                            <  2:2623568/100‑1 (MQ=255)
                                                   gCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAATAGCGGAATCAGCCCGGCGATATATTAATTTTTCGTTTTTGCTAAAACACCAATCAAc                                          <  2:1974406/100‑1 (MQ=255)
                                                   gCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAAc                                          >  1:2251747/1‑100 (MQ=255)
                                                   gCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAAc                                          <  1:1208254/100‑1 (MQ=255)
                                                    cTGGCAATAATATAATTCGCTATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACa                                         <  1:1814967/100‑1 (MQ=255)
                                                    cTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACa                                         >  1:2777190/1‑100 (MQ=255)
                                                     tGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAg                                        <  2:2532058/100‑1 (MQ=255)
                                                      ggCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGc                                       <  2:1858709/100‑1 (MQ=255)
                                                        caacaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCAc                                     >  1:1536872/1‑100 (MQ=255)
                                                         aacaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCAGGCGATATAATAATTTTTCGTTTTTGCTAAAACATCAATCAACAGCAc                                     <  2:2851224/99‑1 (MQ=255)
                                                          acaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGc                                       <  2:300529/96‑1 (MQ=255)
                                                          acaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTa                                   <  2:1210029/100‑1 (MQ=255)
                                                          acaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTa                                   >  1:2434254/1‑100 (MQ=255)
                                                           caaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTAc                                  >  2:2325405/1‑100 (MQ=255)
                                                            aaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTAcc                                 <  1:1241202/100‑1 (MQ=255)
                                                             atatAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCa                                >  1:2447076/1‑100 (MQ=255)
                                                                taATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAgcg                             <  1:2307319/100‑1 (MQ=255)
                                                                taATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAgcg                             >  1:2416442/1‑100 (MQ=255)
                                                                taATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAgcg                             >  1:671316/1‑100 (MQ=255)
                                                                   tGCGCGATAAGCGCAAAGGTCGGGAAACGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCAc                          >  1:283966/1‑100 (MQ=255)
                                                                   tGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCAc                          <  1:2213218/100‑1 (MQ=255)
                                                                   tGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCAc                          <  2:901056/100‑1 (MQ=255)
                                                                        gATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAgc                     <  1:85870/100‑1 (MQ=255)
                                                                            aGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTACAACACCAATCAACAGCACTACCAGCGCACCGAGCGCgg                 >  2:1647337/1‑100 (MQ=255)
                                                                              cgcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCTGCGCACCGAGCGCGGc                <  2:411245/99‑1 (MQ=255)
                                                                              cgcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCt               <  1:806453/100‑1 (MQ=255)
                                                                               gcAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCtt              >  1:1900201/1‑100 (MQ=255)
                                                                                 aaaGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTg            <  1:991100/100‑1 (MQ=255)
                                                                                 aaaGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTg            <  2:1569689/100‑1 (MQ=255)
                                                                                 aaaGGTCGGGAAAAGCGGAATCAGACCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTg            <  1:990267/100‑1 (MQ=255)
                                                                                   aGGTCGGGAAAAGCGGATTTTGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGAt          <  1:1980713/100‑1 (MQ=255)
                                                                                   aGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGAt          <  1:2764607/100‑1 (MQ=255)
                                                                                    ggTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGAtt         >  1:2237154/1‑100 (MQ=255)
                                                                                     gTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACCGCACTACCAGCGCACCGAGCGCGGCTTTGATTa        >  2:2204093/1‑100 (MQ=255)
                                                                                     gTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATAAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTa        >  1:1614012/1‑100 (MQ=255)
                                                                                       ctggAAAAGCGGAATCAGTCCGTCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTATCAGCGCACCGAGCGCGGCTTTGATTAcc      <  1:1346782/98‑1 (MQ=255)
                                                                                       cGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTAcc      <  1:1138357/100‑1 (MQ=255)
                                                                                       cGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTAcc      <  2:2163915/100‑1 (MQ=255)
                                                                                        gggAAAAGCGGAATCAGCCCGGCGATATAATATTTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCa     <  1:557907/100‑1 (MQ=255)
                                                                                        gggAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCa     >  1:291497/1‑100 (MQ=255)
                                                                                         tgAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTAcc      <  2:1024284/97‑1 (MQ=255)
                                                                                         ggAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAg    >  2:1941550/1‑100 (MQ=255)
                                                                                          gAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTTTTATTATTTTAATAAACAGCAATACCACCGCACCGAGCGCGGCATTGATTACCAg    >  2:1500022/1‑99 (MQ=255)
                                                                                           aaaaGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGcc  <  2:1193336/100‑1 (MQ=255)
                                                                                           aaaaGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGcc  <  1:1747117/100‑1 (MQ=255)
                                                                                                |                                                                                              
CGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGCC  >  NC_000913/1681836‑1682026

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: