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,0580CA96.9% 622.9 / 8.9 207D281Y (GAT→TAT) rbbAribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein
Reads supporting (aligned to +/- strand):  ref base C (5/2);  new base A (90/110);  total (95/112)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 212 total aligned reads displayed.

GCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCTTC  >  CP009273/3621916‑3622202
                                                                                                                                              |                                                                                                                                                
gCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATagcgcg                                                                                                                                             <  1:2319464/148‑1 (MQ=255)
gCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATagcgcg                                                                                                                                             <  2:2615109/148‑1 (MQ=255)
 cGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAgcgcgc                                                                                                                                            >  2:727350/1‑148 (MQ=255)
 cGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAgcgcgc                                                                                                                                            >  2:2538345/1‑148 (MQ=255)
  gggCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCt                                                                                                                                           <  1:1292775/148‑1 (MQ=255)
  gggCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCt                                                                                                                                           <  2:2969291/148‑1 (MQ=255)
  gggCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCt                                                                                                                                           <  2:2969179/148‑1 (MQ=255)
       gcagTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGAt                                                                                                                                      >  1:2332081/1‑148 (MQ=255)
         agTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATgg                                                                                                                                    <  2:3369940/148‑1 (MQ=255)
           tCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCa                                                                                                                                  >  1:1733552/1‑148 (MQ=255)
           tCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCa                                                                                                                                  >  2:167969/1‑148 (MQ=255)
                tGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAAtctc                                                                                                                             >  1:1764446/1‑148 (MQ=255)
                 gAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATctct                                                                                                                            >  2:198222/1‑148 (MQ=255)
                  aGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTg                                                                                                                           <  2:213271/148‑1 (MQ=255)
                     aTTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGt                                                                                                                        >  1:416774/1‑148 (MQ=255)
                      ttttCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGtt                                                                                                                       <  1:2469951/148‑1 (MQ=255)
                        ttCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGtttt                                                                                                                     >  1:606420/1‑148 (MQ=255)
                        ttCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGtttt                                                                                                                     <  2:185164/148‑1 (MQ=255)
                          cATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCa                                                                                                                   <  2:1807907/148‑1 (MQ=255)
                          cATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCa                                                                                                                   <  2:2519117/148‑1 (MQ=255)
                            tggtggtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAgg                                                                                                                 >  2:3109918/1‑148 (MQ=255)
                             ggtggtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGt                                                                                                                >  1:1310308/1‑148 (MQ=255)
                             ggtggtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGt                                                                                                                <  1:1656378/148‑1 (MQ=255)
                                ggtggATTTCCCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGa                                                                                                             <  2:2720185/148‑1 (MQ=255)
                                ggtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGa                                                                                                             >  1:2477684/1‑148 (MQ=255)
                                  tggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATa                                                                                                           <  1:1044595/148‑1 (MQ=255)
                                     aTTTACCGCAGCCGTTCGAACCAAGAAAACCACAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACgg                                                                                                        <  2:3297882/148‑1 (MQ=255)
                                      tttACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGt                                                                                                       >  1:3153789/1‑148 (MQ=255)
                                       ttACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTg                                                                                                      <  1:567589/148‑1 (MQ=255)
                                          ccGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGa                                                                                                   >  1:1748090/1‑148 (MQ=255)
                                            gCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATc                                                                                                 <  2:1466134/148‑1 (MQ=255)
                                             cAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCa                                                                                                <  2:3449952/148‑1 (MQ=255)
                                             cAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCa                                                                                                >  2:1427895/1‑148 (MQ=255)
                                               gCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACt                                                                                              >  2:698009/1‑148 (MQ=255)
                                                    tCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGc                                                                                         >  2:3061900/1‑148 (MQ=255)
                                                     cGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGcc                                                                                        <  2:3103342/148‑1 (MQ=255)
                                                       aaCCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTg                                                                                      >  1:2180654/1‑148 (MQ=255)
                                                       aaCCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTg                                                                                      >  1:3379114/1‑148 (MQ=255)
                                                        aCCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGa                                                                                     <  2:2857796/148‑1 (MQ=255)
                                                            aGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATgcg                                                                                 >  2:1620108/1‑148 (MQ=255)
                                                              aaaaCCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGcc                                                                               <  1:1802304/148‑1 (MQ=255)
                                                               aaaCCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCt                                                                              <  2:3107833/148‑1 (MQ=255)
                                                                aaCCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCAGGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTg                                                                             >  2:3092018/1‑148 (MQ=255)
                                                                 aCCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGACCCAAAACGCAGGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTgg                                                                            >  2:567382/1‑148 (MQ=255)
                                                                   cAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCg                                                                          >  2:399521/1‑148 (MQ=255)
                                                                   cAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCg                                                                          >  2:932169/1‑148 (MQ=255)
                                                                      aaaTCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTg                                                                       <  2:2529472/148‑1 (MQ=255)
                                                                               cgcgTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCgg                                                              <  1:952810/148‑1 (MQ=255)
                                                                               cgcgTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCgg                                                              <  2:647516/148‑1 (MQ=255)
                                                                               cgcgTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCgg                                                              <  1:1444107/148‑1 (MQ=255)
                                                                                 cgTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTa                                                            >  2:2358525/1‑148 (MQ=255)
                                                                                  gTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTaa                                                           <  2:2474605/148‑1 (MQ=255)
                                                                                  gTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTaa                                                           <  1:2444200/148‑1 (MQ=255)
                                                                                     gAATGCGGAAATTAACGTGATCACCTGCAACGAAGGAACCAAAACGCATGGTCAGATATCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTTATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAg                                                        <  1:3432076/148‑1 (MQ=255)
                                                                                     gAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAg                                                        >  2:2769633/1‑148 (MQ=255)
                                                                                     gAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAg                                                        <  1:2611032/148‑1 (MQ=255)
                                                                                     gAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAg                                                        >  2:1164830/1‑148 (MQ=255)
                                                                                      aaTGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGa                                                       <  1:722232/148‑1 (MQ=255)
                                                                                      aaTGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGa                                                       <  1:2376595/148‑1 (MQ=255)
                                                                                      aaTGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGa                                                       <  1:1368975/148‑1 (MQ=255)
                                                                                       aTGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGAt                                                      <  2:1377423/148‑1 (MQ=255)
                                                                                        tGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGAtt                                                     >  1:2520176/1‑148 (MQ=255)
                                                                                          cGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTta                                                   <  2:3153789/148‑1 (MQ=255)
                                                                                              aaTTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATaaa                                               <  2:3401586/148‑1 (MQ=255)
                                                                                                ttAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATg                                             <  1:1565236/148‑1 (MQ=255)
                                                                                                ttAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATg                                             >  1:2692831/1‑148 (MQ=255)
                                                                                                 tAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGc                                            <  2:1271975/148‑1 (MQ=255)
                                                                                                   aCGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGctt                                          >  2:794144/1‑148 (MQ=255)
                                                                                                   aCGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGctt                                          >  1:1774418/1‑148 (MQ=255)
                                                                                                   aCGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGctt                                          >  2:679338/1‑148 (MQ=255)
                                                                                                    cGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGcttc                                         <  1:1406650/148‑1 (MQ=255)
                                                                                                    cGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGcttc                                         <  2:2813803/148‑1 (MQ=255)
                                                                                                    cGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGcttc                                         <  1:2947335/148‑1 (MQ=255)
                                                                                                       gATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCttcttc                                      >  1:1010887/1‑148 (MQ=255)
                                                                                                       gATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCttcttc                                      >  1:2764473/1‑148 (MQ=255)
                                                                                                       gATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCttcttc                                      >  1:477830/1‑148 (MQ=255)
                                                                                                       gATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCttcttc                                      <  1:2141602/148‑1 (MQ=255)
                                                                                                        aTCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTcc                                     <  1:2528442/148‑1 (MQ=255)
                                                                                                        aTCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTcc                                     <  1:1060024/148‑1 (MQ=255)
                                                                                                            aCGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCg                                 >  2:3084144/1‑148 (MQ=255)
                                                                                                             cGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGt                                >  1:3011042/1‑148 (MQ=255)
                                                                                                                cAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAgc                             >  2:2920936/1‑148 (MQ=255)
                                                                                                                   cGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCt                          <  2:2997303/148‑1 (MQ=255)
                                                                                                                   cGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCt                          >  2:1569605/1‑148 (MQ=255)
                                                                                                                    gAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCtt                         >  1:647496/1‑148 (MQ=255)
                                                                                                                      aGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCACTCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCATCCAGCGTAGCGCTTTg                       >  1:893721/1‑148 (MQ=255)
                                                                                                                       ggAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGc                      >  1:3469721/1‑148 (MQ=255)
                                                                                                                        gAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCg                     <  1:3058492/148‑1 (MQ=255)
                                                                                                                        gAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCg                     >  1:2260807/1‑148 (MQ=255)
                                                                                                                        gAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCg                     >  1:2260985/1‑148 (MQ=255)
                                                                                                                        gAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCg                     <  1:727350/148‑1 (MQ=255)
                                                                                                                        gAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCg                     <  1:2538345/148‑1 (MQ=255)
                                                                                                                         aaCCAAAACTCTTTGTAATATAGCGCGCTTCGATGGCAATATCTGCGTTTTCATTTTGATACGGTGGGATCTCTACCGCCTGATGCGCCTGGCGTTGCGATTGCGGTAACATATTTATAATTGCTTCTTCCAGCGTAGCGCTTTGCGt                    <  2:1795085/148‑1 (MQ=255)
                                                                                                                           ccAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGttt                  >  2:85059/1‑148 (MQ=255)
                                                                                                                              aaaCGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGTGATCACTACCGCCTTATGCGCCTGGCGTTTCGATTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTtgct               >  2:2002711/1‑148 (MQ=255)
                                                                                                                                     tGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAg        <  2:3458082/148‑1 (MQ=255)
                                                                                                                                       gTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGct      <  1:1219568/148‑1 (MQ=255)
                                                                                                                                          aGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCtt   <  1:1595844/148‑1 (MQ=255)
                                                                                                                                          aGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCtt   <  1:3336994/148‑1 (MQ=255)
                                                                                                                                           gATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCTTc  >  1:2010395/1‑148 (MQ=255)
                                                                                                                                              |                                                                                                                                                
GCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCTTC  >  CP009273/3621916‑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: ‑