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,0580CA97.3% 574.3 / 6.5 184D281Y (GAT→TAT) rbbAribosome‑associated ATPase: ATP‑binding protein/ATP‑binding membrane protein
Reads supporting (aligned to +/- strand):  ref base C (3/2);  new base A (91/88);  total (94/90)
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.64e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 186 total aligned reads displayed.

TGGCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGCTCTTC  >  CP009273/3621914‑3622202
                                                                                                                                                |                                                                                                                                                
tGGCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATagcg                                                                                                                                               <  1:325975/148‑1 (MQ=255)
  gCGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATagcgcg                                                                                                                                             <  2:797663/148‑1 (MQ=255)
   cGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAgcgcgc                                                                                                                                            <  2:16907/148‑1 (MQ=255)
   cGGGCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAgcgcgc                                                                                                                                            <  1:2054347/148‑1 (MQ=255)
     ggCAGCAGTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCtt                                                                                                                                          <  1:3165281/148‑1 (MQ=255)
         gcagTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGAt                                                                                                                                      <  2:3358428/148‑1 (MQ=255)
           agTCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATgg                                                                                                                                    <  2:2191402/148‑1 (MQ=255)
             tCCGGTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCa                                                                                                                                  >  1:3048903/1‑148 (MQ=255)
                 gTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAAtct                                                                                                                              <  2:2198388/148‑1 (MQ=255)
                 gTGAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCACGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGAGAGCGCGCTTCGATGGCAAtct                                                                                                                              >  2:1294200/1‑148 (MQ=255)
                 gTGAGCATGTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAAtct                                                                                                                              <  2:2198384/148‑1 (MQ=255)
                   gAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATctct                                                                                                                            >  1:2420155/1‑148 (MQ=255)
                   gAGCATTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCAGGGTCAGATAGCGCGCTTCGATGGCAATctct                                                                                                                            >  1:255619/1‑148 (MQ=255)
                       aTTTTCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGt                                                                                                                        <  2:3186765/148‑1 (MQ=255)
                        ttttCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGtt                                                                                                                       <  1:637366/148‑1 (MQ=255)
                         tttCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGttt                                                                                                                      >  2:550927/1‑148 (MQ=255)
                         tttCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGttt                                                                                                                      <  1:3420234/148‑1 (MQ=255)
                         tttCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGttt                                                                                                                      >  1:2681270/1‑148 (MQ=255)
                          ttCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGtttt                                                                                                                     >  1:373904/1‑148 (MQ=255)
                           tCATGGTGGTGGATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCACTCTCTGCGTTTTc                                                                                                                    >  1:3140702/1‑148 (MQ=255)
                               ggtggtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGt                                                                                                                <  2:773749/148‑1 (MQ=255)
                               ggtggtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGt                                                                                                                <  2:955299/148‑1 (MQ=255)
                               ggtggtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGt                                                                                                                <  2:774585/148‑1 (MQ=255)
                                 tggtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTg                                                                                                              >  1:1004558/1‑148 (MQ=255)
                                  ggtggATTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGa                                                                                                             >  2:1952549/1‑148 (MQ=255)
                                       aTTTACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACgg                                                                                                        <  2:3037933/148‑1 (MQ=255)
                                        tttACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGt                                                                                                       >  2:149247/1‑148 (MQ=255)
                                          tACCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTgg                                                                                                     <  2:398193/148‑1 (MQ=255)
                                           aCCGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTggg                                                                                                    >  1:773071/1‑148 (MQ=255)
                                            ccGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGa                                                                                                   <  2:2543306/148‑1 (MQ=255)
                                            ccGCAGCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGa                                                                                                   <  2:2543487/148‑1 (MQ=255)
                                                 gCCGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACt                                                                                              >  1:866101/1‑148 (MQ=255)
                                                   cGTTCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTAc                                                                                            <  2:2072464/148‑1 (MQ=255)
                                                      tCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGc                                                                                         >  1:988649/1‑148 (MQ=255)
                                                      tCGAACCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGc                                                                                         >  2:331883/1‑148 (MQ=255)
                                                         aaCCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTg                                                                                      <  1:828114/148‑1 (MQ=255)
                                                         aaCCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTg                                                                                      <  1:3067085/148‑1 (MQ=255)
                                                         aaCCAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTg                                                                                      >  2:1776306/1‑148 (MQ=255)
                                                           ccAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGAt                                                                                    >  1:338309/1‑148 (MQ=255)
                                                           ccAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGAt                                                                                    >  1:1451095/1‑148 (MQ=255)
                                                           ccAAGAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCAGGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGAt                                                                                    >  2:655426/1‑148 (MQ=255)
                                                               gAAAACCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATgcgc                                                                                >  2:70784/1‑148 (MQ=255)
                                                                aaaaCCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGcc                                                                               >  1:186898/1‑148 (MQ=255)
                                                                  aaCCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTg                                                                             >  1:2060350/1‑148 (MQ=255)
                                                                  aaCCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTg                                                                             >  1:2124876/1‑148 (MQ=255)
                                                                  aaCCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTg                                                                             >  1:2447813/1‑148 (MQ=255)
                                                                  aaCCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTg                                                                             >  1:2447841/1‑148 (MQ=255)
                                                                   aCCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTgg                                                                            >  2:188465/1‑148 (MQ=255)
                                                                   aCCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTgg                                                                            >  1:219533/1‑148 (MQ=255)
                                                                   aCCAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTgg                                                                            >  2:243633/1‑148 (MQ=255)
                                                                    ccAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGc                                                                           >  2:1135013/1‑148 (MQ=255)
                                                                     cAAAAATCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCg                                                                          <  2:919194/148‑1 (MQ=255)
                                                                      aaaaaTCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGt                                                                         <  2:1581080/148‑1 (MQ=255)
                                                                        aaaTCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTg                                                                       >  1:968838/1‑148 (MQ=255)
                                                                        aaaTCTCCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTg                                                                       >  1:901751/1‑148 (MQ=255)
                                                                            ctcCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCt                                                                   <  2:2442107/148‑1 (MQ=255)
                                                                             tcCCCGCGTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCtt                                                                  <  1:614914/148‑1 (MQ=255)
                                                                                 cgcgTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCgg                                                              >  1:139236/1‑148 (MQ=255)
                                                                                 cgcgTGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCgg                                                              <  1:3231714/148‑1 (MQ=255)
                                                                                     tGGAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAAc                                                          >  2:494517/1‑148 (MQ=255)
                                                                                      ggAATGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACa                                                         >  1:1320437/1‑148 (MQ=255)
                                                                                         aTGCGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGAt                                                      <  2:2577150/148‑1 (MQ=255)
                                                                                           gCGGAACTTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGAttt                                                    <  2:3513216/148‑1 (MQ=255)
                                                                                            cGGAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTta                                                   <  1:471814/148‑1 (MQ=255)
                                                                                              gAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTtata                                                 >  1:2963744/1‑148 (MQ=255)
                                                                                              gAAATTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATGtata                                                 >  2:3445425/1‑148 (MQ=255)
                                                                                               aaaTTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATaa                                                <  2:2647730/148‑1 (MQ=255)
                                                                                                 aTTAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAAt                                              >  1:948143/1‑148 (MQ=255)
                                                                                                  ttAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATg                                             >  1:1188143/1‑148 (MQ=255)
                                                                                                   tAACGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGc                                            <  2:2201555/148‑1 (MQ=255)
                                                                                                    aaCGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGct                                           <  2:362778/148‑1 (MQ=255)
                                                                                                     aCGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGctt                                          <  1:667767/148‑1 (MQ=255)
                                                                                                     aCGTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGctt                                          <  1:667918/148‑1 (MQ=255)
                                                                                                       gTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGcttct                                        <  1:100833/148‑1 (MQ=255)
                                                                                                       gTGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGcttct                                        <  1:493855/148‑1 (MQ=255)
                                                                                                        tGATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGcttctt                                       >  1:1768871/1‑148 (MQ=255)
                                                                                                         gATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCttcttc                                      <  2:2700531/148‑1 (MQ=255)
                                                                                                         gATCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCttcttc                                      <  2:753894/148‑1 (MQ=255)
                                                                                                          aTCAACGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTcc                                     >  1:868726/1‑148 (MQ=255)
                                                                                                             aaCGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGc                                  >  1:3541609/1‑148 (MQ=255)
                                                                                                             aaCGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGc                                  >  1:3541037/1‑148 (MQ=255)
                                                                                                              aCGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCg                                 <  1:2311561/148‑1 (MQ=255)
                                                                                                              aCGGCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCg                                 >  1:1286477/1‑148 (MQ=255)
                                                                                                                 gCAACGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAg                              >  1:1090248/1‑148 (MQ=255)
                                                                                                                   aaCGAAGGAACCAAAACGCATGGTCAGATAGGGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAgcg                            <  2:1832600/148‑1 (MQ=255)
                                                                                                                    aCGAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAgcgc                           <  2:139236/148‑1 (MQ=255)
                                                                                                                      gAAGGAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCtt                         <  2:3281660/148‑1 (MQ=255)
                                                                                                                          gAACCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCTCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCg                     <  2:277339/148‑1 (MQ=255)
                                                                                                                           aaCCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGt                    <  2:528247/148‑1 (MQ=255)
                                                                                                                           aaCCAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGt                    <  2:2821272/148‑1 (MQ=255)
                                                                                                                             ccAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGttt                  >  2:2230916/1‑148 (MQ=255)
                                                                                                                              cAAAACGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTtg                 >  2:3238887/1‑148 (MQ=255)
                                                                                                                               aaaaCGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTtgc                >  1:2032928/1‑148 (MQ=255)
                                                                                                                                 aaCGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTtgctg              >  1:1190380/1‑148 (MQ=255)
                                                                                                                                 aaCGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTtgctg              >  1:2307560/1‑148 (MQ=255)
                                                                                                                                  aCGCATGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTtgctgc             <  1:1068914/148‑1 (MQ=255)
                                                                                                                                     cATGGTCAGATCGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGt          <  2:939139/148‑1 (MQ=255)
                                                                                                                                       tGGTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAg        <  2:3510723/148‑1 (MQ=255)
                                                                                                                                         gTCAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGct      >  1:1257641/1‑148 (MQ=255)
                                                                                                                                           cAGATAGCGCGCTTCGATGGCAATCTCTGCGTTTTCAGGTTGATACGGTGGGATCACTACCGCCTGATGCGCCTGGCGTTGCGCTTGCGGTAACAGATTTATAAATGCTTCTTCCAGCGTAGCGCTTTGCGTTTGCTGCCGTAGctct    <  1:2778534/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: