Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 726,087 G→T 100% R350L (CGT→CTT)  rhsC → Rhs family putative polymorphic toxin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273726,0870GT98.8% 555.2 / ‑1.8 171R350L (CGT→CTT) rhsCRhs family putative polymorphic toxin
Reads supporting (aligned to +/- strand):  ref base G (2/0);  new base T (87/82);  total (89/82)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-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 171 total aligned reads displayed.

TACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAG  >  CP009273/725942‑726230
                                                                                                                                                 |                                                                                                                                               
tACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTc                                                                                                                                               <  1:824324/148‑1 (MQ=255)
tACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTc                                                                                                                                               <  1:2870781/148‑1 (MQ=255)
tACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTc                                                                                                                                               <  2:2828577/148‑1 (MQ=255)
 aCCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTca                                                                                                                                              <  1:1662138/148‑1 (MQ=255)
  cccGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGACCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTcac                                                                                                                                             <  1:1616924/148‑1 (MQ=255)
   ccGGAGAATTTACCGGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTcaca                                                                                                                                            <  1:1512399/148‑1 (MQ=255)
     ggAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACg                                                                                                                                          >  2:2757624/1‑148 (MQ=255)
     ggAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACg                                                                                                                                          >  2:2757670/1‑148 (MQ=255)
       agaATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACggg                                                                                                                                        >  2:851299/1‑148 (MQ=255)
         aaTTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGcc                                                                                                                                      >  2:2628029/1‑148 (MQ=255)
         aaTTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGcc                                                                                                                                      >  1:1337598/1‑148 (MQ=255)
         aaTTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGcc                                                                                                                                      >  2:2627893/1‑148 (MQ=255)
           tttCCCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGCTGATAAATACCGGGGCCGGAGGGTGGCGCACCTTCACACGGGCCGa                                                                                                                                    <  1:1236495/148‑1 (MQ=255)
           tttACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGa                                                                                                                                    >  1:2701189/1‑148 (MQ=255)
               ccTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCgg                                                                                                                                <  1:2340470/148‑1 (MQ=255)
                 tGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGaa                                                                                                                              >  2:1023168/1‑148 (MQ=255)
                 tGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGaa                                                                                                                              >  2:1023180/1‑148 (MQ=255)
                   ccgcgccgcTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAAt                                                                                                                            <  1:3000916/148‑1 (MQ=255)
                     gcgccgcTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCt                                                                                                                          <  1:2124386/148‑1 (MQ=255)
                     gcgccgcTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCt                                                                                                                          <  1:2121747/148‑1 (MQ=255)
                      cgccgcTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTg                                                                                                                         >  2:2140308/1‑148 (MQ=255)
                       gccgcTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGt                                                                                                                        >  1:1973818/1‑148 (MQ=255)
                             ggTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGt                                                                                                                  <  1:1447266/148‑1 (MQ=255)
                             ggTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGt                                                                                                                  <  2:1024101/148‑1 (MQ=255)
                               tGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTa                                                                                                                >  1:950194/1‑148 (MQ=255)
                               tGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTa                                                                                                                >  2:3230260/1‑148 (MQ=255)
                               tGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTa                                                                                                                >  2:3229531/1‑148 (MQ=255)
                                  gcTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGa                                                                                                             <  2:2530163/148‑1 (MQ=255)
                                    tATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACa                                                                                                           >  1:1915823/1‑148 (MQ=255)
                                     aTGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAg                                                                                                          >  2:3506279/1‑148 (MQ=255)
                                     aTGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAg                                                                                                          <  1:600654/148‑1 (MQ=255)
                                        gCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGGGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGa                                                                                                       <  1:680137/148‑1 (MQ=255)
                                             aCGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGAcgggc                                                                                                  >  1:2083915/1‑148 (MQ=255)
                                                  ccGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTg                                                                                             >  2:1578741/1‑148 (MQ=255)
                                                        cgAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGaa                                                                                       >  2:595905/1‑148 (MQ=255)
                                                                   cggTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAAccc                                                                            >  2:1473681/1‑148 (MQ=255)
                                                                     gtgtATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCgg                                                                          >  1:70116/1‑148 (MQ=255)
                                                                        tATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAg                                                                       >  2:1763454/1‑148 (MQ=255)
                                                                         aTGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAgg                                                                      >  2:1103625/1‑148 (MQ=255)
                                                                          tGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGc                                                                     <  1:2934734/148‑1 (MQ=255)
                                                                          tGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGc                                                                     <  2:2677674/148‑1 (MQ=255)
                                                                          tGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGc                                                                     >  2:2896251/1‑148 (MQ=255)
                                                                            aCCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCtt                                                                   >  2:1492679/1‑148 (MQ=255)
                                                                             ccGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTa                                                                  <  1:2101542/148‑1 (MQ=255)
                                                                                cagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGt                                                               <  2:381546/148‑1 (MQ=255)
                                                                                cagcaATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGt                                                               <  2:1073491/148‑1 (MQ=255)
                                                                                    aaTACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTaca                                                           <  2:376084/148‑1 (MQ=255)
                                                                                     aTACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTacac                                                          >  1:1823581/1‑148 (MQ=255)
                                                                                       acacAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGt                                                        >  2:183365/1‑148 (MQ=38)
                                                                                       acacAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGt                                                        >  2:2968671/1‑148 (MQ=38)
                                                                                        cacaGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTa                                                       >  2:1734766/1‑148 (MQ=38)
                                                                                          caGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATc                                                     >  1:3196545/1‑148 (MQ=37)
                                                                                          caGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATc                                                     >  1:1451827/1‑148 (MQ=37)
                                                                                           aGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCa                                                    <  2:622402/148‑1 (MQ=37)
                                                                                            ggTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAg                                                   <  1:3305845/148‑1 (MQ=37)
                                                                                            ggTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAg                                                   <  1:3306579/148‑1 (MQ=37)
                                                                                                cgcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATg                                               <  1:2236091/148‑1 (MQ=37)
                                                                                                cgcAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATg                                               <  1:131435/148‑1 (MQ=37)
                                                                                                  catcTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCTACCGGAAATCTTTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATgag                                             <  1:1198934/145‑1 (MQ=255)
                                                                                                  cAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATgag                                             >  1:2382132/1‑148 (MQ=37)
                                                                                                  cAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATgag                                             <  1:2475330/148‑1 (MQ=37)
                                                                                                   aGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATgaga                                            <  2:2232581/148‑1 (MQ=37)
                                                                                                   aGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATgaga                                            <  2:2232582/148‑1 (MQ=37)
                                                                                                   aGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATgaga                                            <  2:844270/148‑1 (MQ=37)
                                                                                                    gCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGaa                                           <  1:895526/148‑1 (MQ=37)
                                                                                                    gCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGaa                                           <  1:895462/148‑1 (MQ=37)
                                                                                                      tttACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAgg                                         >  2:1991401/1‑148 (MQ=38)
                                                                                                      tttACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAgg                                         <  1:1023180/148‑1 (MQ=38)
                                                                                                      tttACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAgg                                         <  2:1777503/148‑1 (MQ=38)
                                                                                                      tttACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAgg                                         <  1:1023168/148‑1 (MQ=38)
                                                                                                        tacttacGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGAc                                       <  1:3407124/148‑1 (MQ=255)
                                                                                                          cttacGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCg                                     <  1:1382746/148‑1 (MQ=255)
                                                                                                           ttacGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGc                                    >  1:2512237/1‑148 (MQ=255)
                                                                                                           ttacGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGc                                    <  1:2838558/148‑1 (MQ=255)
                                                                                                            tacGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGca                                   >  2:3300656/1‑148 (MQ=255)
                                                                                                            tacGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGca                                   >  1:378020/1‑148 (MQ=255)
                                                                                                            tacGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGca                                   >  1:378013/1‑148 (MQ=255)
                                                                                                              cGATGATAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatc                                 >  2:2295545/1‑148 (MQ=255)
                                                                                                                 tgatAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcacc                              >  2:1205167/1‑148 (MQ=255)
                                                                                                                 tgatAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcacc                              >  2:2955166/1‑148 (MQ=255)
                                                                                                                 tgatAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcacc                              <  1:2589944/148‑1 (MQ=255)
                                                                                                                 tgatAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcacc                              <  2:2770096/148‑1 (MQ=255)
                                                                                                                  gatAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcacca                             >  2:928782/1‑148 (MQ=255)
                                                                                                                  gatAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcacca                             >  1:1153860/1‑148 (MQ=255)
                                                                                                                  gatAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcacca                             >  2:2467462/1‑148 (MQ=255)
                                                                                                                   atAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccat                            >  2:246148/1‑148 (MQ=255)
                                                                                                                    tAAATACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccatc                           >  1:543249/1‑148 (MQ=255)
                                                                                                                     aaaTACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccatca                          <  1:2627893/148‑1 (MQ=255)
                                                                                                                     aaaTACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccatca                          <  1:2628053/148‑1 (MQ=255)
                                                                                                                     aaaTACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGcatcaccatca                          >  2:131746/1‑148 (MQ=255)
                                                                                                                        tACCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCg                       >  1:1546398/1‑148 (MQ=255)
                                                                                                                         aCCGGGGCCGGATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGa                      >  1:2135547/1‑148 (MQ=255)
                                                                                                                                   gATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAc            >  2:2480628/1‑148 (MQ=255)
                                                                                                                                   gATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAc            >  1:2330141/1‑148 (MQ=255)
                                                                                                                                   gATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAc            <  1:2048065/148‑1 (MQ=255)
                                                                                                                                   gATGGTGGCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAc            >  2:1891784/1‑148 (MQ=255)
                                                                                                                                      ggtggCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAccgc         >  1:1710953/1‑148 (MQ=255)
                                                                                                                                      ggtggCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAccgc         <  1:91513/148‑1 (MQ=255)
                                                                                                                                       gtggCGCACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAAccgcc        >  1:1138726/1‑148 (MQ=255)
                                                                                                                                            gcACCTTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGaa   <  1:1293059/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
TACCCGGAGAATTTACCTGCCGCGCCGCTGGTGCGCTATGGCTGGACGCCCCGCGGCGAACTGGCGGCGGTGTATGACCGCAGCAATACACAGGTGCGCAGCTTTACTTACGATGATAAATACCGGGGCCGGATGGTGGCGCACCGTCACACGGGCCGACCGGAAATCTGTTACCGTTACGACAGCGACGGGCGGGTGACAGAACAGCTAAACCCGGCAGGCTTAAGTTACACGTATCAGTATGAGAAGGACCGCATCACCATCACCGACAGCCTGAACCGCCGTGAAG  >  CP009273/725942‑726230

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: ‑