Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,181,290 A→T 16.5% E290D (GAA→GAT)  ygiC → ATP‑Grasp family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,181,2900AT16.5% 476.0 / 77.0 193E290D (GAA→GAT) ygiCATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand):  ref base A (87/73);  new base T (14/18);  total (102/91)
Fisher's exact test for biased strand distribution p-value = 3.33e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

Only 100 of 196 total aligned reads displayed.

ACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGC  >  NC_000913/3181191‑3181388
                                                                                                   |                                                                                                  
aCCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGaa                                                                                                    >  2:105131/1‑100 (MQ=255)
 ccGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAg                                                                                                   >  2:4250273/1‑100 (MQ=255)
   ggCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgat                                                                                                 <  1:3162965/100‑1 (MQ=255)
   ggCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgat                                                                                                 >  2:1329897/1‑100 (MQ=255)
   ggCGTGGAAGAGCATTATATCCAACAAGGCACTTCTGCCGCTAGTGTGGGAGATGTTCACGAATCACCCGAGCCTGCTGCGCGCTTATTTTGCGGAAgat                                                                                                 <  1:1808088/100‑1 (MQ=255)
     cGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatgn                                                                                               <  2:4366032/100‑2 (MQ=255)
     cGTGGAAGAGCATTATCTCCAACAAGGCACATCTACCGCTCCTGTGCGAGCTGTTCACGCATCACCCGCACCTGCTGCCCGCTTATTTTGCGGAAgatga                                                                                               <  1:1298647/100‑1 (MQ=255)
       tGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatc                                                                                             <  1:1266827/100‑1 (MQ=255)
        ggAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGCGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatca                                                                                            >  1:2007313/1‑100 (MQ=255)
        ggAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGCATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatca                                                                                            >  1:2820155/1‑100 (MQ=255)
        ggAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatca                                                                                            >  2:558992/1‑100 (MQ=255)
           agagCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGATATGTTCCCGAATCACCCGAACCTCCTCCCCGCCTCTTTTGCTGAAAATGagca                                                                                            >  1:1895573/1‑97 (MQ=255)
           agagCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACTCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATcc                                                                                         <  1:637964/100‑1 (MQ=255)
           agagCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATcc                                                                                         >  2:2627634/1‑100 (MQ=255)
             agCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGc                                                                                       <  1:4173233/100‑1 (MQ=255)
             agCATTATCTCCAACAAGGCACTTCTACCGCTACTGGGGGAGATGTTGCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGc                                                                                       <  1:3021163/100‑1 (MQ=255)
             agCAATATCTCCAACATGGCCCTGCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGATCCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGc                                                                                       <  1:1992650/100‑1 (MQ=255)
               cATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCaa                                                                                     >  2:507547/1‑100 (MQ=255)
                aTTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaa                                                                                    <  2:920357/100‑1 (MQ=255)
                aTTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCaaa                                                                                    <  2:3262912/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                                                                                  <  2:1592843/100‑1 (MQ=255)
                  tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATg                                                                                  <  1:1606144/100‑1 (MQ=255)
                   aTCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaata                                                                                 >  1:1564316/1‑97 (MQ=38)
                   aTCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaat                                                                                  <  2:1564321/99‑3 (MQ=39)
                    tctcCAACAAGGCACTTCTACCCCTACTGTTGGAGATTTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaat                                                                                  <  2:2937762/98‑3 (MQ=37)
                     ctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaa                                                                               >  1:751669/1‑100 (MQ=255)
                     ctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaa                                                                               >  1:2369675/1‑100 (MQ=255)
                     ctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGAGCATCCGCAAATGGaa                                                                               >  2:4333133/1‑100 (MQ=255)
                      tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaagtaagc                                                                              <  1:3276002/100‑7 (MQ=38)
                      tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagc                                                                              <  1:1511769/100‑7 (MQ=37)
                      tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagc                                                                              <  1:4383467/100‑7 (MQ=38)
                      tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagc                                                                              <  2:250571/100‑7 (MQ=37)
                      tcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaa                                                                              >  2:2359781/1‑100 (MQ=255)
                       ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagcg                                                                             <  1:4263219/100‑8 (MQ=37)
                       ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGaaaa                                                                             >  2:2776575/1‑100 (MQ=255)
                        caacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagcgg                                                                            <  1:3680921/100‑9 (MQ=37)
                        caacaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagcgg                                                                            <  2:1984636/100‑9 (MQ=37)
                          acaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAata                                                                          <  1:3827690/100‑1 (MQ=255)
                            aaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtg                                                                        >  1:4089155/1‑100 (MQ=255)
                            aaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtg                                                                        <  2:764409/100‑1 (MQ=255)
                             aGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtgt                                                                       <  2:2842831/100‑1 (MQ=255)
                                cACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGt                                                                    >  2:3107750/1‑100 (MQ=255)
                                cACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGt                                                                    >  1:3957736/1‑100 (MQ=255)
                                cACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGt                                                                    >  1:1306561/1‑100 (MQ=255)
                                 aCTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATTTGTGGtt                                                                   <  2:2309402/100‑1 (MQ=255)
                                 aCTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCAACAACAAATAGAAAAATAtgtg                                                                      <  1:3332189/97‑1 (MQ=255)
                                  cTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTa                                                                  >  1:719430/1‑100 (MQ=255)
                                    tCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTaaa                                                                >  1:2234253/1‑100 (MQ=255)
                                     cTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAc                                                               <  1:935691/100‑1 (MQ=255)
                                      tACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                              <  1:566597/100‑1 (MQ=255)
                                      tACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                              >  1:425137/1‑100 (MQ=255)
                                      tACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                              <  2:2667875/100‑1 (MQ=255)
                                        ccGCTACTGTGGGAGATGATAACGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAAATGATCATCAGCAAATGGAAAAATATGTGGTTAAACCg                                                             >  1:256126/1‑99 (MQ=255)
                                          gCTACTGTGGGAGATGTTCCCGAATCAGCCCAACCTGCTGGCCGCTTATTTTGCGGAAGATGATCATCAGCAAAGGGAAACATATGTGGCTAAACCGAtc                                                          <  1:821017/100‑1 (MQ=255)
                                             aCTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtcttc                                                       <  2:948356/100‑1 (MQ=255)
                                                   ggAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGt                                                 >  1:2034102/1‑100 (MQ=255)
                                                   ggAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTATGCGGAAGATGATCATCCGAAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGt                                                 >  1:2398188/1‑100 (MQ=255)
                                                    gagaTGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTg                                                <  2:122840/100‑1 (MQ=255)
                                                     agaTGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGa                                               >  1:3402823/1‑100 (MQ=255)
                                                     agaTGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGa                                               <  1:2440218/100‑1 (MQ=255)
                                                        tGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                            >  2:3302916/1‑100 (MQ=255)
                                                        tGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                            >  1:4138384/1‑100 (MQ=255)
                                                         gTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATCGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAgtc                                           <  2:2969960/100‑3 (MQ=255)
                                                              cGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCaaa                                      >  1:2664713/1‑100 (MQ=255)
                                                                aaTCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACg                                    <  2:3607884/100‑1 (MQ=255)
                                                                   cACCCGAACCTGCTGCCCGCTTATCTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgtgt                                 >  2:3417186/1‑100 (MQ=255)
                                                                     cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCg                               >  1:3501117/1‑100 (MQ=255)
                                                                     cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCGTCTCCCGTGAAGGCGCAAACGTGTCg                               <  2:2679323/100‑1 (MQ=255)
                                                                       cGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGAt                             >  2:1211963/1‑100 (MQ=255)
                                                                           ccttttgcCTGCTTATTTTGCGTAAGTTGATTATCCGCAATTGGAAAAATATGTGGTTAACCCTATCTTCTCCCTTGATGGCGTAAACgtg                                  <  1:2862325/86‑1 (MQ=255)
                                                                            ttgctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTg                        <  1:92532/99‑1 (MQ=255)
                                                                            ctgctgCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTg                        <  2:1137921/100‑1 (MQ=255)
                                                                            ctgctgCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTg                        <  2:3741451/100‑1 (MQ=255)
                                                                             tgctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTga                       <  1:842124/100‑1 (MQ=255)
                                                                             tgctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTga                       >  1:1675311/1‑100 (MQ=255)
                                                                             tgctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTg                        <  2:3927033/99‑1 (MQ=255)
                                                                              gctgcCCGCTTATTTTGCGGAAGATGCTCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCGCCCGTGAAGGCGCAAACGTGTCGATCATTgag                      <  1:2418220/100‑1 (MQ=255)
                                                                                tgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGaa                    >  1:783814/1‑100 (MQ=255)
                                                                                tgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGaa                    >  1:1763333/1‑100 (MQ=255)
                                                                                    cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                >  2:378898/1‑100 (MQ=255)
                                                                                    cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                <  2:2883613/100‑1 (MQ=255)
                                                                                    cGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGc                >  2:1803708/1‑100 (MQ=255)
                                                                                    cGCTTATTTTACGGAAGATGATCATCCCCAAATGGCAAAATATGTGGTTAAACCGATCTGCTCCCGTGAAGGCACAATCGTGTCGATCATTGAGAAc                   <  1:3005833/97‑1 (MQ=255)
                                                                                        ccatttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCaaaa            >  2:2183584/4‑100 (MQ=38)
                                                                                         aTTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc           >  2:1690996/1‑100 (MQ=255)
                                                                                          atttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAcc          >  2:2293952/2‑100 (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      <  1:3460782/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:4232388/1‑100 (MQ=255)
                                                                                                 gAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAag   >  2:742993/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
ACCGGCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGC  >  NC_000913/3181191‑3181388

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