Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,275,142 C→A 100% intergenic (‑1/+199) yjcH ← / ← acs DUF485 family inner membrane protein/acetyl‑CoA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,275,1420CA99.7% 1125.9 / ‑4.5 343intergenic (‑1/+199)yjcH/acsDUF485 family inner membrane protein/acetyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (190/152);  total (190/153)
Fisher's exact test for biased strand distribution p-value = 4.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01
Rejected as polymorphism: E-value score below prediction cutoff.
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 356 total aligned reads displayed.

GCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTC  >  CP009273/4274995‑4275288
                                                                                                                                                   |                                                                                                                                                  
gccGGGCGCGATGGCGATCTGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATaac                                                                                                                                                  <  2:4903842/150‑1 (MQ=255)
gccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATaac                                                                                                                                                  <  2:3199964/150‑1 (MQ=255)
 ccGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATaaca                                                                                                                                                 <  2:3099938/150‑1 (MQ=255)
  cGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTGTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAg                                                                                                                                                >  2:5168178/1‑150 (MQ=255)
     gcgcGAAGGAGATCATTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGAtt                                                                                                                                             <  1:3316567/150‑1 (MQ=255)
         gAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTcc                                                                                                                                         <  1:3710318/150‑1 (MQ=255)
          aaGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCa                                                                                                                                        <  2:2118816/150‑1 (MQ=255)
           aGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCag                                                                                                                                       >  2:177627/1‑150 (MQ=255)
            ggCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCaga                                                                                                                                      >  2:2520968/1‑150 (MQ=255)
                aTCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTa                                                                                                                                  >  1:2334512/1‑150 (MQ=255)
                aTCAGTAAAATAAAGCCGATACAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTa                                                                                                                                  >  2:2721418/1‑150 (MQ=255)
                  cAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAAt                                                                                                                                >  1:4387167/1‑150 (MQ=255)
                  cAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAAt                                                                                                                                <  1:351091/150‑1 (MQ=255)
                      aaaaTAATGCCGATATAAACTGCCAGCATAATACTCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAg                                                                                                                            <  1:710626/150‑1 (MQ=255)
                        aaTAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTAGCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAggg                                                                                                                          >  2:3802843/1‑150 (MQ=255)
                         aTAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGa                                                                                                                         >  1:235209/1‑150 (MQ=255)
                           aaaGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGAtt                                                                                                                       <  1:3437624/150‑1 (MQ=255)
                             aGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGt                                                                                                                     <  2:392830/150‑1 (MQ=255)
                               ccGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGttt                                                                                                                   <  2:2978504/150‑1 (MQ=255)
                                cGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGtttt                                                                                                                  >  2:4597438/1‑150 (MQ=255)
                                 gATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTa                                                                                                                 <  1:4378800/150‑1 (MQ=255)
                                  atatAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTCGGGATTGTTTTaa                                                                                                                >  2:5197879/1‑150 (MQ=255)
                                   tataAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAAt                                                                                                               >  1:2396758/1‑150 (MQ=255)
                                   tataAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAAt                                                                                                               <  2:3128697/150‑1 (MQ=255)
                                   tataAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGGTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAAt                                                                                                               <  2:4406772/150‑1 (MQ=255)
                                    ataAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAAtt                                                                                                              >  1:2231164/1‑150 (MQ=255)
                                    ataAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATGGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAAtt                                                                                                              <  2:3729198/150‑1 (MQ=255)
                                      aaaCTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGACATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTcc                                                                                                            >  2:2671468/1‑150 (MQ=255)
                                       aaCTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTccc                                                                                                           <  2:1615331/150‑1 (MQ=255)
                                            ccAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTTGGGATTGTTTTAATTCCCGCTcc                                                                                                      >  1:3582778/1‑150 (MQ=255)
                                            ccAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTcc                                                                                                      <  1:3902251/150‑1 (MQ=255)
                                             cAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTccc                                                                                                     >  2:1714875/1‑150 (MQ=255)
                                             cAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTccc                                                                                                     >  2:3590890/1‑150 (MQ=255)
                                               gCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCtt                                                                                                   >  1:3376099/1‑150 (MQ=255)
                                                 ataataATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTAt                                                                                                 >  2:5293230/1‑150 (MQ=255)
                                                   aataatCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATgg                                                                                               <  2:1554792/150‑1 (MQ=255)
                                                    ataatCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATggg                                                                                              <  2:4774262/150‑1 (MQ=255)
                                                       atCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGgaga                                                                                           <  1:5168933/150‑1 (MQ=255)
                                                        tCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGaa                                                                                          <  2:3997839/150‑1 (MQ=255)
                                                        tCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGaa                                                                                          <  1:4841227/150‑1 (MQ=255)
                                                         cGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAg                                                                                         >  2:2453110/1‑150 (MQ=255)
                                                         cGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAg                                                                                         >  2:2867493/1‑150 (MQ=255)
                                                         cGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAg                                                                                         >  1:2705527/1‑150 (MQ=255)
                                                             aGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTa                                                                                     >  2:3244203/1‑150 (MQ=255)
                                                             aGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTa                                                                                     >  1:90840/1‑150 (MQ=255)
                                                              ggATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAACTAGTGCCATTCATAACAGCTTCTCCGGAGGTCCTGTAGGGATTGTTTTAATTCCCGCTCCCGGGTGGGAGAAGGTTaa                                                                                    >  2:4332331/1‑150 (MQ=255)
                                                              ggATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTaa                                                                                    >  2:2932025/1‑150 (MQ=255)
                                                               gATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGTAGAAGGTTAAc                                                                                   <  1:2732546/150‑1 (MQ=255)
                                                               gATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAAc                                                                                   <  1:42528/150‑1 (MQ=255)
                                                                aTGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACg                                                                                  <  2:644277/150‑1 (MQ=255)
                                                                aTGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTGATGGGAGAAGGTTAACg                                                                                  >  2:2095185/1‑150 (MQ=255)
                                                                aTGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGGAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACg                                                                                  >  1:4680695/1‑150 (MQ=255)
                                                                    tggCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCg                                                                              <  1:1313486/150‑1 (MQ=255)
                                                                      gCAAACCGGTGCCGATTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCggg                                                                            <  1:3003950/150‑1 (MQ=255)
                                                                        aaaCCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTa                                                                          <  2:1070312/150‑1 (MQ=255)
                                                                             gTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGTGAGAAGGTTAACGCTCGGGTAACCCt                                                                     <  2:423312/150‑1 (MQ=255)
                                                                              ttGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCtt                                                                    <  1:3069130/150‑1 (MQ=255)
                                                                               tGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTg                                                                   >  2:1242985/1‑150 (MQ=255)
                                                                                gCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGc                                                                  <  1:4070714/150‑1 (MQ=255)
                                                                                  cGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCg                                                                >  1:289136/1‑150 (MQ=255)
                                                                                  cGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCg                                                                >  1:3273009/1‑150 (MQ=255)
                                                                                   gTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGa                                                               >  1:1602789/1‑150 (MQ=255)
                                                                                   gTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCACTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGa                                                               >  1:1688930/1‑150 (MQ=255)
                                                                                     ttttCGACTATCTCCCTGAATTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAAt                                                             <  2:1674870/150‑1 (MQ=255)
                                                                                         cGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAg                                                         >  2:873402/1‑150 (MQ=255)
                                                                                         cGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTGATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAg                                                         >  2:340169/1‑150 (MQ=255)
                                                                                           aCTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGc                                                       <  2:745672/150‑1 (MQ=255)
                                                                                           aCTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGc                                                       >  2:5262468/1‑150 (MQ=255)
                                                                                                    cTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAgg                                              <  2:3489592/150‑1 (MQ=255)
                                                                                                    cTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAgg                                              >  1:252407/1‑150 (MQ=255)
                                                                                                       aaaTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGt                                           >  1:737351/1‑150 (MQ=255)
                                                                                                       aaaTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGt                                           >  1:5274916/1‑150 (MQ=255)
                                                                                                        aaTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTGATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGtt                                          <  1:2051838/150‑1 (MQ=255)
                                                                                                        aaTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGtt                                          <  2:2425196/150‑1 (MQ=255)
                                                                                                        aaTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGtt                                          <  2:197780/150‑1 (MQ=255)
                                                                                                        aaTGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGtt                                          >  2:4622470/1‑150 (MQ=255)
                                                                                                          tGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTa                                        >  2:5061753/1‑150 (MQ=255)
                                                                                                             gcATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTAcgc                                     >  1:4952608/1‑150 (MQ=255)
                                                                                                             gcAGTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTAcgc                                     <  1:4630463/150‑1 (MQ=255)
                                                                                                              cATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTAcgcc                                    >  1:1160780/1‑150 (MQ=255)
                                                                                                               agtGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGACTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTAcgccg                                   <  2:3278139/148‑1 (MQ=255)
                                                                                                               aTTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTAcgccg                                   >  2:3815043/1‑150 (MQ=255)
                                                                                                               aTTGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTAcgccg                                   <  2:501600/150‑1 (MQ=255)
                                                                                                                ttGTCTTCTATCCGCTGATAAATAGTACCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTAcgccgc                                  <  2:1863556/150‑1 (MQ=255)
                                                                                                                 tGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCa                                 >  1:1708931/1‑150 (MQ=255)
                                                                                                                 tGTCTTCTATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCa                                 >  1:58312/1‑150 (MQ=255)
                                                                                                                   tcttctATCCGCTGATAAATAGTGCCTTTCATAACAGTTTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATc                               <  1:3226243/150‑1 (MQ=255)
                                                                                                                    cttctATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATcc                              <  2:4886463/150‑1 (MQ=255)
                                                                                                                     ttctATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                             <  1:1041445/150‑1 (MQ=255)
                                                                                                                     ttctATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                             <  1:2289147/150‑1 (MQ=255)
                                                                                                                     ttctATCCGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCg                             <  2:4108875/150‑1 (MQ=255)
                                                                                                                           ccGCTGATAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGcaatc                       <  2:1265313/150‑1 (MQ=255)
                                                                                                                                 aTAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGcc                 >  1:1469493/1‑150 (MQ=255)
                                                                                                                                 aTAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGcc                 <  2:2255645/150‑1 (MQ=255)
                                                                                                                                  tAAATAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCt                >  1:1754614/1‑150 (MQ=255)
                                                                                                                                   aaaTAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTg               >  1:2003203/1‑150 (MQ=255)
                                                                                                                                   aaaTAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTg               >  2:3454849/1‑150 (MQ=255)
                                                                                                                                      tAGTGCCATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATg            <  1:1830078/150‑1 (MQ=255)
                                                                                                                                            cATTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGc      >  2:5141243/1‑150 (MQ=255)
                                                                                                                                             aTTCATAACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCt     >  1:4198511/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCCGGGCGCGAAGGCGATCAGTAAAATAAAGCCGATATAAACTGCCAGCATAATAATCGACAGGATGGTGGCAAACCGTTGCCGTTTTTCGACTAACTCCCTGAAATGCGCATTGTCTTCTATCCGCTGATAAATAGTGCCATTCATCACAGATTCTCCAGAGGTAATGTAGGGATTGTTTTAATTCCCGCTCCCTTATGGGAGAAGGTTAACGCTCGGGTAACCCTTGCCGAATGTAGGCCGGATAAGGCGTTTACGCCGCATCCGGCAATCAATGCCTGATGCGACGCTGTC  >  CP009273/4274995‑4275288

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: