Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 378,179 A→T 100% pseudogene (150/282 nt) yaiX ← pseudogene, interrupted by IS2A, acetyltransferase homolog; nonfunctional; interruped by IS2; putative transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273378,1790AT96.1% 580.5 / 10.3 181pseudogene (150/282 nt)yaiXpseudogene, interrupted by IS2A, acetyltransferase homolog; nonfunctional; interruped by IS2; putative transferase
Reads supporting (aligned to +/- strand):  ref base A (4/3);  new base T (82/92);  total (86/95)
Fisher's exact test for biased strand distribution p-value = 7.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 191 total aligned reads displayed.

AGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATAA  >  CP009273/378034‑378322
                                                                                                                                                 |                                                                                                                                               
aGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAAtttt                                                                                                                                               >  1:884704/1‑148 (MQ=255)
 ggggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATtttt                                                                                                                                              <  1:808842/148‑1 (MQ=255)
 ggggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATtttt                                                                                                                                              <  2:535893/148‑1 (MQ=255)
  gggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAAtttttt                                                                                                                                             <  1:1386976/148‑1 (MQ=255)
   ggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATtttttt                                                                                                                                            <  1:2241679/148‑1 (MQ=255)
   ggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATtttttt                                                                                                                                            <  2:1709039/148‑1 (MQ=255)
      aaaTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTtctc                                                                                                                                         <  1:2902596/148‑1 (MQ=255)
      aaaTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTtctc                                                                                                                                         <  1:112498/148‑1 (MQ=255)
        aTCCAATTACCTATCAGGCAGTTTGCGCCTATCACGACCGTACCCTGAATATCGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCAGCTGCAGTTTCATCAATAATTACCTGGTCTGCTATCTGAATTTTTTTCTCtt                                                                                                                                       <  2:616115/148‑1 (MQ=255)
        aTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCtt                                                                                                                                       <  1:2336101/148‑1 (MQ=255)
        aTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCtt                                                                                                                                       <  1:770912/148‑1 (MQ=255)
         tCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCttt                                                                                                                                      >  1:1552767/1‑148 (MQ=255)
          ccAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTa                                                                                                                                     <  2:1825468/148‑1 (MQ=255)
           cAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTaa                                                                                                                                    >  1:2784471/1‑148 (MQ=255)
            aaTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAAc                                                                                                                                   <  2:1007130/148‑1 (MQ=255)
              ttACCTATCAGGGAGTTTGCGCCAATCTGTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGACGATAACTACTTCTCCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACAt                                                                                                                                 <  1:1609767/148‑1 (MQ=255)
              ttACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACAt                                                                                                                                 <  1:1269595/148‑1 (MQ=255)
               tACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACAtt                                                                                                                                >  1:2918765/1‑148 (MQ=255)
               tACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACAtt                                                                                                                                >  2:364176/1‑148 (MQ=255)
                aCCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTa                                                                                                                               <  2:691435/148‑1 (MQ=255)
                         gCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGa                                                                                                                      >  1:1026972/1‑148 (MQ=255)
                         gCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGa                                                                                                                      >  1:1026980/1‑148 (MQ=255)
                           aGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATa                                                                                                                    <  2:634636/148‑1 (MQ=255)
                            gTTTGCGACAATCACTACCTGACACTGTATAACGGCAACATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAg                                                                                                                   <  1:1302431/148‑1 (MQ=255)
                             tttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGt                                                                                                                  <  2:468608/148‑1 (MQ=255)
                              ttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGtt                                                                                                                 <  2:231607/148‑1 (MQ=255)
                              ttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGtt                                                                                                                 <  1:1459837/148‑1 (MQ=255)
                               tGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGttt                                                                                                                <  1:2175084/148‑1 (MQ=255)
                               tGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGttt                                                                                                                >  2:2223646/1‑148 (MQ=255)
                                gcgcCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGtttt                                                                                                               >  1:2884463/1‑148 (MQ=255)
                                 cgcCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGttttt                                                                                                              <  1:2461978/148‑1 (MQ=255)
                                 cgcCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGttttt                                                                                                              <  1:3057594/148‑1 (MQ=255)
                                  gcCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTc                                                                                                             <  1:1752145/148‑1 (MQ=255)
                                     aaTCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGt                                                                                                          >  1:797930/1‑148 (MQ=255)
                                       tCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTc                                                                                                        >  2:2729439/1‑148 (MQ=255)
                                       tCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTc                                                                                                        <  2:1592377/148‑1 (MQ=255)
                                       tCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTc                                                                                                        >  2:1123962/1‑148 (MQ=255)
                                       tCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTc                                                                                                        >  2:1006143/1‑148 (MQ=255)
                                                gACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAg                                                                                               >  2:1774083/1‑148 (MQ=255)
                                                  cccTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGaa                                                                                             <  2:1077379/148‑1 (MQ=255)
                                                   ccTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAg                                                                                            <  2:435671/148‑1 (MQ=255)
                                                       aaTAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCTGTTCCTGTGACAGAAGGTTg                                                                                        >  1:3213250/1‑148 (MQ=255)
                                                           aCGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGaa                                                                                    >  1:3021536/1‑148 (MQ=255)
                                                           aCGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGaa                                                                                    <  2:2119582/148‑1 (MQ=255)
                                                           aCGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGaa                                                                                    <  1:2075240/148‑1 (MQ=255)
                                                              gCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATaa                                                                                 >  1:711658/1‑148 (MQ=255)
                                                              gCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATaa                                                                                 <  1:3008424/148‑1 (MQ=255)
                                                                 ccATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCg                                                                              >  2:2122717/1‑148 (MQ=255)
                                                                 ccATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCg                                                                              >  1:60056/1‑148 (MQ=255)
                                                                 ccATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCg                                                                              >  1:705284/1‑148 (MQ=255)
                                                                 ccATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCg                                                                              >  1:235930/1‑148 (MQ=255)
                                                                   aTGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGt                                                                            <  1:1989445/148‑1 (MQ=255)
                                                                     gACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTg                                                                          <  1:2945743/148‑1 (MQ=255)
                                                                       cAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGAc                                                                        >  2:1645541/1‑148 (MQ=255)
                                                                        aaaTACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACg                                                                       >  1:1427705/1‑148 (MQ=255)
                                                                         aaTACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGc                                                                      <  2:1963794/148‑1 (MQ=255)
                                                                          aTACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCt                                                                     >  2:1397992/1‑148 (MQ=255)
                                                                               ggTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTccac                                                                >  1:3049590/1‑148 (MQ=255)
                                                                               ggTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTccac                                                                >  1:2774434/1‑148 (MQ=255)
                                                                                    tCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTg                                                           <  2:3100372/148‑1 (MQ=255)
                                                                                        gcCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGa                                                       <  2:935414/148‑1 (MQ=255)
                                                                                        gcCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGa                                                       <  1:1644689/148‑1 (MQ=255)
                                                                                         ccGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGaa                                                      <  2:2975873/148‑1 (MQ=255)
                                                                                         ccGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGaa                                                      <  1:376425/148‑1 (MQ=255)
                                                                                            aTAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAg                                                   <  2:2229623/148‑1 (MQ=255)
                                                                                             tAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGc                                                  >  2:1946386/1‑148 (MQ=255)
                                                                                               acaacTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCgg                                                >  2:244191/1‑148 (MQ=255)
                                                                                                   cTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCTGTTGACGCTTCAACCAGTGCGGAATAGCgggtgg                                            >  1:2380472/1‑148 (MQ=255)
                                                                                                   cTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCgggtgg                                            >  2:2842414/1‑148 (MQ=255)
                                                                                                   cTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCgggtgg                                            >  2:311244/1‑148 (MQ=255)
                                                                                                      cACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGtt                                         <  2:2213602/148‑1 (MQ=255)
                                                                                                      cACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGtt                                         <  2:2213632/148‑1 (MQ=255)
                                                                                                       aCCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGttt                                        >  1:1318115/1‑148 (MQ=255)
                                                                                                              gTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCgcgc                                 >  2:2264224/1‑148 (MQ=255)
                                                                                                                ttCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTg                               >  2:188337/1‑148 (MQ=255)
                                                                                                                   atcaATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAgg                            <  2:1552767/148‑1 (MQ=255)
                                                                                                                    tcaATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGt                           >  2:3024135/1‑148 (MQ=255)
                                                                                                                    tcaATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGt                           <  2:884704/148‑1 (MQ=255)
                                                                                                                    tcaATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGt                           >  2:1948021/1‑148 (MQ=255)
                                                                                                                     caATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGtt                          >  2:1454650/1‑148 (MQ=255)
                                                                                                                       ataatTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTcg                        >  1:270115/1‑148 (MQ=255)
                                                                                                                        taatTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTcgc                       >  1:2739132/1‑148 (MQ=255)
                                                                                                                         aatTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTcgcg                      >  1:1919124/1‑148 (MQ=255)
                                                                                                                          atTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGt                     >  2:1389753/1‑148 (MQ=255)
                                                                                                                           ttACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTc                    >  1:1193468/1‑148 (MQ=255)
                                                                                                                               cTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGt                <  1:1107993/148‑1 (MQ=255)
                                                                                                                               cTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGt                <  1:460021/148‑1 (MQ=255)
                                                                                                                                 ggTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTaa              >  1:569659/1‑148 (MQ=255)
                                                                                                                                   tCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTaaaa            <  2:2602540/148‑1 (MQ=255)
                                                                                                                                      gCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAg         <  1:2223646/148‑1 (MQ=255)
                                                                                                                                       cAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAgg        >  2:1518691/1‑148 (MQ=255)
                                                                                                                                         aTCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCa      <  2:2370988/148‑1 (MQ=255)
                                                                                                                                         aTCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCa      <  1:2842414/148‑1 (MQ=255)
                                                                                                                                         aTCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCa      <  2:711658/148‑1 (MQ=255)
                                                                                                                                         aTCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCa      >  2:550361/1‑148 (MQ=255)
                                                                                                                                          tCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCaa     <  1:944256/148‑1 (MQ=255)
                                                                                                                                          tCTGAATTTTTTACTCTATAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCaa     <  1:1740304/148‑1 (MQ=255)
                                                                                                                                            tGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTCAAAAAGGCAATa   >  2:1699718/1‑148 (MQ=255)
                                                                                                                                            tGAATTTTTGTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATa   <  2:3311921/148‑1 (MQ=255)
                                                                                                                                             gAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATaa  <  2:2632920/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
AGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATAA  >  CP009273/378034‑378322

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: