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.1% 680.4 / ‑1.8 213D281Y (GAT→TAT) rbbAribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (102/109);  total (103/110)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-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 224 total aligned reads displayed.

TGGCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCTTC  >  CP009273/3621914‑3622202
                                                                                                                                                |                                                                                                                                                
tGGCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATagcg                                                                                                                                               <  1:2440551/148‑1 (MQ=255)
tGGCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATagcg                                                                                                                                               <  1:2125934/148‑1 (MQ=255)
tGGCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATagcg                                                                                                                                               <  2:1541044/148‑1 (MQ=255)
 ggCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATagcgc                                                                                                                                              >  1:3321979/1‑148 (MQ=255)
  gCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATagcgcg                                                                                                                                             >  1:3226560/1‑148 (MQ=255)
  gCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATagcgcg                                                                                                                                             <  2:449783/148‑1 (MQ=255)
   cGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAgcgcgc                                                                                                                                            <  1:428945/148‑1 (MQ=255)
   cGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAgcgcgc                                                                                                                                            <  1:2806647/148‑1 (MQ=255)
       cagcagTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCg                                                                                                                                        <  1:664992/148‑1 (MQ=255)
       cagcagTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCg                                                                                                                                        <  2:2953635/148‑1 (MQ=255)
        agcagTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGa                                                                                                                                       >  1:1462477/1‑148 (MQ=255)
         gcagTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGACTGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGAt                                                                                                                                      >  1:532921/1‑148 (MQ=255)
         gcagTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGAt                                                                                                                                      >  2:3761830/1‑148 (MQ=255)
            gTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGc                                                                                                                                   <  1:1569212/148‑1 (MQ=255)
            gTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGc                                                                                                                                   >  1:3359426/1‑148 (MQ=255)
             tCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCa                                                                                                                                  >  1:1499877/1‑148 (MQ=255)
               cGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAAt                                                                                                                                <  2:1530961/148‑1 (MQ=255)
               cGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAAt                                                                                                                                <  2:3237597/148‑1 (MQ=255)
                 gTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAAtct                                                                                                                              <  2:855283/148‑1 (MQ=255)
                   gAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATctct                                                                                                                            >  2:3035262/1‑148 (MQ=255)
                   gAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATctct                                                                                                                            <  1:1516783/148‑1 (MQ=255)
                   gAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATctct                                                                                                                            <  1:2158496/148‑1 (MQ=255)
                    aGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTg                                                                                                                           >  1:166469/1‑148 (MQ=255)
                      cATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAATAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCg                                                                                                                         <  2:1090195/148‑1 (MQ=255)
                      cATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCg                                                                                                                         <  1:2245677/148‑1 (MQ=255)
                      cATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCg                                                                                                                         >  1:1576122/1‑148 (MQ=255)
                       aTTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGt                                                                                                                        <  2:2493047/148‑1 (MQ=255)
                       aTTTTCATGGTGGTGGATTAACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGt                                                                                                                        <  1:2541546/148‑1 (MQ=255)
                        ttttCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGtt                                                                                                                       <  1:2652297/148‑1 (MQ=255)
                           tCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTc                                                                                                                    <  1:3283621/148‑1 (MQ=255)
                              tggtggtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAgg                                                                                                                 <  1:92122/148‑1 (MQ=255)
                                  ggtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGa                                                                                                             <  1:1763824/148‑1 (MQ=255)
                                      gATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACg                                                                                                         <  2:2077359/148‑1 (MQ=255)
                                       aTTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACgg                                                                                                        >  1:2187679/1‑148 (MQ=255)
                                        tttACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGt                                                                                                       >  2:2129618/1‑148 (MQ=255)
                                               cAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCa                                                                                                <  2:3466028/148‑1 (MQ=255)
                                               cAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCa                                                                                                <  2:3734981/148‑1 (MQ=255)
                                                 gCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACt                                                                                              >  2:772443/1‑148 (MQ=255)
                                                  ccGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTa                                                                                             <  2:115431/148‑1 (MQ=255)
                                                         aaCCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTg                                                                                      >  1:2379658/1‑148 (MQ=255)
                                                            cAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATg                                                                                   >  1:1968666/1‑148 (MQ=255)
                                                                aaaaCCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGcc                                                                               <  2:519778/148‑1 (MQ=255)
                                                                   aCCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTgg                                                                            >  2:1047998/1‑148 (MQ=255)
                                                                     cAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCg                                                                          <  2:2022328/148‑1 (MQ=255)
                                                                     cAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCg                                                                          <  2:1973316/148‑1 (MQ=255)
                                                                     cAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCg                                                                          <  2:278403/148‑1 (MQ=255)
                                                                      aaaaaTCTCCCCGCGTGGAAGGCGGAACTTAACGTGATCACCGGCAACGACGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGt                                                                         <  1:1129220/148‑1 (MQ=255)
                                                                       aaaaTCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGtt                                                                        <  1:522236/148‑1 (MQ=255)
                                                                          aTCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTgcg                                                                     >  2:615025/1‑148 (MQ=255)
                                                                           tctcCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCGTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTgcgc                                                                    <  2:3380407/148‑1 (MQ=255)
                                                                            ctcCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCt                                                                   <  2:3065307/148‑1 (MQ=255)
                                                                              ccccGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTg                                                                 >  1:891040/1‑148 (MQ=255)
                                                                              ccccGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTg                                                                 <  2:2817002/148‑1 (MQ=255)
                                                                               cccGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGc                                                                >  1:135365/1‑148 (MQ=255)
                                                                                 cgcgTGGACTGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCgg                                                              >  1:310194/1‑148 (MQ=255)
                                                                                 cgcgTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCgg                                                              >  2:466911/1‑148 (MQ=255)
                                                                                  gcgTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGt                                                             <  1:3425630/148‑1 (MQ=255)
                                                                                   cgTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTa                                                            >  1:2289000/1‑148 (MQ=255)
                                                                                     tGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAAc                                                          >  1:3026038/1‑148 (MQ=255)
                                                                                      ggAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACa                                                         <  2:243074/148‑1 (MQ=255)
                                                                                       gAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAg                                                        >  2:23351/1‑148 (MQ=255)
                                                                                       gAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAg                                                        <  2:3359426/148‑1 (MQ=255)
                                                                                        aaTGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGa                                                       >  1:880780/1‑148 (MQ=255)
                                                                                           gCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGAttt                                                    >  1:1543439/1‑148 (MQ=255)
                                                                                           gCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGAttt                                                    <  2:3795864/148‑1 (MQ=255)
                                                                                           gCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGAttt                                                    <  2:1506665/148‑1 (MQ=255)
                                                                                           gCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGAttt                                                    >  1:3168257/1‑148 (MQ=255)
                                                                                               aaaTTAACGTGATCAACGGCAACGAAGGAACCAAATCGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACATATTTATaa                                                <  1:13778/148‑1 (MQ=255)
                                                                                               aaaTTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATaa                                                >  1:95344/1‑148 (MQ=255)
                                                                                                   tAACGTGATCAACGGCAACGCAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGc                                            <  2:2517565/148‑1 (MQ=255)
                                                                                                   tAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGc                                            >  2:1051411/1‑148 (MQ=255)
                                                                                                   tAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGc                                            >  1:1670432/1‑148 (MQ=255)
                                                                                                    aaCGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGct                                           <  2:2648319/148‑1 (MQ=255)
                                                                                                     aCGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGctt                                          <  1:1020484/148‑1 (MQ=255)
                                                                                                     aCGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGctt                                          <  2:165704/148‑1 (MQ=255)
                                                                                                     aCGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGctt                                          <  2:1366759/148‑1 (MQ=255)
                                                                                                     aCGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGctt                                          <  1:2931653/148‑1 (MQ=255)
                                                                                                        tGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGcttctt                                       >  2:1887169/1‑148 (MQ=255)
                                                                                                         gATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCttcttc                                      <  1:3557603/148‑1 (MQ=255)
                                                                                                               cGGCAACGAAGGAACCAACACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGt                                <  2:3436646/148‑1 (MQ=255)
                                                                                                                 gCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAg                              >  2:921701/1‑148 (MQ=255)
                                                                                                                    aCGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAgcgc                           >  1:1081700/1‑148 (MQ=255)
                                                                                                                     cGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCt                          >  1:1822987/1‑148 (MQ=255)
                                                                                                                     cGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCt                          >  2:1406082/1‑148 (MQ=255)
                                                                                                                     cGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCt                          <  2:3794800/148‑1 (MQ=255)
                                                                                                                      gAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCtt                         <  2:1360460/148‑1 (MQ=255)
                                                                                                                        aGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTg                       >  1:1229439/1‑148 (MQ=255)
                                                                                                                            aCCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTACCAGCGTAGCGCTTTGCGtt                   <  2:1576122/148‑1 (MQ=255)
                                                                                                                               aaaaCGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTtgc                >  1:590940/1‑148 (MQ=255)
                                                                                                                                 aaCGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTtgctg              >  1:1048939/1‑148 (MQ=255)
                                                                                                                                  aCGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTtgctgc             >  2:22584/1‑148 (MQ=255)
                                                                                                                                   cGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGcc            >  2:12053/1‑148 (MQ=255)
                                                                                                                                   cGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGcc            >  2:2039914/1‑148 (MQ=255)
                                                                                                                                       tGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAg        <  1:2529154/148‑1 (MQ=255)
                                                                                                                                        ggTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGc       <  2:532921/148‑1 (MQ=255)
                                                                                                                                         gTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGct      <  1:3498239/148‑1 (MQ=255)
                                                                                                                                         gTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGct      <  1:2872578/148‑1 (MQ=255)
                                                                                                                                          tCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGctc     <  2:2449590/148‑1 (MQ=255)
                                                                                                                                           cAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGctct    >  2:2951813/1‑148 (MQ=255)
                                                                                                                                             gATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCTTc  <  2:3321962/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                
TGGCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCTTC  >  CP009273/3621914‑3622202

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: