Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,622,058 C→A 100% 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,0580CA99.0% 983.5 / 0.1 306D281Y (GAT→TAT) rbbAribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein
Reads supporting (aligned to +/- strand):  ref base C (0/3);  new base A (149/153);  total (150/156)
Fisher's exact test for biased strand distribution p-value = 2.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-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 313 total aligned reads displayed.

CGCTGGCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCTTCGG  >  CP009273/3621911‑3622204
                                                                                                                                                   |                                                                                                                                                  
  cTGGCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATagcgc                                                                                                                                                >  2:1036173/1‑150 (MQ=255)
     gCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCt                                                                                                                                             <  1:2513213/150‑1 (MQ=255)
      cGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGAGAgcgcgcgt                                                                                                                                            >  1:2271104/1‑148 (MQ=255)
       gggCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTc                                                                                                                                           <  1:1527637/150‑1 (MQ=255)
          cagcagTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGAt                                                                                                                                        >  1:4940763/1‑150 (MQ=255)
          cagcagTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGAt                                                                                                                                        <  2:3806977/150‑1 (MQ=255)
              agTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCa                                                                                                                                    >  2:405304/1‑150 (MQ=255)
               gTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCaa                                                                                                                                   <  1:563057/150‑1 (MQ=255)
                tCCGGTGATCATTTTCATGGTTGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAAt                                                                                                                                  <  1:3670219/150‑1 (MQ=255)
                 ccGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAAtc                                                                                                                                 <  2:238784/150‑1 (MQ=255)
                   ggTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAAtctc                                                                                                                               >  2:2171822/1‑150 (MQ=255)
                    gTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATctct                                                                                                                              <  1:3294710/150‑1 (MQ=255)
                    gTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAAAGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGAAAGCGCGCTTCGATGGCAATctct                                                                                                                              >  1:2917199/1‑150 (MQ=255)
                         cATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGtt                                                                                                                         <  1:855759/150‑1 (MQ=255)
                          aTTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGttt                                                                                                                        <  2:2698700/150‑1 (MQ=255)
                            tttCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTc                                                                                                                      <  1:2136265/150‑1 (MQ=255)
                               cATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAgg                                                                                                                   >  1:4697806/1‑150 (MQ=255)
                               cATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAgg                                                                                                                   <  2:1886617/150‑1 (MQ=255)
                                   gtggtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGa                                                                                                               >  1:1383453/1‑150 (MQ=255)
                                     ggtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATa                                                                                                             <  1:1931905/150‑1 (MQ=255)
                                      gtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATAc                                                                                                            >  2:1830056/1‑150 (MQ=255)
                                      gtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATAc                                                                                                            <  2:4288469/150‑1 (MQ=255)
                                      gtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATAc                                                                                                            <  2:4169825/150‑1 (MQ=255)
                                       tggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACg                                                                                                           <  1:3529678/150‑1 (MQ=255)
                                        ggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACgg                                                                                                          >  1:286050/1‑150 (MQ=255)
                                        ggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACgg                                                                                                          <  1:4489321/150‑1 (MQ=255)
                                        ggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACgg                                                                                                          >  1:3711348/1‑150 (MQ=255)
                                          aTTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTg                                                                                                        >  2:3232094/1‑150 (MQ=255)
                                            ttACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTggg                                                                                                      <  2:521663/150‑1 (MQ=255)
                                            ttACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTggg                                                                                                      >  1:3744679/1‑150 (MQ=255)
                                             tACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGa                                                                                                     >  1:4351565/1‑150 (MQ=255)
                                                 gCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCAc                                                                                                 <  1:1635061/150‑1 (MQ=255)
                                                   aGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTACCGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTa                                                                                               <  2:4275605/150‑1 (MQ=255)
                                                    gCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTAc                                                                                              >  1:4682179/1‑150 (MQ=255)
                                                    gCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTAc                                                                                              <  2:717563/150‑1 (MQ=255)
                                                     ccGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTAcc                                                                                             <  1:405304/150‑1 (MQ=255)
                                                      cGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCg                                                                                            >  1:670789/1‑150 (MQ=255)
                                                          cGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTg                                                                                        >  1:530932/1‑150 (MQ=255)
                                                          cGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTg                                                                                        <  1:837492/150‑1 (MQ=255)
                                                           gAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGa                                                                                       >  1:1203388/1‑150 (MQ=255)
                                                             aCCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATg                                                                                     >  2:659998/1‑150 (MQ=255)
                                                              ccAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATgc                                                                                    <  2:4228694/150‑1 (MQ=255)
                                                              ccAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATgc                                                                                    >  1:5079044/1‑150 (MQ=255)
                                                              ccAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATgc                                                                                    >  2:2572388/1‑150 (MQ=255)
                                                              ccAAGAAAACCAAAAATATCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTAAGATAGCTCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATgc                                                                                    >  1:2432697/1‑150 (MQ=255)
                                                               cAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATgcg                                                                                   <  2:3976710/150‑1 (MQ=255)
                                                                aaGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATgcgc                                                                                  <  1:1574392/150‑1 (MQ=255)
                                                                aaGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATgcgc                                                                                  <  2:1499856/150‑1 (MQ=255)
                                                                 aGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGcc                                                                                 >  1:3033519/1‑150 (MQ=255)
                                                                  gAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCt                                                                                >  1:3923351/1‑150 (MQ=255)
                                                                      aCCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCg                                                                            >  1:4759436/1‑150 (MQ=255)
                                                                          aaaaTCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTgc                                                                        <  1:3855217/150‑1 (MQ=255)
                                                                           aaaTCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTgcg                                                                       <  2:1171331/150‑1 (MQ=255)
                                                                               ctcCCCGCGTGGAATGCGGAAATTAACGTGATCAACTGCAACGAAGGAACGCAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCTTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTg                                                                   <  2:3265288/150‑1 (MQ=255)
                                                                                 ccccGCGTGGAATGCGGAACTGATCGGGATCATCGGCAACGAAGGAACCATAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGTTGCGCCTGGCGTTGCGCTTGCg                                                                 <  1:5212328/150‑1 (MQ=255)
                                                                                  cccGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCgg                                                                <  2:3873214/150‑1 (MQ=255)
                                                                                    cgcgTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTa                                                              >  2:4403973/1‑150 (MQ=255)
                                                                                     gcgTGGAATGCGGAAATTAACGTGATCACCGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTaa                                                             <  2:3358196/150‑1 (MQ=255)
                                                                                      cgTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAAc                                                            <  2:3710351/150‑1 (MQ=255)
                                                                                       gTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACa                                                           >  2:3920221/1‑150 (MQ=255)
                                                                                       gTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACa                                                           >  1:4206863/1‑150 (MQ=255)
                                                                                        tGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAg                                                          >  2:4686135/1‑150 (MQ=255)
                                                                                             tGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTta                                                     >  1:24648/1‑150 (MQ=255)
                                                                                              gCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTtat                                                    <  2:2722768/150‑1 (MQ=255)
                                                                                                ggAAATTAACGTGATCAACGGCAACGAAGTAAACAAAACGCATGGAAACATAGCGCTCAACGATGGCAATCTCTGCGTGTTCAGGTTGATACGGTGGTAACACAACCGCCTGATGCGCCTGTAGATGCGCTTGCGGTAACAGATTTATaa                                                  >  2:2322535/1‑150 (MQ=255)
                                                                                                ggAAATTAACGTGATCAACGGCAACGAAGGAACCCAACCGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATaa                                                  <  2:5043890/150‑1 (MQ=255)
                                                                                                    aTTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGc                                              <  2:4908473/150‑1 (MQ=255)
                                                                                                    aTTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGc                                              >  2:582385/1‑150 (MQ=255)
                                                                                                       aaCGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGcttc                                           >  2:5143995/1‑150 (MQ=255)
                                                                                                        aCGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGcttct                                          <  2:5313781/150‑1 (MQ=255)
                                                                                                         cGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGcttctt                                         >  1:3448149/1‑150 (MQ=255)
                                                                                                           tGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTcc                                       >  2:3544254/1‑150 (MQ=255)
                                                                                                              tCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGc                                    >  1:3925003/1‑150 (MQ=255)
                                                                                                                 aCGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTa                                 >  2:2957997/1‑150 (MQ=255)
                                                                                                                    gCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAgcg                              <  2:2048069/150‑1 (MQ=255)
                                                                                                                      aaCGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCt                            >  1:4606343/1‑150 (MQ=255)
                                                                                                                      aaCGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCt                            >  1:3050885/1‑150 (MQ=255)
                                                                                                                      aaCGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCt                            <  2:1638777/150‑1 (MQ=255)
                                                                                                                       aCGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCtt                           <  2:4581991/150‑1 (MQ=255)
                                                                                                                        cGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCttt                          <  1:147580/150‑1 (MQ=255)
                                                                                                                        cGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCttt                          <  1:3148168/150‑1 (MQ=255)
                                                                                                                        cGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCttt                          >  1:1733556/1‑150 (MQ=255)
                                                                                                                         gAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTg                         >  2:5067840/1‑150 (MQ=255)
                                                                                                                           aGGAACCAAAACGCATGGTCAGTTAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCg                       <  1:1036173/150‑1 (MQ=255)
                                                                                                                           aGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCg                       >  2:3094031/1‑150 (MQ=255)
                                                                                                                              aaCCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGttt                    <  1:941218/150‑1 (MQ=255)
                                                                                                                               aCCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCTCTTTGCGTTtg                   <  1:2731131/150‑1 (MQ=255)
                                                                                                                               aCCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTtg                   <  2:3056824/150‑1 (MQ=255)
                                                                                                                               aCCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTtg                   <  1:262276/150‑1 (MQ=255)
                                                                                                                                  aaaaCGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTtgctg                >  1:4571111/1‑150 (MQ=255)
                                                                                                                                    aaCGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGcc              <  1:3647716/150‑1 (MQ=255)
                                                                                                                                      cGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGt            >  2:4516583/1‑150 (MQ=255)
                                                                                                                                        cATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAg          <  2:1556347/150‑1 (MQ=255)
                                                                                                                                         aTGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGc         <  2:7949/150‑1 (MQ=255)
                                                                                                                                         aTGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGc         >  2:726425/1‑150 (MQ=255)
                                                                                                                                         aTGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGc         <  2:3180711/150‑1 (MQ=255)
                                                                                                                                           ggTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGctc       >  1:3200288/1‑150 (MQ=255)
                                                                                                                                              cAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCTTc    >  1:3087084/1‑150 (MQ=255)
                                                                                                                                                gATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCTTCgg  >  1:3102658/1‑150 (MQ=255)
                                                                                                                                                gATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCTTCgg  <  2:1451844/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CGCTGGCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCTTCGG  >  CP009273/3621911‑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: