Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,464,207 T→C 6.3% Q473R (CAG→CGG)  algX ← alginate biosynthesis protein AlgX

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,464,2070TC 6.3% 259.3 / 4.4 90Q473R (CAG→CGG) algXalginate biosynthesis protein AlgX
Reads supporting (aligned to +/- strand):  ref base T (41/42);  new base C (2/4);  total (44/46)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCG  >  NC_002947/1464060‑1464353
                                                                                                                                                   |                                                                                                                                                  
gAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgc                                                                                                                                                  >  1:148771/1‑150 (MQ=255)
  gCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCTTCACTGGCCTGCCTGCGCGGTCTgcgc                                                                                                                                                >  1:719602/1‑150 (MQ=255)
    ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAg                                                                                                                                              >  2:887294/1‑150 (MQ=255)
    ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAg                                                                                                                                              <  2:1915642/150‑1 (MQ=255)
    ttttCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAg                                                                                                                                              <  2:628526/150‑1 (MQ=255)
       tCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCa                                                                                                                                           <  1:1412733/150‑1 (MQ=255)
        cGATCCCCTGGTAGTAACCCTGAGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCag                                                                                                                                          <  2:2330268/150‑2 (MQ=255)
         gATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACg                                                                                                                                         >  2:716924/1‑150 (MQ=255)
           tCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTTGCCTGCCTGCGCGGTCGGCGCAGGCACTGc                                                                                                                                       >  2:334652/1‑150 (MQ=255)
              ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                                                    >  2:700322/1‑150 (MQ=255)
              ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                                                    <  1:1758710/150‑1 (MQ=255)
              ccTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGaa                                                                                                                                    >  2:2400061/1‑150 (MQ=255)
               cTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgc                                                                                                                                                  >  1:1388299/1‑135 (MQ=255)
               cTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTgc                                                                                                                                                  <  2:1388299/135‑1 (MQ=255)
               cTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGacg                                                                                                                                   >  2:293060/1‑148 (MQ=255)
                 ggTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAgtt                                                                                                                                 <  2:976972/150‑1 (MQ=255)
                  gtagtaACCCTGAGGCGGTACCAGCGCGGCGTTCTCGTCACCACCAAACAAGTCCAACGCGAGCAGGGCGGGGCTGGACATTCGTTTGAACAGCGTTTTGCGACGGCTCACTGTCCTGCCGGCGTTTTTTTCGCAGTCACGGCGAAgttg                                                                                                                                >  2:1105553/1‑150 (MQ=255)
                   tagtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAgttgt                                                                                                                               <  1:1784179/150‑1 (MQ=255)
                   tagtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAgttgt                                                                                                                               >  1:651346/1‑150 (MQ=255)
                     gtaACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCCCCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGttgttg                                                                                                                             <  1:1707629/150‑1 (MQ=255)
                          ccTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGttt                                                                                                                        >  2:751736/1‑150 (MQ=255)
                           cTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTg                                                                                                                       >  2:1723915/1‑150 (MQ=255)
                            tGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGc                                                                                                                      <  2:148771/150‑1 (MQ=255)
                                gcggTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCagag                                                                                                                  <  1:2400061/150‑1 (MQ=255)
                                   tacgcAGCGCGACGTGCGTGGCAGCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCtt                                                                                                               <  2:1957836/146‑1 (MQ=255)
                                   gTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCGCTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCtt                                                                                                               >  1:613710/1‑150 (MQ=255)
                                    tgccAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTg                                                                                                              <  2:1345378/148‑1 (MQ=255)
                                    tACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTg                                                                                                              >  1:687293/1‑150 (MQ=255)
                                     aCCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCTTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTgg                                                                                                             <  2:719602/150‑1 (MQ=255)
                                      ccAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGc                                                                                                            <  1:363054/150‑1 (MQ=255)
                                         gcgcGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTc                                                                                                         >  2:1966754/1‑150 (MQ=255)
                                           gcggcgTGCGCGGCAACACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGa                                                                                                       <  1:284349/150‑1 (MQ=255)
                                                 tGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCt                                                                                                 <  1:1763864/150‑1 (MQ=255)
                                                  gcgcGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTg                                                                                                <  1:264462/150‑1 (MQ=255)
                                                   cgcgGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTgg                                                                                               <  1:700322/150‑1 (MQ=255)
                                                        caccacCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGGTGGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTGCGACCTTCTGGGTGcc                                                                                          >  2:1823033/1‑150 (MQ=255)
                                                          ccaccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGc                                                                                        <  2:687293/150‑1 (MQ=255)
                                                             ccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGAcc                                                                                                     <  2:738877/134‑1 (MQ=255)
                                                             ccaAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGAcc                                                                                                     >  1:738877/1‑134 (MQ=255)
                                                                  cAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGc                                                                                >  2:1752628/1‑150 (MQ=255)
                                                                   aaGGACAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCCGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGcc                                                                               <  2:1187305/150‑1 (MQ=255)
                                                                    aGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGccc                                                                              <  2:2257114/150‑1 (MQ=255)
                                                                       ccAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTg                                                                           >  2:1569098/1‑150 (MQ=255)
                                                                        cAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGt                                                                          >  1:901234/1‑150 (MQ=255)
                                                                            gcgAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCag                                                                                                                                          <  2:1758920/82‑2 (MQ=38)
                                                                            gcgAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCAGACCGCGCAGGCag                                                                                                                                          >  1:1758920/1‑81 (MQ=38)
                                                                                gCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGaa                                                                  <  1:751736/150‑1 (MQ=255)
                                                                                 cAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAAc                                                                 <  1:528018/150‑1 (MQ=255)
                                                                                    ggCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTCGGCTTCGACCTTCTGGGTGCCGCTTTCCgg                                                                                 <  2:1229889/131‑1 (MQ=255)
                                                                                    ggCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTCGGCTTCGACCTTCTGGGTGCCGCTTTCCgg                                                                                 >  1:1229889/1‑131 (MQ=255)
                                                                                     gCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGa                                                             <  1:1723915/150‑1 (MQ=255)
                                                                                     gCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGa                                                             <  2:1448645/150‑1 (MQ=255)
                                                                                       ggggCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGc                                                           <  2:651346/150‑1 (MQ=255)
                                                                                                aaTTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgctgg                                                  <  1:1105553/150‑1 (MQ=255)
                                                                                                 aTTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGGAAGAGCGTCGTCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCgctggc                                                 >  2:664627/1‑150 (MQ=255)
                                                                                                 aTTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCTGAATGACCTCGAACGCGAGCAGGCgctggc                                                 <  1:664627/150‑1 (MQ=255)
                                                                                                     gTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGc                                                                                                                 >  1:2440316/1‑82 (MQ=255)
                                                                                                     gTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGc                                                                                                                 <  2:2440316/82‑1 (MQ=255)
                                                                                                      tttGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCGGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGcc                                            >  1:1833996/1‑150 (MQ=255)
                                                                                                      tttGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGc                                                                                >  2:124542/1‑114 (MQ=255)
                                                                                                      tttGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGc                                                                                <  1:124542/114‑1 (MQ=255)
                                                                                                            aaGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTCAGCTGCGACCTTCTGCGTTCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTcc                                      <  1:887294/150‑1 (MQ=255)
                                                                                                              gggTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTc                                    <  1:41522/150‑1 (MQ=255)
                                                                                                                   aTGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCtt                               >  2:1777639/1‑150 (MQ=255)
                                                                                                                       gTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGt                                                                                             <  2:29309/84‑1 (MQ=255)
                                                                                                                       gTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGt                                                                                             >  1:29309/1‑84 (MQ=255)
                                                                                                                        tCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCgcgc                          >  1:1802162/1‑150 (MQ=255)
                                                                                                                          gCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCAt                        >  1:1167859/1‑150 (MQ=255)
                                                                                                                          gCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCAt                        >  1:1167866/1‑150 (MQ=255)
                                                                                                                          gCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCAt                        >  2:1275986/1‑150 (MQ=255)
                                                                                                                           ccTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCACTGGCTGGCCCAGTCATCGTCTTCGAGCATc                       >  2:546645/1‑150 (MQ=255)
                                                                                                                            cTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCCCGCATCt                      >  2:186818/1‑150 (MQ=255)
                                                                                                                               aCTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGa                   >  2:1447520/1‑150 (MQ=255)
                                                                                                                               aCTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCATCGTCTTCGCGCATCTGGa                   >  2:2212562/1‑150 (MQ=255)
                                                                                                                                   gcctgcctGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGCGATAGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGTATCTGGAAGAc               >  2:1796492/1‑150 (MQ=255)
                                                                                                                                     gtgcctgcGCAGACCGCGCAGGCACGGCGAAGATGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGa             >  1:1774425/2‑150 (MQ=255)
                                                                                                                                      tgcctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGaa            <  2:988347/150‑1 (MQ=255)
                                                                                                                                       gcctgcGCGGTCTGCGCAGGCCCGTCCAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAg           <  1:526778/150‑1 (MQ=255)
                                                                                                                                        cctgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGc          >  2:822699/1‑150 (MQ=255)
                                                                                                                                          tgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCgg        <  1:334652/150‑1 (MQ=255)
                                                                                                                                          tgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAACAGGCGctggctgg                                              >  1:1697717/1‑112 (MQ=255)
                                                                                                                                          tgcGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAACAGGCGctggctgg                                              <  2:1697717/112‑1 (MQ=255)
                                                                                                                                            cgcgGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGcc      >  2:1701648/1‑150 (MQ=255)
                                                                                                                                               ggTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCGGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTc   <  2:1833996/150‑1 (MQ=255)
                                                                                                                                               ggTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTc   <  1:1275986/150‑1 (MQ=255)
                                                                                                                                                gTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTcgc                            >  2:385314/1‑124 (MQ=255)
                                                                                                                                                gTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTcgc                            <  1:385314/124‑1 (MQ=255)
                                                                                                                                                gTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAg           <  2:1028193/141‑1 (MQ=255)
                                                                                                                                                gTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAg           >  1:1028193/1‑141 (MQ=255)
                                                                                                                                                gTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCg  <  1:1966754/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GAGCTTTTCGATCCCCTGGTAGTAACCCTGCGGCGGTACCAGCGCGGCGTGCGCGGCACCACCAAACAAGGCCAACGCGAGCAGGGCGGGGCTGGAAATTCGTTTGAAAAGGGTCATGGGTCGCCTCACTGGCCTGCCTGCGCGGTCTGCGCAGGCACGGCGAAGTTGTTGCGTTTGCAGAGCTTGGCTTCGACCTTCTGGGTGCCGCTTTCCGGCCCCTGTACCTCGAACGCGAGCAGGCGCTGGCTGGCCCAGTCCTCGTCTTCGCGCATCTGGAAGACGAAGCGGCCGTCG  >  NC_002947/1464060‑1464353

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: