Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,181,294 G→C 15.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,2940GC15.8% 365.7 / 55.4 132D292H (GAT→CAT) ygiCATP‑Grasp family ATPase
Reads supporting (aligned to +/- strand):  ref base G (60/51);  new base C (8/13);  total (68/64)
Fisher's exact test for biased strand distribution p-value = 2.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

Only 100 of 140 total aligned reads displayed.

GCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGG  >  NC_000913/3181195‑3181393
                                                                                                   |                                                                                                   
  gTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatgat                                                                                                   >  1:359440/1‑100 (MQ=255)
  gTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatgat                                                                                                   >  1:2712850/1‑100 (MQ=255)
  gTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatgat                                                                                                   <  2:888214/100‑1 (MQ=255)
  gTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatgat                                                                                                   <  2:2688821/100‑1 (MQ=255)
  gTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCAAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAgatgat                                                                                                   <  2:944297/100‑1 (MQ=255)
    ggAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATTTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatca                                                                                                 <  1:2167703/100‑1 (MQ=255)
    ggAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatca                                                                                                 <  2:1276456/100‑1 (MQ=255)
     gAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatcat                                                                                                >  2:1418248/1‑100 (MQ=255)
      nagagCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGatcatc                                                                                               >  1:2585706/2‑100 (MQ=255)
       agagCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATcc                                                                                              >  1:2082559/1‑100 (MQ=255)
       agagCATTATCTCCAACAAGGCACTTCTACCGCTACTATGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATcc                                                                                              >  2:1588690/1‑100 (MQ=255)
         agCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGc                                                                                            <  2:2700715/100‑1 (MQ=255)
         agCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGc                                                                                            <  2:1954664/100‑1 (MQ=255)
          gCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCa                                                                                           >  1:2560539/1‑100 (MQ=255)
          gCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCa                                                                                           >  1:1708437/1‑100 (MQ=255)
          gCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCa                                                                                           >  2:2684496/1‑100 (MQ=255)
           gatcttctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaa                                                                                          <  2:1664811/95‑1 (MQ=255)
           cATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCaa                                                                                          <  2:1423332/100‑1 (MQ=255)
            aTTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaa                                                                                         <  2:2393950/100‑1 (MQ=255)
            aTTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCaaa                                                                                         >  2:318167/1‑100 (MQ=255)
             ttATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaa                                                                                        <  2:571791/100‑2 (MQ=39)
             ttATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaa                                                                                        <  2:7104/100‑2 (MQ=39)
              tctctcCAACAAGGCACTTCTACCGCTACGGGGGGAGATGTTCCCGAATCGCCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaat                                                                                       <  1:1759573/98‑3 (MQ=37)
              tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaat                                                                                       <  1:1707746/100‑3 (MQ=39)
              tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATg                                                                                       <  1:470539/100‑1 (MQ=255)
              tATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGATTATTTTGCGGATGATCTTCTTCCGCAAAtt                                                                                       <  2:1961553/100‑2 (MQ=255)
               aTCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATgg                                                                                      >  1:886005/1‑100 (MQ=255)
               aTCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATgg                                                                                      <  1:2157105/100‑1 (MQ=255)
               aTCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATgg                                                                                      >  2:1628116/1‑100 (MQ=255)
                tctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaagatc                                                                                     >  1:1664806/1‑96 (MQ=255)
                tctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGa                                                                                     >  2:2378343/1‑100 (MQ=255)
                 ctcCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAACTGGaa                                                                                    >  2:1126603/1‑100 (MQ=255)
                   ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagcg                                                                                  <  2:841287/100‑8 (MQ=37)
                   ccAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGATGATCATCTTCCGCaaaataagcg                                                                                  <  2:839641/100‑8 (MQ=37)
                      acaaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAata                                                                               <  2:1013846/100‑1 (MQ=255)
                        aaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtg                                                                             <  1:2902675/100‑1 (MQ=255)
                        aaGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtg                                                                             <  2:1273960/100‑1 (MQ=255)
                         aGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATAtgt                                                                            >  1:2634412/1‑100 (MQ=255)
                            cACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGt                                                                         >  1:194632/1‑100 (MQ=255)
                             aCTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGtt                                                                        >  1:419745/1‑100 (MQ=255)
                                 cTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAc                                                                    <  1:2208431/100‑1 (MQ=255)
                                 cTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAc                                                                    <  2:313547/100‑1 (MQ=255)
                                 cTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGCGGTTAAAc                                                                    <  1:1628112/100‑1 (MQ=255)
                                  tACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAAcc                                                                   <  2:968397/100‑1 (MQ=255)
                                    ccGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGa                                                                 >  1:346176/1‑100 (MQ=255)
                                      gCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtc                                                               <  1:2561514/100‑1 (MQ=255)
                                          cTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGAtcttct                                                           >  2:1353597/1‑100 (MQ=255)
                                           tgtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTtctc                                                          >  1:1357905/1‑100 (MQ=255)
                                           tgtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTtctc                                                          >  2:587505/1‑100 (MQ=255)
                                            gtgGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTcc                                                         <  2:1926117/100‑1 (MQ=255)
                                              gggAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCg                                                       >  1:1470072/1‑100 (MQ=255)
                                              gggAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCg                                                       >  1:1505527/1‑100 (MQ=255)
                                                gagaTGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTg                                                     <  2:2792551/100‑1 (MQ=255)
                                                  gaTGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGaa                                                   <  1:2730502/100‑1 (MQ=255)
                                                   aTGTTGCCGAATCACCCGAACCTGCTGCGCGATTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAg                                                  <  1:2758866/100‑1 (MQ=255)
                                                    tGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAgg                                                 <  2:2197342/100‑1 (MQ=255)
                                                     gTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgc                                                <  1:1012524/100‑1 (MQ=255)
                                                     gTTCCCGAATCACCCGAACCTGCTGACCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTGAACCGATCTTCTCCCGTGAAgg                                                 <  1:1212467/99‑1 (MQ=255)
                                                      ttCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcg                                               <  1:56983/100‑1 (MQ=255)
                                                       tCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcgc                                              >  2:250299/1‑100 (MQ=255)
                                                       tCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGACCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcgc                                              >  1:1618930/1‑100 (MQ=255)
                                                        cccGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGgcgc                                              >  1:786103/1‑99 (MQ=255)
                                                        cccGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCa                                             <  2:2826549/100‑1 (MQ=255)
                                                        cccGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCa                                             >  2:889959/1‑100 (MQ=255)
                                                         ccGAATCACCCGAACCTGCTACCCGCTTATTTTGCAGAAGATGATCATCCGCAACTGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCaa                                            >  2:1732688/1‑100 (MQ=255)
                                                            aaTCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACg                                         <  2:1083510/100‑1 (MQ=255)
                                                             aTCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgt                                        >  1:99096/1‑100 (MQ=255)
                                                              tCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGGAGGCGCAAACgtg                                       >  2:678123/1‑100 (MQ=255)
                                                               cACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACgtgt                                      >  1:234673/1‑100 (MQ=255)
                                                                aCCCGAACCTGCTGCACGCTTATTTTGCGGAAGATGATAAAAAAAAAAAAAAAAAATATGTGGTTAAACCGATCTTCTCTCGTGAAGGCGCAAACGTGTc                                     <  1:1782910/100‑1 (MQ=255)
                                                                 cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCg                                    >  1:1804834/1‑100 (MQ=255)
                                                                 cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCg                                    >  2:1032253/1‑100 (MQ=255)
                                                                 cccGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCg                                    >  1:401460/1‑100 (MQ=255)
                                                                       ccTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCAtt                              >  1:1898802/1‑100 (MQ=255)
                                                                        ctgctgCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTg                             <  2:166421/100‑1 (MQ=255)
                                                                        ctgctgCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTg                             <  1:2032409/100‑1 (MQ=255)
                                                                        ctgctgCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTg                             <  2:234675/100‑1 (MQ=255)
                                                                         tgctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTga                            <  2:903742/100‑1 (MQ=255)
                                                                          gctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTga                            >  1:1922334/1‑99 (MQ=255)
                                                                           ctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATCTGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTgaga                          >  1:886813/1‑100 (MQ=255)
                                                                           ctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTgaga                          <  1:645954/100‑1 (MQ=255)
                                                                           ctgcCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTgaga                          >  2:347078/1‑100 (MQ=255)
                                                                              cccGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACg                       <  1:1468352/100‑1 (MQ=255)
                                                                                   ttATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCaaa                  <  2:2156913/100‑1 (MQ=255)
                                                                                   ttAGTTTGCGAAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCaaa                  <  2:209007/100‑1 (MQ=255)
                                                                                    tATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCaaaa                 >  1:244574/1‑100 (MQ=255)
                                                                                    tATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCaaaa                 >  1:1877468/1‑100 (MQ=255)
                                                                                    ccatttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCaaaa                 >  1:1217255/4‑100 (MQ=38)
                                                                                     catttGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAAc                >  1:1084579/3‑100 (MQ=38)
                                                                                     aTTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTTGATCATTGAGAACGGCAAAAc                >  1:1496618/1‑100 (MQ=255)
                                                                                      atttGCGGATGATCATCCTCCGCAAATGGAAAACTATGTGGTTAAACCGATCTTCTCCCGAGAAGACGCAAACGTAACGATCATTAAGAACGGCaaaa                 >  1:1741402/2‑98 (MQ=255)
                                                                                       tttGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCa              >  2:2338728/1‑100 (MQ=255)
                                                                                       tttGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCa              <  2:1576636/100‑1 (MQ=255)
                                                                                        ttGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACATGTCGATCATTGAGAACGGCAAAACCAt             >  2:2448783/1‑100 (MQ=255)
                                                                                         tGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAtt            >  2:672915/1‑100 (MQ=255)
                                                                                         tGCGGATGATCATCTTCCGCAAATGGAAAAATATGTGGTTAAAACGATCTTCTCCCGTGAAGGCGCAAACGTGGCGATCATTGAGAACGGCAAAACCAt             >  1:1079029/1‑99 (MQ=255)
                                                                                         tGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCAtt            <  2:2487052/100‑1 (MQ=255)
                                                                                            ggAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGaa         <  2:173065/100‑1 (MQ=255)
                                                                                                  tgatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCg   >  1:267113/1‑100 (MQ=255)
                                                                                                   gatCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCgg  >  1:2851398/1‑100 (MQ=255)
                                                                                                   |                                                                                                   
GCGTGGAAGAGCATTATCTCCAACAAGGCACTTCTACCGCTACTGTGGGAGATGTTCCCGAATCACCCGAACCTGCTGCCCGCTTATTTTGCGGAAGATGATCATCCGCAAATGGAAAAATATGTGGTTAAACCGATCTTCTCCCGTGAAGGCGCAAACGTGTCGATCATTGAGAACGGCAAAACCATTGAAGCAGCGG  >  NC_000913/3181195‑3181393

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑