Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,181,294 G→C 17.8% D292H (GAT→CAT)  ygiC → ATP‑Grasp family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,181,2940GC17.8% 463.3 / 86.4 183D292H (GAT→CAT) ygiCATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand):  ref base G (83/67);  new base C (15/18);  total (98/85)
Fisher's exact test for biased strand distribution p-value = 3.39e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.52e-01

Only 100 of 195 total aligned reads displayed.

GCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGG  >  NC_000913/3181195‑3181393
                                                                                                   |                                                                                                   
 cGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatgn                                                                                                    <  2:4366032/100‑2 (MQ=255)
 cGTGGAAGAGCATTATCTCCAACAAGGCACATCTACCGCTCCTGTGCGAGCTGTTCACGCATCACCCGCACCTGCTGCCCGCTTATTTTGCGGAAgatga                                                                                                    <  1:1298647/100‑1 (MQ=255)
  gTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatgat                                                                                                   <  2:3474718/100‑1 (MQ=255)
   tGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatc                                                                                                  <  1:1266827/100‑1 (MQ=255)
    ggAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatca                                                                                                 >  2:558992/1‑100 (MQ=255)
       agagCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACTCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATcc                                                                                              <  1:637964/100‑1 (MQ=255)
       agagCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATcc                                                                                              >  2:2627634/1‑100 (MQ=255)
       agagCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATcc                                                                                              <  1:2033777/100‑1 (MQ=255)
         agCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGc                                                                                            <  1:4173233/100‑1 (MQ=255)
         agCATTATCTCCAACAAGGCACTTCTACCGCTACTGGGGGAGATGTTGCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGc                                                                                            <  1:3021163/100‑1 (MQ=255)
          gCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCa                                                                                           >  2:3277957/1‑100 (MQ=255)
           cATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCaa                                                                                          >  2:507547/1‑100 (MQ=255)
            aTTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaa                                                                                         <  2:920357/100‑1 (MQ=255)
            aTTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCaaa                                                                                         >  2:996330/1‑100 (MQ=255)
             ttATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaa                                                                                        <  2:1450886/100‑2 (MQ=39)
              tATCTCCAACACGGCACTTCTCCCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATg                                                                                       <  2:3506912/100‑1 (MQ=255)
              tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaat                                                                                       <  1:1260299/100‑3 (MQ=38)
              tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATg                                                                                       <  1:1606144/100‑1 (MQ=255)
              tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATg                                                                                       <  2:191372/100‑1 (MQ=255)
              tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATg                                                                                       <  2:1592843/100‑1 (MQ=255)
               atttcCAACAAGGCACGTGTACAGCGATTGTGGGTGATGTTCTCGAAGCACCCGAACCTGCTGCTAGCTTATTTTGCGGAAGATGATCATCCGCAAATgg                                                                                      >  1:918701/4‑100 (MQ=255)
               aTCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaata                                                                                      >  1:1564316/1‑97 (MQ=38)
               aTCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaat                                                                                       <  2:1564321/99‑3 (MQ=39)
                tctcCAACAAGGCACTTCTACCCCTACTGTTGGAGATTTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaat                                                                                       <  2:2937762/98‑3 (MQ=37)
                 ctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataag                                                                                    <  1:2987924/100‑6 (MQ=38)
                 ctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaa                                                                                    >  1:2839504/1‑100 (MQ=255)
                 ctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaa                                                                                    >  1:751669/1‑100 (MQ=255)
                 ctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGAGCATCCGCAAATGGaa                                                                                    >  2:4333133/1‑100 (MQ=255)
                  tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagc                                                                                   <  1:1511769/100‑7 (MQ=37)
                  tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagc                                                                                   <  2:250571/100‑7 (MQ=37)
                  tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagc                                                                                   <  1:4383467/100‑7 (MQ=38)
                  tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaa                                                                                   >  2:2359781/1‑100 (MQ=255)
                  tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaa                                                                                   >  2:2281433/1‑100 (MQ=255)
                   ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagcg                                                                                  <  1:2458624/100‑8 (MQ=37)
                   ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagcg                                                                                  <  1:4263219/100‑8 (MQ=37)
                   ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaaa                                                                                  >  2:2776575/1‑100 (MQ=255)
                   ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCACAAATGGaaaa                                                                                  <  1:2168117/100‑1 (MQ=255)
                    caacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagcgg                                                                                 <  1:3680921/100‑9 (MQ=37)
                        aaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtg                                                                             >  1:4089155/1‑100 (MQ=255)
                        aaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtg                                                                             <  2:764409/100‑1 (MQ=255)
                         aGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtgt                                                                            <  2:2842831/100‑1 (MQ=255)
                            cACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGt                                                                         >  1:1306561/1‑100 (MQ=255)
                             aCTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATTTGTGGtt                                                                        <  2:2309402/100‑1 (MQ=255)
                              cTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTa                                                                       >  1:719430/1‑100 (MQ=255)
                              cTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTa                                                                       <  1:3219340/100‑1 (MQ=255)
                                 cTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAc                                                                    <  1:935691/100‑1 (MQ=255)
                                  tACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                                   >  1:425137/1‑100 (MQ=255)
                                  tACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                                   <  2:2667875/100‑1 (MQ=255)
                                  tACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                                   <  1:566597/100‑1 (MQ=255)
                                   aCCGCTACTGTGGGAGAAGTTACCGAATCACCCGAACCTCCTGCCCGCTTATTTTGTGAAAGACGATTATCCGCAAATGGAAAAATATGGGGTTCAAcc                                                                   >  1:1998395/1‑99 (MQ=255)
                                    ccGCTACTGTGGGAGATGATAACGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAAATGATCATCAGCAAATGGAAAAATATGTGGTTAAACCg                                                                  >  1:256126/1‑99 (MQ=255)
                                         aCTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtcttc                                                            <  2:948356/100‑1 (MQ=255)
                                            gtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTcc                                                         >  1:784158/1‑100 (MQ=255)
                                            gtgGGAGATGTGCCTGAATCACCCGTTCCTGCTGCACGGTTATTTTGCGGAAGATGATCATCCGCAAATGGAAATATATGTGGTTATGCCGAtcttct                                                           <  1:345657/98‑1 (MQ=255)
                                              gggAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATATTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCg                                                       >  1:2403611/1‑100 (MQ=255)
                                               ggAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGt                                                      <  2:1987291/100‑1 (MQ=255)
                                                    tGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                                 >  1:4138384/1‑100 (MQ=255)
                                                     gTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgc                                                >  1:2285686/1‑100 (MQ=255)
                                                       tCCCGAATCACCCGAAAATGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcg                                               >  1:1863047/1‑99 (MQ=255)
                                                          cGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCaaa                                           >  1:2664713/1‑100 (MQ=255)
                                                            aaTCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACg                                         <  2:3607884/100‑1 (MQ=255)
                                                                aCCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgtgt                                      >  1:1924597/1‑99 (MQ=255)
                                                                 cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCg                                    >  1:3501117/1‑100 (MQ=255)
                                                                 cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCGTCTCCCGTGAAGGCGCAAACGTGTCg                                    <  2:2679323/100‑1 (MQ=255)
                                                                  ncGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGa                                   >  2:1780268/2‑100 (MQ=255)
                                                                        ctgctgCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTg                             <  2:3741451/100‑1 (MQ=255)
                                                                         tgctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTga                            >  1:1675311/1‑100 (MQ=255)
                                                                         tgctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTg                             <  2:3927033/99‑1 (MQ=255)
                                                                          gctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTgag                           <  2:1206371/100‑1 (MQ=255)
                                                                          gctgcCCGCTTATTTTGCGGAAGATGACCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTgag                           <  1:2390486/100‑1 (MQ=255)
                                                                            tgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGaa                         >  1:1763333/1‑100 (MQ=255)
                                                                               ccGCTTATTTTGCGGAAGATGATCATCCGCAAATGTAAAAATAGGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACg                       <  2:3096310/99‑1 (MQ=255)
                                                                                cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                     >  2:1218644/1‑100 (MQ=255)
                                                                                cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                     <  2:2883613/100‑1 (MQ=255)
                                                                                cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                     >  2:378898/1‑100 (MQ=255)
                                                                                 gCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCa                    >  1:3314262/1‑100 (MQ=255)
                                                                                   ttATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCaaa                  >  2:4247186/1‑100 (MQ=255)
                                                                                   tccatttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCaaa                  >  2:1263547/5‑100 (MQ=38)
                                                                                   tccatttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCaaa                  >  2:3129084/5‑100 (MQ=38)
                                                                                     catttGCGGATGTTCATCTTCCGCAACTGGAAAAATATGTGGTTAAACCGATCCTCTCCCGTGAACGCGCACACGTGTCGATCATTGCGAACGGCAAAAc                >  1:3968023/3‑100 (MQ=255)
                                                                                     catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc                >  2:1068791/3‑100 (MQ=39)
                                                                                     catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc                >  1:824238/3‑100 (MQ=39)
                                                                                     catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc                >  1:3476987/3‑100 (MQ=39)
                                                                                     aTTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc                <  1:275181/100‑1 (MQ=255)
                                                                                     aTTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc                >  2:1690996/1‑100 (MQ=255)
                                                                                      ntttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAcc               >  1:2438468/2‑100 (MQ=255)
                                                                                      atttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAcc               >  2:2293952/2‑100 (MQ=39)
                                                                                      atttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAAAGTGTCGATCATTGAGAACGGCAAAAc                >  1:3367453/2‑99 (MQ=39)
                                                                                       tttGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCa              <  2:814176/100‑1 (MQ=255)
                                                                                        ttGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAt             <  2:852549/100‑1 (MQ=255)
                                                                                         tGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAtt            >  2:409610/1‑100 (MQ=255)
                                                                                         tGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAtt            <  1:1380480/100‑1 (MQ=255)
                                                                                          gCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTg           <  2:3314278/100‑1 (MQ=255)
                                                                                          gCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTg           <  1:1706779/100‑1 (MQ=255)
                                                                                           cGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGa          >  2:1761770/1‑100 (MQ=255)
                                                                                            ggAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGa          >  1:3144621/1‑99 (MQ=255)
                                                                                             gAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAag        >  2:742993/1‑100 (MQ=255)
                                                                                             gAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAag        <  1:3862128/100‑1 (MQ=255)
                                                                                                 atgatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAagcagc    >  2:4213852/1‑100 (MQ=255)
                                                                                                  tgatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGACGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCg   >  1:3756846/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGG  >  NC_000913/3181195‑3181393

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