Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,622,058 C→A 85.2% D281Y (GAT→TAT)  rbbA ← ribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,622,0580CA85.2% 737.4 / 147.7 325D281Y (GAT→TAT) rbbAribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein
Reads supporting (aligned to +/- strand):  ref base C (27/21);  new base A (121/156);  total (148/177)
Fisher's exact test for biased strand distribution p-value = 1.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.69e-01

Only 100 of 330 total aligned reads displayed.

TCGCTGGCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCTTCGG  >  CP009273/3621910‑3622204
                                                                                                                                                    |                                                                                                                                                  
  gCTGGCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATagc                                                                                                                                                  <  2:3559376/149‑1 (MQ=255)
        gggCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCtt                                                                                                                                            <  1:3472603/149‑1 (MQ=255)
         ggCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTc                                                                                                                                           >  1:2053366/1‑149 (MQ=255)
          gcagcaGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCg                                                                                                                                          <  2:1910953/149‑1 (MQ=255)
             gcagTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATg                                                                                                                                       >  1:1467858/1‑149 (MQ=255)
                      tGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCCATctct                                                                                                                              >  2:4468982/1‑149 (MQ=255)
                       gAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTg                                                                                                                             >  1:3203429/1‑149 (MQ=255)
                       gAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTg                                                                                                                             <  1:2276919/149‑1 (MQ=255)
                           aTTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGtt                                                                                                                         >  2:838721/1‑149 (MQ=255)
                           aTTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGtt                                                                                                                         <  1:2785985/149‑1 (MQ=255)
                            ttttCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGttt                                                                                                                        <  1:2921652/149‑1 (MQ=255)
                               tCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCa                                                                                                                     <  2:4653941/149‑1 (MQ=255)
                                  tggtggtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGt                                                                                                                  <  1:185490/149‑1 (MQ=255)
                                   ggtggtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGtt                                                                                                                 <  1:1917888/149‑1 (MQ=255)
                                     tggtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGa                                                                                                               <  2:2022634/149‑1 (MQ=255)
                                         ggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACg                                                                                                           <  2:895115/149‑1 (MQ=255)
                                           atgtACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGt                                                                                                         <  2:2276290/146‑1 (MQ=255)
                                            tttACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTg                                                                                                        <  2:1830263/149‑1 (MQ=255)
                                             ttACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTgg                                                                                                       <  1:1666003/149‑1 (MQ=255)
                                              tACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTggg                                                                                                      >  1:381999/1‑149 (MQ=255)
                                                ccGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGAt                                                                                                    >  2:4723822/1‑149 (MQ=255)
                                                  gCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCa                                                                                                  >  2:1055832/1‑149 (MQ=255)
                                                      ccGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTAc                                                                                              <  2:2538683/149‑1 (MQ=255)
                                                      ccGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTAc                                                                                              >  1:609304/1‑149 (MQ=255)
                                                      ccGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTAc                                                                                              <  2:3903365/149‑1 (MQ=255)
                                                        gTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCg                                                                                            >  1:3954713/1‑149 (MQ=255)
                                                           cGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCt                                                                                         >  2:2159060/1‑149 (MQ=255)
                                                           cGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCt                                                                                         >  2:689163/1‑149 (MQ=255)
                                                             aaCCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGa                                                                                       <  1:2123131/149‑1 (MQ=255)
                                                             aaCCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGa                                                                                       >  1:1136953/1‑149 (MQ=255)
                                                              aCCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCACCGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGAt                                                                                      <  2:3330432/149‑1 (MQ=255)
                                                               ccAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATg                                                                                     <  2:4145792/149‑1 (MQ=255)
                                                                cAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATgc                                                                                    >  2:3644593/1‑149 (MQ=255)
                                                                 aaGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATgcg                                                                                   >  2:4396563/1‑149 (MQ=255)
                                                                   gAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGcc                                                                                 <  1:4182052/149‑1 (MQ=255)
                                                                   gAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGcc                                                                                 <  1:1003792/149‑1 (MQ=255)
                                                                   gAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGcc                                                                                 <  1:879297/149‑1 (MQ=255)
                                                                    aaaaCCATAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCt                                                                                >  1:2759033/1‑149 (MQ=255)
                                                                    aaaaCCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCt                                                                                <  1:3081483/149‑1 (MQ=255)
                                                                     aaaCCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTg                                                                               <  1:785591/149‑1 (MQ=255)
                                                                       aCCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGc                                                                             <  1:2187816/149‑1 (MQ=255)
                                                                         cAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGt                                                                           >  1:2595738/1‑149 (MQ=255)
                                                                           aaaaTCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTg                                                                         >  1:4378706/1‑149 (MQ=255)
                                                                               tctcCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCt                                                                     >  1:4018383/1‑149 (MQ=255)
                                                                               tctcCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCt                                                                     >  1:1435885/1‑149 (MQ=255)
                                                                               tctcCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCt                                                                     <  2:1306397/149‑1 (MQ=255)
                                                                               tctcCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTGCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCt                                                                     <  2:1926416/149‑1 (MQ=255)
                                                                                  ccccGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGc                                                                  >  2:599735/1‑149 (MQ=255)
                                                                                  ccccGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGc                                                                  >  1:1625649/1‑149 (MQ=255)
                                                                                   cccGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCg                                                                 >  2:1401072/1‑149 (MQ=255)
                                                                                    ccGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCgg                                                                <  2:885543/149‑1 (MQ=255)
                                                                                     cgcgTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGt                                                               <  1:3699549/149‑1 (MQ=255)
                                                                                     cgcgTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGt                                                               <  1:702913/149‑1 (MQ=255)
                                                                                     cgcgTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGt                                                               <  1:4321249/149‑1 (MQ=255)
                                                                                      gcgTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTa                                                              <  2:2217812/149‑1 (MQ=255)
                                                                                      gcgTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTa                                                              <  2:3231714/149‑1 (MQ=255)
                                                                                      gcgTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTa                                                              <  1:111300/149‑1 (MQ=255)
                                                                                      gcgTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTa                                                              >  1:988786/1‑149 (MQ=255)
                                                                                        gTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAAc                                                            >  1:2456364/1‑149 (MQ=255)
                                                                                        gTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAAc                                                            >  1:4558878/1‑149 (MQ=255)
                                                                                            aaTGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGAt                                                        >  1:3030111/1‑149 (MQ=255)
                                                                                            aaTGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGAt                                                        <  1:2961426/149‑1 (MQ=255)
                                                                                            aaTGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGAt                                                        >  2:1978217/1‑149 (MQ=255)
                                                                                                 ggAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTtata                                                   >  1:547582/1‑149 (MQ=255)
                                                                                                 ggAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTtata                                                   <  1:1593374/149‑1 (MQ=255)
                                                                                                  gAAATTAACGTGATCAACTGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATaa                                                  <  1:2708317/149‑1 (MQ=255)
                                                                                                  gAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATaa                                                  <  1:3608361/149‑1 (MQ=255)
                                                                                                   aaaTTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATaaa                                                 <  2:3217014/149‑1 (MQ=255)
                                                                                                   aaaTTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATaaa                                                 <  1:746453/149‑1 (MQ=255)
                                                                                                    aaTTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAAt                                                >  1:1668009/1‑149 (MQ=255)
                                                                                                     aTTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATg                                               <  2:3125531/149‑1 (MQ=255)
                                                                                                      ttAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGc                                              <  2:2405932/149‑1 (MQ=255)
                                                                                                         aCGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGcttc                                           >  1:4763252/1‑149 (MQ=255)
                                                                                                          cGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGcttct                                          <  2:2388821/149‑1 (MQ=255)
                                                                                                          cGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGcttct                                          >  1:1688024/1‑149 (MQ=255)
                                                                                                             gATCAACGGCAACGAAGGAACCAACACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTcc                                       <  2:754404/149‑1 (MQ=255)
                                                                                                             gATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTcc                                       >  2:4230826/1‑149 (MQ=255)
                                                                                                             gATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTcc                                       >  2:3065762/1‑149 (MQ=255)
                                                                                                                  aCGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCACTCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGAGTTATAAATGCTTCTTCCAGCGt                                  >  1:887379/1‑149 (MQ=255)
                                                                                                                   cGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTa                                 <  1:3914494/149‑1 (MQ=255)
                                                                                                                   cGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTa                                 <  1:3306853/149‑1 (MQ=255)
                                                                                                                     gCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAgc                               >  2:3677523/1‑149 (MQ=255)
                                                                                                                     gCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAgc                               <  2:2754361/149‑1 (MQ=255)
                                                                                                                       aaCGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAgcgc                             >  1:176890/1‑149 (MQ=255)
                                                                                                                         cGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCtt                           <  2:295851/149‑1 (MQ=255)
                                                                                                                            aGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGc                        <  2:4260120/149‑1 (MQ=255)
                                                                                                                              gAACCACAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGt                      <  2:4100348/149‑1 (MQ=255)
                                                                                                                              gAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGt                      >  2:1687822/1‑149 (MQ=255)
                                                                                                                                aCCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGttt                    >  2:2861267/1‑149 (MQ=255)
                                                                                                                                aCCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGttt                    >  2:1332855/1‑149 (MQ=255)
                                                                                                                                aCCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGttt                    >  2:2336620/1‑149 (MQ=255)
                                                                                                                                aCCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGttt                    <  1:2315480/149‑1 (MQ=255)
                                                                                                                                  cAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTtgc                  <  1:1849708/149‑1 (MQ=255)
                                                                                                                                    aaaCGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTtgctg                >  2:3732141/1‑149 (MQ=255)
                                                                                                                                     aaCGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTtgctgc               >  2:2622410/1‑149 (MQ=255)
                                                                                                                                     aaCGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTtgctgc               >  1:3403363/1‑149 (MQ=255)
                                                                                                                                               cAGATAGCGCGCTTCGATTGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCtt     <  1:2380511/149‑1 (MQ=255)
                                                                                                                                               cAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCtt     <  1:2718870/149‑1 (MQ=255)
                                                                                                                                                aGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCTTc    >  2:1744296/1‑149 (MQ=255)
                                                                                                                                                  aTCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCTTCgg  >  2:4623584/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCGCTGGCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCTTCGG  >  CP009273/3621910‑3622204

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: