Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 354,812 C→T intergenic (+220/‑110) prpE → / → codB propionate‑‑CoA ligase/cytosine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913354,8120CT98.2% 555.2 / ‑0.6 163intergenic (+220/‑110)prpE/codBpropionate‑‑CoA ligase/cytosine transporter
Reads supporting (aligned to +/- strand):  ref base C (1/2);  new base T (77/83);  total (78/85)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.28e-01

Only 100 of 170 total aligned reads displayed.

CGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAA  >  NC_000913/354717‑354907
                                                                                               |                                                                                               
cGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATtcat                                                                                             >  1:313252/1‑100 (MQ=255)
 gACACCATGCTCAGATGCCTGATGCGCCGCTGATGCGTCTTATCACGCCTACCCACTGTTTTTACAGCAATAATTTTTCCCCCACCTTTTTGcattcatt                                                                                            <  2:3715745/100‑1 (MQ=37)
  acacCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCattcattc                                                                                           <  1:515492/100‑1 (MQ=255)
   cacCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCAttcattca                                                                                          <  2:2824913/100‑1 (MQ=255)
     ccATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCata                                                                                        <  2:3295561/100‑1 (MQ=255)
     ccATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCata                                                                                        >  1:1295480/1‑100 (MQ=255)
     ccATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCata                                                                                        >  1:4017122/1‑100 (MQ=255)
     ccATGCTCAGAGGCCTGATGCGGCGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCata                                                                                        <  2:3894156/100‑1 (MQ=38)
        tGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATaa                                                                                     >  1:3309180/1‑100 (MQ=255)
          cTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATaaaa                                                                                   >  2:356992/1‑100 (MQ=255)
             aGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAat                                                                                <  2:652339/100‑1 (MQ=255)
             aGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAat                                                                                <  2:3061852/100‑1 (MQ=255)
              gATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAata                                                                               >  1:492959/1‑100 (MQ=255)
              gATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAata                                                                               <  1:2551134/100‑1 (MQ=255)
              gATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATGCATTCATATGAAAAata                                                                               >  1:4204239/1‑100 (MQ=38)
                 gCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAatatat                                                                            <  1:1066556/100‑1 (MQ=255)
                  ccTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCTCCTTTTTGCATTCATTCATATAAAAAATATAtt                                                                           <  2:1089265/100‑1 (MQ=255)
                  ccTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATAtt                                                                           >  1:2207091/1‑100 (MQ=255)
                  ccTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATAtt                                                                           >  2:2437854/1‑100 (MQ=255)
                   cTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATAttt                                                                          >  1:668812/1‑100 (MQ=255)
                    tGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTc                                                                         >  1:1871280/1‑100 (MQ=255)
                      aTGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTccc                                                                       >  2:1191929/1‑100 (MQ=255)
                      aTGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTccc                                                                       <  2:485075/100‑1 (MQ=255)
                       tGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTcccc                                                                      >  1:1383150/1‑100 (MQ=255)
                       tGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTAATTAATATAAAAAATATATTTcccc                                                                      <  2:306361/100‑1 (MQ=255)
                        gCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTccccn                                                                     <  2:1042146/100‑2 (MQ=255)
                               tGACGCGTCTTATCAGGCCTACCCGCTGTTTTTACACGGATAATTTTTCCCCCACTTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAAc                                                              <  2:1547115/100‑1 (MQ=255)
                               tGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAAc                                                              <  1:550141/100‑1 (MQ=255)
                                gACGCGTCTTATCAGGCCTACCCACTTTTTTTACACCGATAATTTTTCCCCCCCTTTTTTGCATTCATTCAAATAAAGAATATATTTCCCCACGAAAGCg                                                             >  1:4154772/1‑100 (MQ=255)
                                gACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACg                                                             <  1:4345956/100‑1 (MQ=255)
                                gACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTCCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACg                                                             <  2:2502932/100‑1 (MQ=255)
                                  cgcgTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGAt                                                           <  1:2773833/100‑1 (MQ=255)
                                   gcgTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTTTTTTTTTTTTTATAAAAAATATATTTCCCCACGAAAACGAtt                                                          >  2:2414181/1‑100 (MQ=255)
                                    cgTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTg                                                         <  2:3222504/100‑1 (MQ=255)
                                    cgTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTg                                                         <  1:3978689/100‑1 (MQ=255)
                                    cgTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTg                                                         <  2:3079808/100‑1 (MQ=255)
                                     gTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGc                                                        >  1:822241/1‑100 (MQ=255)
                                        ttATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCttt                                                     <  2:1892033/100‑1 (MQ=255)
                                          aTCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCttttt                                                   >  1:3839970/1‑100 (MQ=255)
                                            cAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTAt                                                 >  2:447390/1‑100 (MQ=255)
                                            cAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACATTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTAt                                                 >  2:561949/1‑100 (MQ=255)
                                             aGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATc                                                <  1:3206748/100‑1 (MQ=255)
                                              ggCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCt                                               >  1:3706250/1‑100 (MQ=255)
                                              ggCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCt                                               >  2:3674883/1‑100 (MQ=255)
                                                ncTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTc                                             >  1:1929338/2‑100 (MQ=255)
                                                ccTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTc                                             >  2:1874496/1‑100 (MQ=255)
                                                ccTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTc                                             >  1:2077053/1‑100 (MQ=255)
                                                ccTACCCACTGGTTTTACACCGATGATTTTCCCCCCACCTTTTTGCATTCATTCATATAATAAATATATCTCCCCACGAAAACGATTGCTTTTTATCTTc                                             <  2:3660548/100‑1 (MQ=255)
                                                 cTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTTCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCa                                            <  1:849141/100‑1 (MQ=255)
                                                 cTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCa                                            >  2:4198202/1‑100 (MQ=255)
                                                  tACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAg                                           <  1:1751452/100‑1 (MQ=255)
                                                  tACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAg                                           <  2:4124165/100‑1 (MQ=255)
                                                   aCCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGa                                          >  2:2868229/1‑100 (MQ=255)
                                                     ccACTGTTTTTACACCGATAATTTTTCCCCCGCCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGAtn                                        <  1:3635342/100‑2 (MQ=255)
                                                      cACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGa                                       >  2:261692/1‑100 (MQ=255)
                                                       aCTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGaa                                      <  2:86825/100‑1 (MQ=255)
                                                         tGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATa                                    <  1:2055566/100‑1 (MQ=255)
                                                         tGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATa                                    >  1:2345482/1‑100 (MQ=255)
                                                          gTTTTTACACCGATAATTTTTCCCCCACCTTTTTTTTTTTTTTTTTATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTGTCTTCAGATGAAttt                                   >  2:2120029/1‑98 (MQ=255)
                                                          gTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAg                                   >  1:933178/1‑100 (MQ=255)
                                                             tttACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGATTTTTATCTTCAGATGAATAGAAt                                >  2:3264908/1‑100 (MQ=255)
                                                               tACACCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATgc                              <  2:2917074/100‑1 (MQ=255)
                                                                 cacCTATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATgcgg                            <  2:1201280/100‑1 (MQ=255)
                                                                 cacCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATgcgg                            <  1:1837490/100‑1 (MQ=255)
                                                                 cacCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATgcgg                            >  1:3491378/1‑100 (MQ=255)
                                                                  acCGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATgcggc                           <  2:1387848/100‑1 (MQ=255)
                                                                   ccGATAATTTTTCCCCCACCTTTTTTTTTTTTTTTTTATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGTATgcggcg                          >  2:700792/1‑100 (MQ=255)
                                                                   ccGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATgtggcg                          >  1:1074958/1‑100 (MQ=255)
                                                                   ccGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATgcggcg                          <  1:4313570/100‑1 (MQ=255)
                                                                   ccGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATgcggcg                          >  1:1672403/1‑100 (MQ=255)
                                                                    cGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGcggcgg                         >  2:1329570/1‑100 (MQ=255)
                                                                    cGATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGcggcgg                         >  1:2833954/1‑100 (MQ=255)
                                                                     gATAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGa                        >  1:1882176/1‑100 (MQ=255)
                                                                       tAATTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGAtt                      >  2:3679031/1‑100 (MQ=255)
                                                                         aTTTTTCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGAtttt                    >  2:656947/1‑100 (MQ=255)
                                                                          tttttCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGAttttt                   >  1:3417423/1‑100 (MQ=255)
                                                                           ttttCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGAtttttt                  <  1:1228715/100‑1 (MQ=255)
                                                                             ttCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTgg                <  1:2750584/100‑1 (MQ=255)
                                                                              tCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTggg               <  1:4382579/100‑1 (MQ=255)
                                                                              tCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTggg               >  2:3568987/1‑100 (MQ=255)
                                                                              tCCCCCACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTggg               <  1:4466408/100‑1 (MQ=255)
                                                                                ccccACCTTTTTGCATTGACTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGtt             <  1:2714158/100‑1 (MQ=255)
                                                                                ccccACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGtt             >  1:556577/1‑100 (MQ=255)
                                                                                 cccACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGttt            <  2:273186/100‑1 (MQ=255)
                                                                                  ccACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTc           >  1:4397244/1‑100 (MQ=255)
                                                                                   cACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCa          <  1:470682/100‑1 (MQ=255)
                                                                                   cACCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCa          <  1:1834909/100‑1 (MQ=255)
                                                                                    aCCTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCaa         <  2:1112842/100‑1 (MQ=255)
                                                                                      cTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAAc       <  1:2683597/100‑1 (MQ=255)
                                                                                      cTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAAc       <  2:3262019/100‑1 (MQ=255)
                                                                                      cTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAAc       <  1:3948340/100‑1 (MQ=255)
                                                                                      cTTTTTGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAAc       <  2:4259779/100‑1 (MQ=255)
                                                                                      cTTTTTGCATTCAATCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAAc       <  2:1694729/100‑1 (MQ=255)
                                                                                       tttttGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACa      <  2:1564125/100‑1 (MQ=255)
                                                                                        ttttGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAg     <  1:1412939/100‑1 (MQ=255)
                                                                                        ttttGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAg     <  1:3225214/100‑1 (MQ=255)
                                                                                         tttGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGc    <  2:2713088/100‑1 (MQ=255)
                                                                                           tGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCaa  <  1:3293752/100‑1 (MQ=255)
                                                                                           tGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCaa  >  1:725967/1‑100 (MQ=255)
                                                                                           tGCATTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCaa  >  2:3070612/1‑100 (MQ=255)
                                                                                               |                                                                                               
CGACACCATGCTCAGATGCCTGATGCGACGCTGACGCGTCTTATCAGGCCTACCCACTGTTTTTACACCGATAATTTTTCCCCCACCTTTTTGCACTCATTCATATAAAAAATATATTTCCCCACGAAAACGATTGCTTTTTATCTTCAGATGAATAGAATGCGGCGGATTTTTTGGGTTTCAAACAGCAA  >  NC_000913/354717‑354907

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: