Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,048,588 C→A 7.2% P77P (CCG→CCT yccM ← putative 4Fe‑4S membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,048,5880CA 7.2% 1375.9 / 50.8 355P77P (CCG→CCTyccMputative 4Fe‑4S membrane protein
Reads supporting (aligned to +/- strand):  ref base C (170/159);  new base A (10/16);  total (180/175)
Fisher's exact test for biased strand distribution p-value = 2.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 359 total aligned reads displayed.

GTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTAATAATTGCGTCGCTTTACGCCAGGTCGTG  >  CP009273/1048439‑1048735
                                                                                                                                                     |                                                                                                                                                   
gTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATc                                                                                                                                                     >  2:2988265/1‑150 (MQ=255)
   ccAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGt                                                                                                                                                  >  1:3926826/1‑150 (MQ=255)
    cAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTa                                                                                                                                                 <  1:1528039/150‑1 (MQ=255)
     aaCCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTag                                                                                                                                                <  2:3744168/150‑1 (MQ=255)
        cGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTagcca                                                                                                                                             <  1:185480/150‑1 (MQ=255)
           gCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTagccagcc                                                                                                                                          >  2:3338162/1‑150 (MQ=255)
             aTAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATAGGTAGCCAGCCtt                                                                                                                                        >  2:3833075/1‑150 (MQ=255)
              tAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTc                                                                                                                                       <  2:2038227/150‑1 (MQ=255)
              tAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTc                                                                                                                                       <  1:1737805/150‑1 (MQ=255)
               aaCCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCg                                                                                                                                      >  2:1959892/1‑150 (MQ=255)
                aCCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGa                                                                                                                                     >  1:69897/1‑150 (MQ=255)
                  cATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATa                                                                                                                                   >  1:1003842/1‑150 (MQ=255)
                  cATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATa                                                                                                                                   >  1:549496/1‑150 (MQ=255)
                   aTGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATCTTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGAtcc                                                                                                                                  >  1:3821831/1‑148 (MQ=255)
                     gAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCt                                                                                                                                >  2:4431006/1‑150 (MQ=255)
                      aGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTc                                                                                                                               >  1:2539332/1‑150 (MQ=255)
                        cAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTccc                                                                                                                             >  1:2234975/1‑150 (MQ=255)
                         aaaaGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCg                                                                                                                            >  2:59142/1‑150 (MQ=255)
                          aaaGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCgg                                                                                                                           <  1:4548548/150‑1 (MQ=255)
                             ggCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGTCGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCt                                                                                                                        <  2:31466/150‑1 (MQ=255)
                             ggCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCt                                                                                                                        >  1:2077537/1‑150 (MQ=255)
                                 tttttGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCg                                                                                                                    <  1:4785354/150‑1 (MQ=255)
                                   tttGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGAc                                                                                                                  <  2:2107442/150‑1 (MQ=255)
                                   gtggAGTAAGAGCCGGATGACGATACCAGCGACCAATCCCCGCCTCGCGGCGGCGTGGACGGACGGGACCTGACCTGGTACAAGGCTATATTTCAGAGTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGAc                                                                                                                  <  2:1673697/147‑1 (MQ=255)
                                    ttGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACa                                                                                                                 >  1:2940772/1‑150 (MQ=255)
                                       aGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAAt                                                                                                              >  1:2887160/1‑150 (MQ=255)
                                       aGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATAGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAAt                                                                                                              <  2:1146296/150‑1 (MQ=255)
                                        gTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATg                                                                                                             <  1:4920908/150‑1 (MQ=255)
                                         tAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGt                                                                                                            <  2:983451/150‑1 (MQ=255)
                                           aTAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCg                                                                                                          <  2:3094560/150‑1 (MQ=255)
                                                 tnatgaCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACACATGTCGTGCTAc                                                                                                    >  1:4441300/1‑150 (MQ=255)
                                                      aCGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGcc                                                                                               <  1:3326408/150‑1 (MQ=255)
                                                         aTAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGtt                                                                                            >  1:276839/1‑150 (MQ=255)
                                                            aaaGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCg                                                                                         <  1:1197703/150‑1 (MQ=255)
                                                            aaaGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCg                                                                                         >  1:975446/1‑150 (MQ=255)
                                                            aaaGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCg                                                                                         >  1:3114822/1‑150 (MQ=255)
                                                             aaGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGt                                                                                        <  2:163535/150‑1 (MQ=255)
                                                              aGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTa                                                                                       >  1:4107664/1‑150 (MQ=255)
                                                              aGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTa                                                                                       <  1:2385319/150‑1 (MQ=255)
                                                                cGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATAGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAAt                                                                                     <  1:1670737/150‑1 (MQ=255)
                                                                    cAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAgcg                                                                                 <  1:1302587/150‑1 (MQ=255)
                                                                         aCAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACATATGTCGTGCTACTTGCCGTTTCGTAATAGCGCAccc                                                                            <  2:1259752/150‑1 (MQ=255)
                                                                         aCAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCAccc                                                                            >  1:2465881/1‑150 (MQ=255)
                                                                          cAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCa                                                                           >  1:3712309/1‑150 (MQ=255)
                                                                             cATCGCGGCGGCGTGGAAGGACGGTAACTGACCTGTTACAATGCTATATTGCAGATTCATCAGCCCGGCAATCGGTAGCCCGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGGTTCGTAATAGCGCACCCAATa                                                                        >  2:2059249/1‑150 (MQ=255)
                                                                              aTCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATAc                                                                       >  1:2304685/1‑150 (MQ=255)
                                                                                cgcgGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCa                                                                     <  2:3594004/150‑1 (MQ=255)
                                                                                cgcgGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCa                                                                     <  2:1557092/150‑1 (MQ=255)
                                                                                  cggcggcgTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGa                                                                   <  1:401045/150‑1 (MQ=255)
                                                                                   ggcggcgTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGaa                                                                  >  1:4532215/1‑150 (MQ=255)
                                                                                   ggcggcgTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGaa                                                                  <  1:508595/150‑1 (MQ=255)
                                                                                   ggcggcgTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGaa                                                                  >  1:3389055/1‑150 (MQ=255)
                                                                                      ggcgTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGt                                                               >  1:4638240/1‑150 (MQ=255)
                                                                                       gcgTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTg                                                              <  2:1899404/150‑1 (MQ=255)
                                                                                         gTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGAt                                                            >  1:3209774/1‑150 (MQ=255)
                                                                                         gTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGAt                                                            >  1:4285436/1‑150 (MQ=255)
                                                                                            gacggacgGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGc                                                         >  2:510798/1‑150 (MQ=255)
                                                                                            gacggacgGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATAGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGc                                                         <  2:176484/150‑1 (MQ=255)
                                                                                                gacggTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAata                                                     <  1:30004/150‑1 (MQ=255)
                                                                                                  cggTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAAtata                                                   >  1:1140983/1‑150 (MQ=255)
                                                                                                   ggTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATaa                                                  <  2:4668650/150‑1 (MQ=255)
                                                                                                   ggTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATaa                                                  <  2:3658441/150‑1 (MQ=255)
                                                                                                     tAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAat                                                >  1:3622128/1‑150 (MQ=255)
                                                                                                        cTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCTACAAATTTCGTGATACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATAtt                                             <  1:844117/150‑1 (MQ=255)
                                                                                                        cTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATAtt                                             >  2:4986902/1‑150 (MQ=255)
                                                                                                        cTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATAtt                                             <  1:1577867/150‑1 (MQ=255)
                                                                                                         tGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATAGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTc                                            <  1:1065950/150‑1 (MQ=255)
                                                                                                           aCCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCAt                                          <  1:917838/150‑1 (MQ=255)
                                                                                                                ttACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCa                                     <  2:456807/150‑1 (MQ=255)
                                                                                                                  aCAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCACTACCAGAACGTGATGGCAATATAAATATTCATTGCCAAc                                   <  2:1371025/150‑1 (MQ=255)
                                                                                                                   cAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACa                                  <  2:4082619/150‑1 (MQ=255)
                                                                                                                   cAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACa                                  >  1:3767633/1‑150 (MQ=255)
                                                                                                                   cAAGGCTAGATTTCAGATTCATCAGCCCGGCAATAGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACa                                  <  2:28050/150‑1 (MQ=255)
                                                                                                                    aaGGCTATATTTCAGATTCATCAGCCCGGCAATAGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAg                                 <  1:199861/150‑1 (MQ=255)
                                                                                                                     aGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGt                                <  2:1317787/150‑1 (MQ=255)
                                                                                                                     aGGCTATATTTCAGATTCATCAGCCCGGCAATAGGGAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGt                                <  2:4896195/150‑1 (MQ=255)
                                                                                                                      ggCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGta                               <  1:5026414/150‑1 (MQ=255)
                                                                                                                       gCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGtaa                              <  1:5035217/150‑1 (MQ=255)
                                                                                                                        cTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGtaat                             <  2:1269828/150‑1 (MQ=255)
                                                                                                                        cTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGtaat                             >  1:664376/1‑150 (MQ=255)
                                                                                                                        cTATATTTCAGAGTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGtaat                             <  2:3228059/150‑1 (MQ=255)
                                                                                                                          atatTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGtaataa                           >  2:2515130/1‑150 (MQ=255)
                                                                                                                           tatTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTaataat                          <  1:1388384/150‑1 (MQ=255)
                                                                                                                                cAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTAATAATTGCGt                     >  1:3497475/1‑150 (MQ=255)
                                                                                                                                   aTTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTAATAATTGCGTCGc                  >  2:2735163/1‑150 (MQ=255)
                                                                                                                                    ttCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTAATAATTGCGTCGCt                 <  2:4247722/150‑1 (MQ=255)
                                                                                                                                     tcatcaGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTAATAATTGCGTCGCtt                >  1:3440018/1‑150 (MQ=255)
                                                                                                                                      catcaGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTAATAATTGCGTCGCttt               <  2:3995976/150‑1 (MQ=255)
                                                                                                                                       atcaGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTAATATTTGCGTCGCTTTa              <  2:4765517/150‑1 (MQ=255)
                                                                                                                                       atcaGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTAATAATTGCGTCGCTTTa              >  1:3167413/1‑150 (MQ=255)
                                                                                                                                       atcaGCCCGGCAATAGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTAATAATTGCGTCGCTTTa              <  2:925838/150‑1 (MQ=255)
                                                                                                                                        tcaGCCCGGCAATCGGTAGCCAGACTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGATCGTGATGGCAATATAAATATTCATTGCCAACAGTAATAATTGCGTCGCTTTAc             <  1:3066713/150‑1 (MQ=255)
                                                                                                                                           gCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTAATAATTGCGTCGCTTTACGcc          >  2:702416/1‑150 (MQ=255)
                                                                                                                                            cccGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTAATAATTGCGTCGCTTTACGCCa         <  2:2180993/150‑1 (MQ=255)
                                                                                                                                             ccGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACTAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTAATAATTGCGTCGCTTTACGCCAg        <  1:1787236/150‑1 (MQ=255)
                                                                                                                                             ccGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTAATAATTGCGTCGCTTTACGCCAg        >  1:4173947/1‑150 (MQ=255)
                                                                                                                                              cGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTAATAATTGCGTCGCTTTACGCCAgg       >  2:1069649/1‑150 (MQ=255)
                                                                                                                                              cGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTAATAATTGCGTCGCTTTACGCCAgg       >  2:2549797/1‑150 (MQ=255)
                                                                                                                                              cGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTAATAATTGCGTCGCTTTACGCCAgg       >  1:746892/1‑150 (MQ=255)
                                                                                                                                                   aTCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTAATAATTGCGTCGCTTTACGCCAGGTCGTg  >  2:2332266/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTACCAACCGGGCATAACCATGAGCAAAAGGCCTTTTTGAGTAATAGACTGATGACGATAAAAGCGACCAATAACAGCATCGCGGCGGCGTGGACGGACGGTAACTGACCTGTTACAAGGCTATATTTCAGATTCATCAGCCCGGCAATCGGTAGCCAGCCTTCGATACCTCCCGGTCTGGCGACAAATGTCGTGCTACTTGCCGTTTCGTAATAGCGCACCCAATACCAGAACGTGATGGCAATATAAATATTCATTGCCAACAGTAATAATTGCGTCGCTTTACGCCAGGTCGTG  >  CP009273/1048439‑1048735

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