Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 108,039 C→A 100% intergenic (+119/+97) mutT → / ← yacG nucleoside triphosphate pyrophosphohydrolase, marked preference for dGTP/DNA gyrase inhibitor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273108,0390CA98.0% 912.0 / 7.4 298intergenic (+119/+97)mutT/yacGnucleoside triphosphate pyrophosphohydrolase, marked preference for dGTP/DNA gyrase inhibitor
Reads supporting (aligned to +/- strand):  ref base C (5/1);  new base A (159/132);  total (165/133)
Fisher's exact test for biased strand distribution p-value = 2.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 308 total aligned reads displayed.

CGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCAC  >  CP009273/107892‑108185
                                                                                                                                                   |                                                                                                                                                  
cGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGtaggtagg                                                                                                                                                  <  2:3735323/150‑1 (MQ=255)
cGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGtaggtagg                                                                                                                                                  <  2:4921959/150‑1 (MQ=255)
cGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCGTATCAGGCCTTAAGGGCGTACTCACTCATTGGTCAATGGGTTTGTGGAGGTAgg                                                                                                                                                  >  2:1204554/1‑150 (MQ=25)
  gTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGAt                                                                                                                                                <  2:533162/150‑1 (MQ=255)
   tAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATa                                                                                                                                               <  1:1535266/150‑1 (MQ=255)
      ttGCGAAGCTTACACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAgg                                                                                                                                            <  2:4967051/150‑1 (MQ=255)
      ttGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAgg                                                                                                                                            <  1:2274074/150‑1 (MQ=255)
      ttGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAgg                                                                                                                                            >  1:3455664/1‑150 (MQ=255)
         cGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGt                                                                                                                                         >  1:328618/1‑150 (MQ=255)
         cGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCAGTGATAAATTTGTTTTTGTAGGTAGGATACGGCGt                                                                                                                                         >  2:1050321/1‑150 (MQ=255)
           aaGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTc                                                                                                                                       <  1:3163920/150‑1 (MQ=255)
             gCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCAc                                                                                                                                     >  2:1518430/1‑150 (MQ=255)
             gCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCAc                                                                                                                                     <  2:2724777/150‑1 (MQ=255)
               ttAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCAcgc                                                                                                                                   >  1:4015710/1‑150 (MQ=255)
               ttAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCAcgc                                                                                                                                   >  1:4976963/1‑150 (MQ=255)
                tAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCAcgcc                                                                                                                                  <  1:4938953/150‑1 (MQ=255)
                 aaaCGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCAcgccg                                                                                                                                 >  1:4754383/1‑150 (MQ=255)
                  aaCGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCAcgccgc                                                                                                                                >  1:3309827/1‑150 (MQ=255)
                   aCGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCa                                                                                                                               <  2:726413/150‑1 (MQ=255)
                   aCGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCa                                                                                                                               <  2:4170179/150‑1 (MQ=255)
                   aCGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCa                                                                                                                               <  1:624834/150‑1 (MQ=255)
                      tCTGTAGGTCAGATACGGCGTTGTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCGTATCAGGCCTAACGGGATGACTAAGTCATTGATAAATGAGTGGGTGGAGGTATGCATAGGCGTTCACGCCGCagca                                                                                                                            >  2:787002/1‑147 (MQ=21)
                      tCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATcc                                                                                                                            <  2:1274752/150‑1 (MQ=255)
                       cTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCg                                                                                                                           <  2:2605662/150‑1 (MQ=255)
                          tAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACa                                                                                                                        <  2:2068554/150‑1 (MQ=255)
                          tAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACa                                                                                                                        <  2:297977/150‑1 (MQ=255)
                          tAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACa                                                                                                                        >  2:2475765/1‑150 (MQ=255)
                             gacAGATAAGGCGTTTTCGCAGCATCCGACATTCGCACACGATGCCTGATGCGACTCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATACGACAtat                                                                                                                     >  2:4510590/3‑148 (MQ=255)
                             gTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACTCTGGCGCGTCTTATCAGGCCTAAAGTGATTTCTAACTCATTGATAATTTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACAttt                                                                                                                     <  2:2299617/150‑1 (MQ=255)
                                  aTAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGcaca                                                                                                                >  1:1949123/1‑150 (MQ=255)
                                   tAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACaa                                                                                                               >  1:1749230/1‑150 (MQ=255)
                                     aGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTTGGTAGGATAAGGCGTTCACGCCGCATCCGACCTTTGCACAAGa                                                                                                             >  1:2967399/1‑150 (MQ=255)
                                        cGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGc                                                                                                          >  1:1707354/1‑150 (MQ=255)
                                         gTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGcc                                                                                                         >  1:3188135/1‑150 (MQ=255)
                                         gTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGcc                                                                                                         >  1:820336/1‑150 (MQ=255)
                                         gTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGcc                                                                                                         >  2:201030/1‑150 (MQ=255)
                                            ttCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGa                                                                                                      >  1:1224674/1‑150 (MQ=255)
                                            ttCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGGAGGATGAGGCGTGCACGCCGCCTCCGACCTGTGCACAAGATGCCTGa                                                                                                      >  2:2758186/1‑150 (MQ=38)
                                              cgccgcATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATg                                                                                                    <  2:3804358/150‑1 (MQ=255)
                                              cgccgcATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATg                                                                                                    <  1:3351692/150‑1 (MQ=255)
                                              cgccgcATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATg                                                                                                    <  2:1751790/150‑1 (MQ=255)
                                               gccgcATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGc                                                                                                   <  1:1518430/150‑1 (MQ=255)
                                                ccgcATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCg                                                                                                  <  2:4335431/150‑1 (MQ=255)
                                                   cATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACg                                                                                               >  2:2586982/1‑150 (MQ=255)
                                                   cATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACg                                                                                               >  1:351610/1‑150 (MQ=255)
                                                   cATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACg                                                                                               <  2:812767/150‑1 (MQ=255)
                                                     tCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATTAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCt                                                                                             <  1:4253563/150‑1 (MQ=255)
                                                       cGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGt                                                                                           >  1:3561472/1‑150 (MQ=255)
                                                        gACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTc                                                                                          >  1:2142167/1‑150 (MQ=255)
                                                         aCATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTcc                                                                                         <  1:4244442/150‑1 (MQ=255)
                                                           aTTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATACATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCcgc                                                                                       >  2:3325831/1‑150 (MQ=255)
                                                           aTTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCcgc                                                                                       >  2:4634186/1‑150 (MQ=255)
                                                           aTTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCcgc                                                                                       >  2:904432/1‑150 (MQ=255)
                                                           aTTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCcgc                                                                                       >  1:3718179/1‑150 (MQ=255)
                                                                 acacGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTa                                                                                 <  2:652468/150‑1 (MQ=255)
                                                                 acacGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGAGTTCTCACTCATTGATAAATGTGTTTGTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGTTGCCTGTGGCGACGCTGGCCGCGTcgtg                                                                                 >  2:4395135/1‑147 (MQ=38)
                                                                  cacGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTAt                                                                                <  2:3047170/150‑1 (MQ=255)
                                                                  cacGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTAt                                                                                >  2:2186719/1‑150 (MQ=255)
                                                                   acGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATc                                                                               >  1:3797627/1‑150 (MQ=255)
                                                                      aTGCCTGATGCGACGCTGGCGCGTCTTGTCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAgg                                                                            >  2:2901938/1‑150 (MQ=255)
                                                                         ccTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTTATAAATTTGTTTTTTTAGGTAGGATAAGGCGTTAAAGCCGCATCCGACATTTGCACAAGATGCCTGATGAGACGATGTCCGCGTCTTATCAGGCCt                                                                         >  1:1095310/1‑150 (MQ=255)
                                                                           tGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTTGGATAAGGCGTTAACGCCGCATCCGACATTTGCACAAGATTCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTAc                                                                       >  2:4834419/1‑150 (MQ=255)
                                                                            gATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACg                                                                      >  1:3056097/1‑150 (MQ=255)
                                                                               gCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCAACCGACATTTGAACAATATGCATGATGCGACGAAGTCCGCGTCTTATCAGGCCTACGTGc                                                                   >  2:1162330/1‑150 (MQ=38)
                                                                                cGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCg                                                                  <  1:2475765/150‑1 (MQ=255)
                                                                                cGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCg                                                                  >  2:4291722/1‑150 (MQ=255)
                                                                                 gACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCgg                                                                 >  1:4696121/1‑150 (MQ=255)
                                                                                  aCGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGAGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGc                                                                >  2:2434933/1‑150 (MQ=255)
                                                                                     cTGGCGCTTGTTTTGAGGCCTTACTGTTTTTCTGTCTGATTGATTATTTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATc                                                             <  2:2150299/150‑1 (MQ=21)
                                                                                      tGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCa                                                            <  2:4415707/150‑1 (MQ=255)
                                                                                      tGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCa                                                            <  2:1844937/150‑1 (MQ=255)
                                                                                        gcgcGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGa                                                          <  1:4160795/150‑1 (MQ=255)
                                                                                                 aTCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCa                                                 >  2:5055319/1‑150 (MQ=255)
                                                                                                  tCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCAc                                                <  1:2777747/150‑1 (MQ=255)
                                                                                                  tCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCGACGTGCGGCATCAGACAAATGTCAc                                                >  1:3023800/1‑150 (MQ=255)
                                                                                                   cAGGCCTAAAGGGATTTCAAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCTGCATCAGACAAATGTCACt                                               >  1:3904306/1‑150 (MQ=255)
                                                                                                    aGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTg                                              <  1:4414387/150‑1 (MQ=255)
                                                                                                     ggCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGc                                             <  1:573013/150‑1 (MQ=255)
                                                                                                        cTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCttt                                          >  1:2938968/1‑150 (MQ=255)
                                                                                                         tAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTg                                         <  2:1425421/150‑1 (MQ=255)
                                                                                                          aaaGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTgg                                        >  2:845968/1‑150 (MQ=255)
                                                                                                           aaGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGt                                       <  1:4812088/150‑1 (MQ=255)
                                                                                                            aGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGtt                                      <  2:487714/150‑1 (MQ=255)
                                                                                                            aGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGtt                                      >  2:3127905/1‑150 (MQ=255)
                                                                                                               gATTTCTAACTCATTGATAAATTTTTTTTTGTAGGTATGATAAGGCGTTCACGCCGCATCCGACATTTGCACTAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttctt                                   <  1:2833333/150‑1 (MQ=255)
                                                                                                                aTTTCTCACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGttcttc                                  <  2:3560366/150‑1 (MQ=255)
                                                                                                                     tAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTcc                             >  1:3769983/1‑150 (MQ=255)
                                                                                                                     tAACTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTcc                             >  2:1878868/1‑150 (MQ=255)
                                                                                                                      aaCTCATTGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCa                            >  2:1114122/1‑150 (MQ=255)
                                                                                                                            ttGATAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATc                      >  2:324347/1‑150 (MQ=255)
                                                                                                                                tAAATTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCtt                  <  1:874700/150‑1 (MQ=255)
                                                                                                                                 aaaTTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCttt                 <  1:5104113/150‑1 (MQ=255)
                                                                                                                                   aTTTGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCg               <  2:4015710/150‑1 (MQ=255)
                                                                                                                                    tttGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCgg              <  1:201030/150‑1 (MQ=255)
                                                                                                                                     ttGTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGa             <  2:1618668/150‑1 (MQ=255)
                                                                                                                                       gTTTTTGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGaaa           <  1:3676485/150‑1 (MQ=255)
                                                                                                                                          tttGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGAt        <  1:2714652/150‑1 (MQ=255)
                                                                                                                                           ttGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATc       >  1:1873221/1‑150 (MQ=255)
                                                                                                                                            tGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCg      >  1:2075248/1‑150 (MQ=255)
                                                                                                                                            tGTAGGTAGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCg      >  2:22216/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGGTAATTGCGAAGCTTAAACGTCTGTAGGTCAGATAAGGCGTTTTCGCCGCATCCGACATTCGCACACGATGCCTGATGCGACGCTGGCGCGTCTTATCAGGCCTAAAGGGATTTCTAACTCATTGATAAATTTGTTTTTGTAGGTCGGATAAGGCGTTCACGCCGCATCCGACATTTGCACAAGATGCCTGATGCGACGCTGTCCGCGTCTTATCAGGCCTACGTGCGGCATCAGACAAATGTCACTGCTTTGGTTCTTCGCTCCAGTCATCGCTTTCGGAAAGATCGCCAC  >  CP009273/107892‑108185

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