Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,931,836 T→G 9.2% intergenic (+26/+29) ygdG → / ← fucO Ssb‑binding protein, misidentified as ExoIX/L‑1,2‑propanediol oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,931,8360TG 9.2% 151.5 / 4.1 69intergenic (+26/+29)ygdG/fucOSsb‑binding protein, misidentified as ExoIX/L‑1,2‑propanediol oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (25/36);  new base G (3/4);  total (29/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTGCG  >  NC_000913/2931687‑2931981
                                                                                                                                                     |                                                                                                                                                 
cGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGCCGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAAcgt                                                                                                                                                   >  1:915965/1‑150 (MQ=255)
cGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATTCGCCGCAAAcgt                                                                                                                                                   >  1:143808/1‑150 (MQ=255)
cGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCCGATACGCCGCAAAcgt                                                                                                                                                   >  1:1275535/1‑150 (MQ=255)
cGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGAAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAAcgt                                                                                                                                                   >  1:768540/1‑150 (MQ=255)
 taaaaGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAAcgtc                                                                                                                                                  <  1:1081731/149‑1 (MQ=255)
 gAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAAcgtc                                                                                                                                                  <  1:377384/150‑1 (MQ=255)
  aaaaGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAAcgtcg                                                                                                                                                 <  1:499197/150‑1 (MQ=255)
   aaaGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAAcgtcgt                                                                                                                                                <  2:397790/150‑1 (MQ=255)
         gcgcAAAAAATTAGAAACCCATAAAGAGATGGCGGTTCTGTGCCGCGATATTGCCCGCTTACAAACCGAATTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGACGTTTGCGGCGTATCCgg                                                                                                                                          <  2:933253/150‑1 (MQ=255)
          cgcAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAAAGTCGTATCCGGc                                                                                                                                         >  1:33171/1‑150 (MQ=255)
          cgcAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGACGTTTGCGGCGTATCCGGc                                                                                                                                         <  1:373850/150‑1 (MQ=255)
           gcAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCa                                                                                                                                        <  2:1210702/150‑1 (MQ=255)
                aaaTTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAACCCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCAGCAAACGTCGTATCCTGCATGATc                                                                                                                                   >  1:888709/1‑150 (MQ=255)
                  aTTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATcac                                                                                                                                 <  2:705746/150‑1 (MQ=255)
                          cccATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAgcgc                                                                                                                         <  1:937721/150‑1 (MQ=255)
                           ccATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGCCTAAGCCGCAAACCTCGTATGCGGGAATATAACATCACCGCa                                                                                                                        >  2:428253/1‑150 (MQ=255)
                             aTAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCTGCAAACGTCGTATCCGGCATTATCACATCAGCGCAtt                                                                                                                      <  2:1290115/150‑1 (MQ=255)
                              tAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTTCAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCAttt                                                                                                                     >  1:898817/1‑150 (MQ=255)
                              tAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCAttt                                                                                                                     <  1:464936/150‑1 (MQ=255)
                                aaGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTAc                                                                                                                   >  2:936825/1‑150 (MQ=255)
                                 agagATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTAcc                                                                                                                  <  2:33958/150‑1 (MQ=255)
                                 agagATGGCGGTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTAcc                                                                                                                  <  2:599250/150‑1 (MQ=255)
                                   agaTGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAg                                                                                                                <  1:1254928/150‑1 (MQ=255)
                                     aTGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGc                                                                                                              <  1:276865/150‑1 (MQ=255)
                                                   ccGCGATATTGCCCGCTTCCAAACCGATTTGCATATCGACGGAAATTTACAGCAATTGCGGGTGGTACGGTAACGGCGAGCCCGGTAAGCCGCAAAAGTCGTATACGCCAATATCACAACAGCCCATTAACCAGGCGGTATTGGAAAGca                                                                                                >  2:319614/1‑149 (MQ=255)
                                                   ccGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAACCGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTTTGGTAAAGct                                                                                                >  2:23432/1‑150 (MQ=255)
                                                     gcgATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGctct                                                                                              <  1:574753/150‑1 (MQ=255)
                                                           ttGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGATCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAata                                                                                        <  2:583473/150‑1 (MQ=255)
                                                             gCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGCTACGCCGCCAACCGCGGATACGGCAGTATCACAACAGCCCATTTACCCCGCGGTTATGTAAAGCGCGAGAatagc                                                                                      >  2:450705/1‑148 (MQ=255)
                                                                cGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGCATAAGCCGCAAACGTCGTATACCGCATTATAACATCATCGCATTTACCAGGCGGTATTGTAAAACTCTAAAAAATACta                                                                                   >  2:643984/1‑149 (MQ=255)
                                                                 gCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCa                                                                                  >  1:910610/1‑150 (MQ=255)
                                                                   ttACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAg                                                                                <  1:795599/150‑1 (MQ=255)
                                                                      cAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCCTTATCACATCAGCGCATTTTCCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGt                                                                             >  1:671793/1‑150 (MQ=255)
                                                                      cAAACCGATTTCCATATCGACGGCAATTTACAGCAATTGCGGGTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGt                                                                             <  2:336434/150‑1 (MQ=255)
                                                                         aCCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATAGGCCGCAAACGGCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGc                                                                          <  2:882738/150‑1 (MQ=255)
                                                                          ccGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCt                                                                         >  1:21562/1‑150 (MQ=255)
                                                                          ccGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGAAAACGTCGTATCCGGCATTATCCCCTAAGCGCATTTACCCGACTGTATGGAAAAGCTCTACAAAATCCTCAAGAGATGca                                                                         >  1:155011/1‑149 (MQ=255)
                                                                            taTTTGCATATCGACGGCAATTTACAGCAATTTCGTTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTTCCAGGAGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGGTTc                                                                       <  1:189935/149‑1 (MQ=255)
                                                                                 gCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCCGATACGCCGCAAACGTCGTATCCGGCATTATGACATCAGCGCATTTACCAGGCGGTATGGTTAACCTCTACAATATCCTCAAGCGTTTCTTAAcgcg                                                                  >  1:228901/1‑150 (MQ=255)
                                                                                   atatCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCggg                                                                <  2:145481/150‑1 (MQ=255)
                                                                                   atatCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCggg                                                                <  1:733821/150‑1 (MQ=255)
                                                                                      tCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTg                                                             <  2:803438/150‑1 (MQ=255)
                                                                                         aCGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCGACAATATCCTCAAGCGTTGCTTCACGCGGGTTGccc                                                          >  1:693878/1‑149 (MQ=255)
                                                                                          cGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCCGATACGCCGCAAACGTCGTATACGGGATTATCCCATCAGCCCATTTACCCGGCGGTATGGTAAAGCTCTACAATATGCTCAAGCGTTTCTTAACGCGGGtttgcac                                                         >  2:1235585/1‑145 (MQ=255)
                                                                                           ggCAATTTACAGCAATTGCGGTTGGTAGGGTAACGGCGAGACGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCAcc                                                        <  2:48800/150‑1 (MQ=255)
                                                                                               aTTTCCAGCAAATGCGGTGGGGACGGAAACGGCGCGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCTGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTa                                                    <  1:1014115/150‑1 (MQ=255)
                                                                                               aTTTACATCAATTGCGGTTAGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTa                                                    <  2:63611/150‑1 (MQ=255)
                                                                                               aTTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCCGAAACGCCGCAAACCTCGTTTACGGCAATATAAAAAAAACGCATTTAACAAGCGGGATTGTAAACCAAGACAATATACACAAACCTTTGTTAACGCGGGTTGCCACCCGTa                                                    >  2:1207525/1‑150 (MQ=255)
                                                                                                  tACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACaa                                                 >  1:915631/1‑150 (MQ=255)
                                                                                                  tACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATACTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACaa                                                 >  2:30834/1‑150 (MQ=255)
                                                                                                         aTTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGACGTATACGGCATTATCACATCAGCGCATTTTCCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAATcatcat                                          >  1:846984/1‑150 (MQ=255)
                                                                                                           tGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATcc                                        <  1:580953/150‑1 (MQ=255)
                                                                                                           tGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATcc                                        <  2:470929/150‑1 (MQ=255)
                                                                                                           tGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATcc                                        <  2:466238/150‑1 (MQ=255)
                                                                                                            gCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCTGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCa                                       <  1:722132/150‑1 (MQ=255)
                                                                                                            gCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCa                                       <  1:483930/150‑1 (MQ=255)
                                                                                                             cGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAg                                      <  2:181982/150‑1 (MQ=255)
                                                                                                                ttGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTgc                                   <  1:874547/150‑1 (MQ=255)
                                                                                                                  ggTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTgccg                                 <  2:300947/150‑1 (MQ=255)
                                                                                                                    taaggTAACGGCGAGCCTGATACGCAGCAAACGTAGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTgccgcc                               <  2:860985/147‑1 (MQ=255)
                                                                                                                     aCGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCt                              >  2:751977/1‑150 (MQ=255)
                                                                                                                     aCGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCt                              >  1:23697/1‑150 (MQ=255)
                                                                                                                            cGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAg                       <  2:1210793/150‑1 (MQ=255)
                                                                                                                                gAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGcc                   <  2:910610/150‑1 (MQ=255)
                                                                                                                                 aGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCg                  <  1:533309/150‑1 (MQ=255)
                                                                                                                                  gCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCgg                 <  2:143808/150‑1 (MQ=255)
                                                                                                                                  gCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCgg                 >  1:397721/1‑150 (MQ=255)
                                                                                                                                  gCCGGATACGACGTTTGCGGCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCgg                 >  1:1114825/1‑150 (MQ=255)
                                                                                                                                   ccGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGATGCCACCGGTCCAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGa                >  1:193625/1‑150 (MQ=255)
                                                                                                                                   ccGGATACGACGTTTGCGGCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACCACCATCATCCAGTGCCGCCTGCGCCAGTGCCGGa                >  1:956174/1‑150 (MQ=255)
                                                                                                                                    cGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGaa               >  1:1057361/1‑150 (MQ=255)
                                                                                                                                      gATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTAGGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATg             <  2:158492/150‑1 (MQ=255)
                                                                                                                                         aCGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCt          <  2:985400/150‑1 (MQ=255)
                                                                                                                                              gcAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCACGCGTTGCTTCGCGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGGCTTCCtt     >  2:528051/1‑150 (MQ=255)
                                                                                                                                               cAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTg    <  1:42532/150‑1 (MQ=255)
                                                                                                                                               cAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTg    >  1:101094/1‑150 (MQ=255)
                                                                                                                                                 aaCGGCGTATCCGGCATTATCACATCAGCGCATTTCCCAGGCGGTAGGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGGTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTGCg  <  1:863319/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CGAAAAGTGGCGCAAAAAATTAGAAACCCATAAAGAGATGGCGTTTCTGTGCCGCGATATTGCCCGCTTACAAACCGATTTGCATATCGACGGCAATTTACAGCAATTGCGGTTGGTACGGTAACGGCGAGCCGGATACGCCGCAAACGTCGTATCCGGCATTATCACATCAGCGCATTTACCAGGCGGTATGGTAAAGCTCTACAATATCCTCAAGCGTTGCTTCACGCGGGTTGCCACCGGTACAAACATCATCCAGTGCCGCCTGCGCCAGTGCCGGAATGTCTTCCTTGCG  >  NC_000913/2931687‑2931981

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