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,9320TC98.3% 388.3 / ‑3.3 113P33P (CCA→CCG) ydgCGlpM family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base C (55/56);  total (57/56)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 113 total aligned reads displayed.

CGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGCCC  >  NC_000913/1681836‑1682027
                                                                                                |                                                                                               
cGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGAAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGa                                                                                              <  1:1448963/100‑1 (MQ=255)
 gACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGaa                                                                                             <  2:1886214/100‑1 (MQ=255)
   ccACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATc                                                                                           >  2:3383988/1‑100 (MQ=255)
      caTACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGc                                                                                        >  2:834862/1‑100 (MQ=255)
         aCTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCg                                                                                     >  2:1470761/1‑100 (MQ=255)
          cTAAAAATGATGGTTGCGCGTAAGGCTTCAATTCCGCGTTCGCTGGCAACAATATAATGCGCTATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCg                                                                                     >  1:2278268/1‑99 (MQ=255)
           tAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGc                                                                                   <  1:566094/100‑1 (MQ=255)
           tAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGc                                                                                   <  2:712573/100‑1 (MQ=255)
             aaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGa                                                                                 >  1:1049767/1‑100 (MQ=255)
             aaaaTGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGAAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGa                                                                                 >  2:2809845/1‑100 (MQ=255)
                atgatgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGatat                                                                              <  1:253753/100‑1 (MQ=255)
                atgatgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGatat                                                                              <  2:2693184/100‑1 (MQ=255)
                 tgatgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGAtata                                                                             <  2:558719/100‑1 (MQ=255)
                 tgatgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGAtata                                                                             <  1:1896120/100‑1 (MQ=255)
                   atgGTTGCTCGTAAGGCTTCGAGGCCGCGTTCGCTGGAAACAAGATAATGCGCGAGAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGCTataa                                                                            <  1:2326243/99‑1 (MQ=255)
                    tgGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATataata                                                                          >  1:1216665/1‑100 (MQ=255)
                     ggTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATAtaataa                                                                         >  2:1597149/1‑100 (MQ=255)
                     ggTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATAtaataa                                                                         >  2:2810315/1‑100 (MQ=255)
                     ggTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATAtaataa                                                                         <  2:2716993/100‑1 (MQ=255)
                      gTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATaataat                                                                        >  2:1269321/1‑100 (MQ=255)
                        taagcgTAAGGCTTCAATGCCTCTTTCGCTAGAAACAATATAATGCACGATAAGCGCAAAGTTTGTGAAAAGCGGAATCAGCCCAGCGATATAATAAttt                                                                      <  1:1099386/97‑1 (MQ=255)
                        tGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttt                                                                      >  2:1156594/1‑100 (MQ=255)
                        tGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttt                                                                      >  2:1236307/1‑100 (MQ=255)
                        tGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttt                                                                      <  2:2484147/100‑1 (MQ=255)
                         gcgcGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAAttt                                                                      <  2:2763600/99‑1 (MQ=255)
                         gcgcGTAAGGCTTCAATGCCGCGTTCGCTGGAAACAATATAAAGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATaataat                                                                        <  1:224546/97‑1 (MQ=255)
                           gcgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTc                                                                   <  2:302433/100‑1 (MQ=255)
                           gcgTAAGGCTACAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTc                                                                   <  1:2226315/100‑1 (MQ=255)
                            cgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTcn                                                                  <  2:1424436/100‑2 (MQ=255)
                            cgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCg                                                                  <  2:1625948/100‑1 (MQ=255)
                            cgTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCg                                                                  <  1:2580867/100‑1 (MQ=255)
                               aaGGCTTCAATGCCGCGTTCGCTGACAACAATACAAAGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGttt                                                               <  1:374462/100‑1 (MQ=255)
                                 ggCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGttttt                                                             >  2:1272396/1‑100 (MQ=255)
                                  gCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg                                                            <  1:3163730/100‑1 (MQ=255)
                                  gCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg                                                            <  2:2614792/100‑1 (MQ=255)
                                  gCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTg                                                            >  2:1658014/1‑100 (MQ=255)
                                   gatctaTGCCGCGTTGGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                           <  2:3333113/95‑1 (MQ=255)
                                   cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGCAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                           >  2:410551/1‑100 (MQ=255)
                                   cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTcagatcgga                                                           >  1:545279/1‑92 (MQ=255)
                                   cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                           >  2:2489497/1‑100 (MQ=255)
                                   cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                           <  2:1542750/100‑1 (MQ=255)
                                   cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                           <  2:1114443/100‑1 (MQ=255)
                                   cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                           >  2:129375/1‑100 (MQ=255)
                                   cTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGc                                                           <  2:3128816/100‑1 (MQ=255)
                                     tCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTa                                                         >  2:2670922/1‑100 (MQ=255)
                                     tCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTa                                                         <  2:2479939/100‑1 (MQ=255)
                                      cAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaa                                                        >  2:2853283/1‑100 (MQ=255)
                                       aaTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaa                                                       <  1:1241719/100‑1 (MQ=255)
                                       aaTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaa                                                       <  1:650937/100‑1 (MQ=255)
                                        aTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGcagatc                                                      >  1:3333085/1‑95 (MQ=255)
                                        aTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaaa                                                      <  1:2705179/100‑1 (MQ=255)
                                        aTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTaaaa                                                      >  2:3379545/1‑100 (MQ=255)
                                        aTGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAAACAGCGCGGCGATATAATAATATTTCGTTTTTGATaaaa                                                      >  1:1379954/1‑100 (MQ=255)
                                         tGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAac                                                     <  1:1101953/100‑1 (MQ=255)
                                         tGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAac                                                     <  2:1398736/100‑1 (MQ=255)
                                         tGCCGCGTTCGCGGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAac                                                     <  1:427/100‑1 (MQ=255)
                                          gCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAaca                                                    >  1:692958/1‑100 (MQ=255)
                                          gCCGCGTTCGATGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAac                                                     <  2:813672/99‑1 (MQ=255)
                                            cgcgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAACACAcc                                                  >  1:2972787/1‑100 (MQ=255)
                                             gcgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCa                                                 >  1:545591/1‑100 (MQ=255)
                                              cgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATACTTTTTCGTTTTTGCTAAAACACCaa                                                <  1:671987/100‑1 (MQ=255)
                                              cgTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCaa                                                <  2:1867681/100‑1 (MQ=255)
                                                ttCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATc                                              >  2:116164/1‑100 (MQ=255)
                                                  cGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCa                                             >  1:3146824/1‑99 (MQ=255)
                                                  cGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCaa                                            >  2:1296114/1‑100 (MQ=255)
                                                  cGCTGGCAACAATATAATGCGCGATAAGAGCAAAGGACGGGAAAACCGGCATCAGCCCGGCGACATCATAATTTTACGTTTTTGCTACAccacca                                                 >  1:2950008/1‑95 (MQ=255)
                                                    cTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACa                                          >  2:3151981/1‑100 (MQ=255)
                                                        caacaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTTTTTTTTTTTTTAACGCCAATAAACAGCa                                       >  1:889845/1‑99 (MQ=255)
                                                        caacaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCAc                                      >  2:2489096/1‑100 (MQ=255)
                                                        caacaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCAc                                      >  1:2198028/1‑100 (MQ=255)
                                                          acaaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACAAATCAACAGCACTa                                    >  2:2689025/1‑100 (MQ=255)
                                                           caaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTAc                                   >  2:1698844/1‑100 (MQ=255)
                                                            aaTATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTAcc                                  <  2:732350/100‑1 (MQ=255)
                                                                 aaTGCGCGATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAgcgc                             >  2:2739250/1‑100 (MQ=255)
                                                                     cgcgATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCg                         >  2:113253/1‑100 (MQ=255)
                                                                     cgcgATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCg                         >  2:3213381/1‑100 (MQ=255)
                                                                       cgATAAGCGCAAAGGTCGGGAACAGCGGAATCAGCACGGCGGTATAATAATGTATCGTTTTTGCGAAAACACGGATCGACAGCACTACTAGCacaca                          >  1:819945/1‑96 (MQ=255)
                                                                       cgATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAg                       >  2:712579/1‑100 (MQ=255)
                                                                        gATAAGCGCAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAgc                      <  1:343273/100‑1 (MQ=255)
                                                                         aTAAGCGCAAAGGTCTGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCa                            <  1:270010/93‑1 (MQ=255)
                                                                             gcgcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGTTAAAACACCAATCAACAGCACTACCAGCCCACCGAGcgcgcc                 <  2:2293238/100‑3 (MQ=255)
                                                                             gcgcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGc                 >  2:1341407/1‑100 (MQ=255)
                                                                             gcgcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGc                 >  2:1135503/1‑100 (MQ=255)
                                                                              cgcAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCt                >  2:1060739/1‑100 (MQ=255)
                                                                                cAAAGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCttt              >  2:3003157/1‑100 (MQ=255)
                                                                                 aaaGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTg             <  1:2608347/100‑1 (MQ=255)
                                                                                  aaGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGa            <  2:881583/100‑1 (MQ=255)
                                                                                   aGGTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGAt           <  2:1613438/100‑1 (MQ=255)
                                                                                    ggTCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGAtt          >  2:3352668/1‑100 (MQ=255)
                                                                                     gTCAGGAAAACCCGAATCAGCCCGGCGCTATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTa         <  1:695564/100‑1 (MQ=255)
                                                                                      tCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTTCTTTTTCTTCAATCAACAGCACTACTAACGCACCGCGCCCCGCTTTGATTAc        >  1:2087215/1‑100 (MQ=255)
                                                                                      tCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTAc        >  2:2401915/1‑100 (MQ=255)
                                                                                      tCGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTAc        >  1:2087135/1‑100 (MQ=255)
                                                                                       cGGGAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTAcc       >  2:450479/1‑100 (MQ=255)
                                                                                        gggAAAAGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCa      <  2:1756738/100‑1 (MQ=255)
                                                                                           aaaaGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGcc   <  1:3180265/100‑1 (MQ=255)
                                                                                            aaaGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGccc  <  2:3196817/100‑1 (MQ=255)
                                                                                            aaaGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGccc  <  2:2824996/100‑1 (MQ=255)
                                                                                            aaaGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGccc  <  1:1401829/100‑1 (MQ=255)
                                                                                             aaGCGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTTTTTTTTTTTTAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAg     >  1:3031286/1‑96 (MQ=255)
                                                                                                |                                                                                               
CGACCACATACTAAAAATGATGGTTGCGCGTAAGGCTTCAATGCCGCGTTCGCTGGCAACAATATAATGCGCGATAAGCGCAAAGGTCGGGAAAAGTGGAATCAGCCCGGCGATATAATAATTTTTCGTTTTTGCTAAAACACCAATCAACAGCACTACCAGCGCACCGAGCGCGGCTTTGATTACCAGCCC  >  NC_000913/1681836‑1682027

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: ‑