Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,679,181 T→G R3R (AGG→CGG)  yhjH ← cyclic‑di‑GMP phosphodiesterase, FlhDC‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,679,1810TG98.6% 742.9 / ‑2.6 210R3R (AGG→CGG) yhjHcyclic‑di‑GMP phosphodiesterase, FlhDC‑regulated
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base G (100/107);  total (101/109)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 210 total aligned reads displayed.

TGCCAGGTGTAAGCACGCTCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCTTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTACCGCCTGGCCGGGGT  >  NC_000913/3679084‑3679277
                                                                                                 |                                                                                                
tGCCAGGTGTAAGCACGCTCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTa                                                                                                <  2:3275999/100‑1 (MQ=255)
 gCCAGGTGTAAGCACGCTCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTat                                                                                               >  1:4492672/1‑100 (MQ=255)
 gCCAGGTGTAAGCACGCTCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTat                                                                                               <  2:2281688/100‑1 (MQ=255)
   cAGGTGTAAGCACGCTCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTatca                                                                                             <  2:3240372/100‑1 (MQ=255)
   cAGGTGTAAGCACGCTCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTatca                                                                                             >  2:2715118/1‑100 (MQ=255)
   cAGGTGTAAGCACGCTCAAAAAGAAACAAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTatca                                                                                             <  2:5017446/100‑1 (MQ=255)
    aGGTGTAAGCACGCTCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTatcat                                                                                            <  2:3318038/100‑1 (MQ=255)
    aGGTGTAAGCACGCTCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTatcat                                                                                            >  1:4199114/1‑100 (MQ=255)
    aGGTGTAAGCACGCTCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTatcat                                                                                            <  1:3155710/100‑1 (MQ=255)
      gtgtAAGCACGCTCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTAtcatca                                                                                          <  1:105074/100‑1 (MQ=255)
        gtAAGCACGCTCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCaaa                                                                                        >  2:746588/1‑100 (MQ=255)
        gtAAGCACGCTCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCaaa                                                                                        >  1:268004/1‑100 (MQ=255)
        gtAAGCACGCTCACACTGCAACCAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTTGCTGGATAACCTGCCGTATCATCAAAg                                                                                       >  1:4659718/1‑100 (MQ=255)
         tAAGCACGCTCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAg                                                                                       <  2:260905/100‑1 (MQ=255)
         tAAGCACGCTCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATACTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAg                                                                                       <  1:3337239/100‑1 (MQ=255)
         tAAGCACGCTCACACTGCAACCAAAAACGCCGCTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGGATCATCAAAg                                                                                       >  1:5000828/1‑100 (MQ=255)
           aGCACGCTCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGAt                                                                                     <  1:2930950/100‑1 (MQ=255)
            gCACGCTCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATg                                                                                    >  2:2946716/1‑100 (MQ=255)
                  tCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCtgt                                                                              >  2:38895/1‑100 (MQ=255)
                  tCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCtgt                                                                              >  1:1962177/1‑100 (MQ=255)
                   cacaCTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCtgtg                                                                             <  2:4203557/100‑1 (MQ=255)
                     cacTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGtt                                                                           <  1:4908595/100‑1 (MQ=255)
                     cacTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGtt                                                                           <  1:4693810/100‑1 (MQ=255)
                     cacTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGtt                                                                           >  2:1621257/1‑100 (MQ=255)
                      acTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGttt                                                                          <  1:509622/100‑1 (MQ=255)
                      acTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGttt                                                                          <  2:3969389/100‑1 (MQ=255)
                      acTGCAACCAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTa                                                                         <  1:2074970/100‑1 (MQ=255)
                       cTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTa                                                                         <  2:3656805/100‑1 (MQ=255)
                       cTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTa                                                                         <  1:2558846/100‑1 (MQ=255)
                       cTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTa                                                                         <  2:771237/100‑1 (MQ=255)
                         gtaaCCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCTTATCATCAAAGATGTCCTGTGTTTaaa                                                                       <  2:5467291/98‑1 (MQ=255)
                         gCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTaaa                                                                       >  2:4478336/1‑100 (MQ=255)
                         gCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTaaa                                                                       <  1:1606630/100‑1 (MQ=255)
                          cAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTaaaa                                                                      >  2:3635094/1‑100 (MQ=255)
                           aaCCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAg                                                                     <  2:3725019/100‑1 (MQ=255)
                           aaCCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAg                                                                     <  2:1072696/100‑1 (MQ=255)
                           aaCCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAg                                                                     >  1:2244879/1‑100 (MQ=255)
                           aaCCAAAAACGCCGTTCCTGCAAGCCCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAg                                                                     <  1:1650597/100‑1 (MQ=255)
                             ccAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTg                                                                   <  1:2800028/100‑1 (MQ=255)
                              cAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGAGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGa                                                                  <  2:3985714/100‑1 (MQ=255)
                               aaaaaCGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGAt                                                                 <  1:1415526/100‑1 (MQ=255)
                               aaaaaCGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGAt                                                                 >  2:1173620/1‑100 (MQ=255)
                               aaaaaCGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGAt                                                                 >  1:2592534/1‑100 (MQ=255)
                               aaaaaCGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGAt                                                                 >  2:345311/1‑100 (MQ=255)
                               aaaaaCGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGAt                                                                 <  2:2567124/100‑1 (MQ=255)
                                     gCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCTTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGcc                                                           <  1:3037678/100‑1 (MQ=255)
                                      ccGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGccc                                                          <  2:2118240/100‑1 (MQ=255)
                                        gTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCgg                                                        >  1:5015983/1‑100 (MQ=255)
                                        gTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCgg                                                        >  1:1829200/1‑100 (MQ=255)
                                           ccTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACt                                                     >  2:313052/1‑100 (MQ=255)
                                            cTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTCATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTGTAAAAGTGATGCTGCCCGGACTc                                                    >  2:457730/1‑100 (MQ=255)
                                            cTGCAAGCTCTCGATGCCTGCTTCCGGGTGGCTTATTAACTGGATAACCTGCCGTATCAGCAAAGATGTTCTGTGGGTAAAAAAGATGCTACCCgcggtc                                                    >  2:4761195/1‑95 (MQ=255)
                                             tGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCg                                                   <  1:2373266/100‑1 (MQ=255)
                                                aaGCGCTCGATGCCTGATGCAGGGATGCTTATTCGCTGGATGACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGGCTCGACa                                                <  2:822669/100‑1 (MQ=255)
                                                  gCTCTCGATGCTTGCTCCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACgaa                                              <  1:4083459/100‑1 (MQ=255)
                                                   ctctCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAg                                             >  1:3704355/1‑100 (MQ=255)
                                                    tctcGATGCTTGCTTCAGGGTTGCATATTAGCTGGATAACCTGCCGTATCATCAGAGATGACCTGTGTTTAAACGTGAGGCTGCCCGGACTCGACAAAGa                                            >  1:5198352/1‑100 (MQ=255)
                                                     ctcGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGAt                                           >  2:1805672/1‑100 (MQ=255)
                                                       cGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTa                                         <  1:120562/100‑1 (MQ=255)
                                                       cGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTa                                         >  1:33892/1‑100 (MQ=255)
                                                         aTGCTTGCTCCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGCGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATc                                       >  2:1526006/1‑100 (MQ=255)
                                                          tgcttgctTCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCg                                      >  2:2737545/1‑100 (MQ=255)
                                                             ttgcttCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGgcg                                   >  2:5213548/1‑100 (MQ=255)
                                                             ttgcttCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGgcg                                   <  1:5519579/100‑1 (MQ=255)
                                                                 ttCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTa                               >  1:4009759/1‑100 (MQ=255)
                                                                  tCAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTaa                              <  2:4621840/100‑1 (MQ=255)
                                                                   cAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAg                             <  1:2334855/100‑1 (MQ=255)
                                                                   cAGGGTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAg                             <  2:2406062/100‑1 (MQ=255)
                                                                     gggTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAggg                           >  1:312506/1‑100 (MQ=255)
                                                                      ggTTGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGc                          <  1:4634582/100‑1 (MQ=255)
                                                                        ttGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAg                        <  2:2023466/100‑1 (MQ=255)
                                                                         tGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGa                       >  2:120007/1‑100 (MQ=255)
                                                                         tGCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGa                       <  1:2283782/100‑1 (MQ=255)
                                                                          gCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGaa                      >  2:1829804/1‑100 (MQ=255)
                                                                          gCTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGaa                      >  2:4792307/1‑100 (MQ=255)
                                                                           cTTATTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAAc                     <  1:5459298/100‑1 (MQ=255)
                                                                              aTTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACttt                  <  2:581694/100‑1 (MQ=255)
                                                                              aTTCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACttt                  >  1:3288721/1‑100 (MQ=255)
                                                                                tCGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTAc                <  1:942859/100‑1 (MQ=255)
                                                                                 cGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTAcc               >  2:621590/1‑100 (MQ=255)
                                                                                 cGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTAcc               >  2:1390767/1‑100 (MQ=255)
                                                                                 cGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTAcc               <  1:555716/100‑1 (MQ=255)
                                                                                 cGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTAcc               >  1:4237221/1‑100 (MQ=255)
                                                                                 cGCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGACAGAACTTTAcc               >  1:658544/1‑100 (MQ=255)
                                                                                  gCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTACCg              >  1:1740953/1‑100 (MQ=255)
                                                                                  gCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTACCg              <  1:1167876/100‑1 (MQ=255)
                                                                                  gCTGGATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTACCg              >  2:509071/1‑100 (MQ=255)
                                                                                     ggATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTACCGCCt           <  2:4457010/100‑1 (MQ=255)
                                                                                      gATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTACCGCCTg          <  2:723677/100‑1 (MQ=255)
                                                                                      gATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTACCGCCTg          >  1:2711250/1‑100 (MQ=255)
                                                                                      gATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTACCGCCTg          >  1:3872033/1‑100 (MQ=255)
                                                                                      gATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTACCGCCTg          <  2:2166192/100‑1 (MQ=255)
                                                                                      gATAACCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTACCGCCTg          <  2:359797/100‑1 (MQ=255)
                                                                                          aCCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTACCGCCTGGCCg      >  1:5339/1‑100 (MQ=255)
                                                                                          aCCTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTACCGCCTGGCCg      >  2:1453361/1‑100 (MQ=255)
                                                                                           ccTGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTACCGCCTGGCCgg     <  1:3616439/100‑1 (MQ=255)
                                                                                             tGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTACCGCCTGGCCgggg   <  1:1085254/100‑1 (MQ=255)
                                                                                             tGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTACCGCCTGGCCgggg   <  1:3247772/100‑1 (MQ=255)
                                                                                             tGCCGTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTACCGCCTGGCCgggg   <  1:3370014/100‑1 (MQ=255)
                                                                                             tGCCGTATCATCAAAGATGGCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTACCGCCTGGCCgggg   <  2:2667206/100‑1 (MQ=255)
                                                                                                 |                                                                                                
TGCCAGGTGTAAGCACGCTCACACTGCAACCAAAAACGCCGTTCCTGCAAGCTCTCGATGCTTGCTTCAGGGTTGCTTATTCGCTGGATAACCTGCCTTATCATCAAAGATGTCCTGTGTTTAAAAGTGATGCTGCCCGGACTCGACAAAGATTATCGGCGCGTAAGGGCAGAACTTTACCGCCTGGCCGGGGT  >  NC_000913/3679084‑3679277

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