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,1790AT97.8% 745.4 / 4.2 228pseudogene (150/282 nt)yaiXpseudogene, interrupted by IS2A, acetyltransferase homolog; nonfunctional; interruped by IS2; putative transferase
Reads supporting (aligned to +/- strand):  ref base A (4/1);  new base T (97/126);  total (101/127)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.84e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 241 total aligned reads displayed.

AGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATAAG  >  CP009273/378034‑378323
                                                                                                                                                 |                                                                                                                                                
aGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATtttt                                                                                                                                               <  1:3375180/149‑1 (MQ=255)
 ggggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAAtttttt                                                                                                                                              >  2:887227/1‑149 (MQ=255)
  gggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATtttttt                                                                                                                                             >  2:1513786/1‑149 (MQ=255)
  gggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATtttttt                                                                                                                                             >  2:1512939/1‑149 (MQ=255)
  gggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATtttttt                                                                                                                                             <  2:2861629/149‑1 (MQ=255)
  gggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATtttttt                                                                                                                                             >  2:3669763/1‑149 (MQ=255)
   ggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTtc                                                                                                                                            <  1:3183658/149‑1 (MQ=255)
     cAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTtctc                                                                                                                                          >  1:4047768/1‑149 (MQ=255)
      aaaTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTctct                                                                                                                                         <  2:3612669/149‑1 (MQ=255)
      aaaTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTctct                                                                                                                                         <  1:496423/149‑1 (MQ=255)
       aaTCCAATTACCTATCAGGCATTTTGCGACTATAACTACCTGACCCTGTATAACGGCAACATGTGAATTACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATTATTACCTGGTCTGCAATCTGAATTTTTTTCTCtt                                                                                                                                        <  2:699694/149‑1 (MQ=255)
       aaTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCtt                                                                                                                                        <  1:3817431/149‑1 (MQ=255)
       aaTCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCtt                                                                                                                                        <  1:345569/149‑1 (MQ=255)
          ccAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTaa                                                                                                                                     <  2:1937410/149‑1 (MQ=255)
           cAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAAc                                                                                                                                    <  1:3563385/149‑1 (MQ=255)
            aagtACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACa                                                                                                                                   <  1:1398062/146‑1 (MQ=255)
            aaTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACa                                                                                                                                   <  1:1313220/149‑1 (MQ=255)
            aaTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACa                                                                                                                                   <  1:2190016/149‑1 (MQ=255)
            aaTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACa                                                                                                                                   <  2:2208959/149‑1 (MQ=255)
            aaTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACa                                                                                                                                   <  1:1664592/149‑1 (MQ=255)
            aaTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACa                                                                                                                                   <  2:3745217/149‑1 (MQ=255)
             aTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACAt                                                                                                                                  >  1:2003147/1‑149 (MQ=255)
               tACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTa                                                                                                                                >  2:343887/1‑149 (MQ=255)
                aCCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTaa                                                                                                                               >  1:647931/1‑149 (MQ=255)
                aCCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTaa                                                                                                                               >  1:2714215/1‑149 (MQ=255)
                aCCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTaa                                                                                                                               >  2:1145761/1‑149 (MQ=255)
                 ccTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAg                                                                                                                              <  1:2382515/149‑1 (MQ=255)
                 ccTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCAACATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAg                                                                                                                              <  2:931911/149‑1 (MQ=255)
                   tATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCt                                                                                                                            <  2:3026944/149‑1 (MQ=255)
                    aTCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCtt                                                                                                                           <  2:1014851/149‑1 (MQ=255)
                     tCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTc                                                                                                                          >  2:1742171/1‑149 (MQ=255)
                       aGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAg                                                                                                                        >  2:755632/1‑149 (MQ=255)
                        gtcAGTTTGCGCCAATCACTACCGGACACTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGa                                                                                                                       <  1:3087196/147‑1 (MQ=255)
                        ggCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGa                                                                                                                       <  1:1869643/149‑1 (MQ=255)
                        ggCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGa                                                                                                                       <  1:4009114/149‑1 (MQ=255)
                          cAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATa                                                                                                                     <  2:3978157/149‑1 (MQ=255)
                           aGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAg                                                                                                                    <  2:1897901/149‑1 (MQ=255)
                             tttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGtt                                                                                                                  >  2:3489585/1‑149 (MQ=255)
                              ttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGttt                                                                                                                 <  1:3278205/149‑1 (MQ=255)
                               tGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGtttt                                                                                                                <  2:2052087/149‑1 (MQ=255)
                                  gcCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCa                                                                                                             <  1:1062172/149‑1 (MQ=255)
                                  gcCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCa                                                                                                             <  2:352344/149‑1 (MQ=255)
                                    cAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGt                                                                                                           >  2:1544065/1‑149 (MQ=255)
                                     aaTCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGtt                                                                                                          >  1:3502590/1‑149 (MQ=255)
                                     aaTCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGtt                                                                                                          <  2:458279/149‑1 (MQ=255)
                                         aCTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCtg                                                                                                      >  1:3141587/1‑149 (MQ=255)
                                          cTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCtgt                                                                                                     <  1:1992481/149‑1 (MQ=255)
                                            aCCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGa                                                                                                   >  1:1273665/1‑149 (MQ=255)
                                                gACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGa                                                                                               >  1:250231/1‑149 (MQ=255)
                                                gACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGa                                                                                               >  1:1408274/1‑149 (MQ=255)
                                                    cTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGt                                                                                           <  1:3956389/149‑1 (MQ=255)
                                                    cTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGt                                                                                           <  1:173940/149‑1 (MQ=255)
                                                       aaTAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTgg                                                                                        >  2:3375373/1‑149 (MQ=255)
                                                         tAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGa                                                                                      >  1:986176/1‑149 (MQ=255)
                                                         tAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGa                                                                                      >  1:2679186/1‑149 (MQ=255)
                                                              gCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAg                                                                                 <  1:122719/149‑1 (MQ=255)
                                                                aCCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCg                                                                               >  1:926130/1‑149 (MQ=255)
                                                                 ccATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCgg                                                                              >  2:3277052/1‑149 (MQ=255)
                                                                  cATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGt                                                                             <  2:4001961/149‑1 (MQ=255)
                                                                   aTGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGtt                                                                            >  2:1966519/1‑149 (MQ=255)
                                                                     gACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGa                                                                          >  1:1928323/1‑149 (MQ=255)
                                                                      aCAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGAc                                                                         <  2:522049/149‑1 (MQ=255)
                                                                       cAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACg                                                                        >  1:446859/1‑149 (MQ=255)
                                                                        aaaTACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGc                                                                       <  2:1351057/149‑1 (MQ=255)
                                                                           tACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTc                                                                    <  2:676797/149‑1 (MQ=255)
                                                                               ggTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTccacc                                                                >  1:849999/1‑149 (MQ=255)
                                                                                   ttCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTg                                                            >  1:1621502/1‑149 (MQ=255)
                                                                                    tCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGc                                                           <  1:2072230/149‑1 (MQ=255)
                                                                                        gcCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGaa                                                       <  2:1360807/149‑1 (MQ=255)
                                                                                           gATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAg                                                    >  2:2444042/1‑149 (MQ=255)
                                                                                           gATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAg                                                    <  1:3707129/149‑1 (MQ=255)
                                                                                                         cTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTc                                      >  1:580611/1‑149 (MQ=255)
                                                                                                             aGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCgcgc                                  <  1:1031890/149‑1 (MQ=255)
                                                                                                              gTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCt                                 >  2:1286540/1‑149 (MQ=255)
                                                                                                                ttCATCAATAATTCCCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGa                               <  1:3606704/149‑1 (MQ=255)
                                                                                                                 tcatcaATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAg                              >  1:2822794/1‑149 (MQ=255)
                                                                                                                  catcaATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAgg                             >  1:3353340/1‑149 (MQ=255)
                                                                                                                   atcaATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGt                            <  2:333273/149‑1 (MQ=255)
                                                                                                                   atcaATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGt                            >  1:2757362/1‑149 (MQ=255)
                                                                                                                    tcaATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGtt                           >  1:579176/1‑149 (MQ=255)
                                                                                                                    tcaATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGtt                           >  1:1521998/1‑149 (MQ=255)
                                                                                                                      aataatTACCTGGTCTGCAATCGGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTcg                         >  2:2962242/1‑149 (MQ=255)
                                                                                                                         aatTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGt                      >  2:159578/1‑149 (MQ=255)
                                                                                                                            tACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCa                   <  1:1156114/149‑1 (MQ=255)
                                                                                                                            tACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCa                   <  2:2849621/149‑1 (MQ=255)
                                                                                                                            tACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCa                   <  1:1782471/149‑1 (MQ=255)
                                                                                                                              ccTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGt                 <  1:2613693/149‑1 (MQ=255)
                                                                                                                               cTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTa                <  2:1429307/149‑1 (MQ=255)
                                                                                                                                   tCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTaaaaa            >  1:2556424/1‑149 (MQ=255)
                                                                                                                                   tCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTaaaaa            >  2:3368437/1‑149 (MQ=255)
                                                                                                                                     tGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAg          >  1:1701011/1‑149 (MQ=255)
                                                                                                                                      gCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAgg         <  1:864618/149‑1 (MQ=255)
                                                                                                                                       cAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGc        >  2:2667824/1‑149 (MQ=255)
                                                                                                                                       cAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGc        <  1:392477/149‑1 (MQ=255)
                                                                                                                                         aTCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCaa      <  1:1544065/149‑1 (MQ=255)
                                                                                                                                         aTCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCaa      >  1:2904542/1‑149 (MQ=255)
                                                                                                                                           cTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATa    <  2:2822794/149‑1 (MQ=255)
                                                                                                                                           cTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATa    <  2:2349149/149‑1 (MQ=255)
                                                                                                                                            tGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATaa   <  1:2653188/149‑1 (MQ=255)
                                                                                                                                             gAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATAAg  >  2:1925048/1‑149 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
AGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATAAG  >  CP009273/378034‑378323

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: