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,1790AT99.0% 657.4 / ‑2.5 202pseudogene (150/282 nt)yaiXpseudogene, interrupted by IS2A, acetyltransferase homolog; nonfunctional; interruped by IS2; putative transferase
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base T (98/102);  total (99/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-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.
Rejected as polymorphism: Polymorphic base substitution creates a homopolymer stretch.

Only 100 of 213 total aligned reads displayed.

AGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATAAG  >  CP009273/378034‑378323
                                                                                                                                                 |                                                                                                                                                
aGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAAtttt                                                                                                                                                <  2:1829072/148‑1 (MQ=255)
aGGGGCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAAtttt                                                                                                                                                >  1:638257/1‑148 (MQ=255)
 ggggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATtttt                                                                                                                                               >  1:2760019/1‑148 (MQ=255)
 ggggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATtttt                                                                                                                                               <  1:2962403/148‑1 (MQ=255)
 ggggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATtttt                                                                                                                                               <  1:3251373/148‑1 (MQ=255)
  gggCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAAtttttt                                                                                                                                              <  2:429573/148‑1 (MQ=255)
    gCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTtc                                                                                                                                            >  2:1301713/1‑148 (MQ=255)
    gCAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTtc                                                                                                                                            <  2:2333983/148‑1 (MQ=255)
     cAAATCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTtct                                                                                                                                           <  1:1796750/148‑1 (MQ=255)
         tCCAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCttt                                                                                                                                       >  2:1753100/1‑148 (MQ=255)
          ccAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTa                                                                                                                                      <  2:745510/148‑1 (MQ=255)
          ccAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTa                                                                                                                                      <  2:2164529/148‑1 (MQ=255)
           cAATTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTaa                                                                                                                                     <  1:1899144/148‑1 (MQ=255)
            aaTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAAc                                                                                                                                    >  1:1845952/1‑148 (MQ=255)
            aaTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAAc                                                                                                                                    <  2:2819809/148‑1 (MQ=255)
             aTTACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCTATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACa                                                                                                                                   <  1:3074442/148‑1 (MQ=255)
               tACCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACAtt                                                                                                                                 <  2:1932004/148‑1 (MQ=255)
                aCCTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTa                                                                                                                                <  2:314588/148‑1 (MQ=255)
                 ccTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTaa                                                                                                                               >  2:3299760/1‑148 (MQ=255)
                 ccTATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACATGCAGTTTCATCAATAATTACCTGGTCTGCAATATGAATTTTTTTCTCTTTAACATTaa                                                                                                                               >  1:1816977/1‑148 (MQ=255)
                   tATCAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGc                                                                                                                             >  1:1342054/1‑148 (MQ=255)
                      cAGGCAGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTc                                                                                                                          <  1:2427920/148‑1 (MQ=255)
                           aGTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATa                                                                                                                     <  2:2954615/148‑1 (MQ=255)
                            gTTTGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAg                                                                                                                    >  2:76449/1‑148 (MQ=255)
                              ttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGtt                                                                                                                  <  1:754778/148‑1 (MQ=255)
                              ttGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATATTTTTCTCTTTAACATTAAGCTTCAGATAGtt                                                                                                                  <  2:1986371/148‑1 (MQ=255)
                               tGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGttt                                                                                                                 <  2:2415145/148‑1 (MQ=255)
                               tGCGCCAATCACTACCGGACCCTGAATAACGGCACCATGACAAAGACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGttt                                                                                                                 <  2:2072816/148‑1 (MQ=255)
                                gcgcCAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGtttt                                                                                                                <  1:1996098/148‑1 (MQ=255)
                                    cAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAg                                                                                                            <  2:1400844/148‑1 (MQ=255)
                                    cAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAg                                                                                                            <  2:1789618/148‑1 (MQ=255)
                                    cAATCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAg                                                                                                            <  2:1080986/148‑1 (MQ=255)
                                     aaTCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGt                                                                                                           <  2:2012852/148‑1 (MQ=255)
                                       tCACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTc                                                                                                         <  2:1188680/148‑1 (MQ=255)
                                        cACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTcc                                                                                                        <  2:2268327/148‑1 (MQ=255)
                                        cACTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTcc                                                                                                        <  2:2205470/148‑1 (MQ=255)
                                          cTACCGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCtg                                                                                                      <  1:2989995/148‑1 (MQ=255)
                                             ccGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGa                                                                                                   <  2:1764816/148‑1 (MQ=255)
                                             ccGGACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGa                                                                                                   >  1:1981080/1‑148 (MQ=255)
                                                gACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAg                                                                                                >  1:3134636/1‑148 (MQ=255)
                                                gACCCTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAg                                                                                                >  1:1048999/1‑148 (MQ=255)
                                                  cccTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGaa                                                                                              <  2:2320397/148‑1 (MQ=255)
                                                  cccTGAATAACGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGaa                                                                                              <  1:1333040/148‑1 (MQ=255)
                                                            cGGCACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAAt                                                                                    >  2:1401826/1‑148 (MQ=255)
                                                               cACCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAg                                                                                 <  2:1584525/148‑1 (MQ=255)
                                                                aCCATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGc                                                                                <  1:3422700/148‑1 (MQ=255)
                                                                  cATGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCgg                                                                              >  2:3227156/1‑148 (MQ=255)
                                                                   aTGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGt                                                                             >  1:775098/1‑148 (MQ=255)
                                                                    tGACAAATACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGtt                                                                            >  2:1544313/1‑148 (MQ=255)
                                                                        aaaTACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACg                                                                        <  1:995342/148‑1 (MQ=255)
                                                                         aaTACGGGTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGc                                                                       <  2:1459075/148‑1 (MQ=255)
                                                                              gggTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTcca                                                                  <  1:419778/148‑1 (MQ=255)
                                                                              gggTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTcca                                                                  >  2:3116789/1‑148 (MQ=255)
                                                                               ggTATTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTccac                                                                 <  1:1618492/148‑1 (MQ=255)
                                                                                  aTTCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAg                                                              >  1:1424092/1‑148 (MQ=255)
                                                                                   ttCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGt                                                             >  2:3184542/1‑148 (MQ=255)
                                                                                    tCGCGCCGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTg                                                            >  2:679063/1‑148 (MQ=255)
                                                                                         ccGATAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGaa                                                       <  2:276242/148‑1 (MQ=255)
                                                                                            aTAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAg                                                    >  1:122994/1‑148 (MQ=255)
                                                                                             tAACAACTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGATTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTACACCAGTGAGGAATAGc                                                   >  2:2796600/1‑148 (MQ=255)
                                                                                              aacaacTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCg                                                  >  2:2131003/1‑148 (MQ=255)
                                                                                                 aacTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCgggt                                               >  1:3508378/1‑148 (MQ=255)
                                                                                                   cTTCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCgggtgg                                             <  2:781100/148‑1 (MQ=255)
                                                                                                   cTTCACCGCCAGTTTCATCAATACTGACCTGCTCTGCAATCTGAATTTGTTTCTCTTTCACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCgggtgg                                             <  2:1849545/148‑1 (MQ=255)
                                                                                                    ttCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCgggtggg                                            >  1:36788/1‑148 (MQ=255)
                                                                                                    ttCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCgggtggg                                            >  1:2494146/1‑148 (MQ=255)
                                                                                                    ttCACCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCgggtggg                                            >  1:1408372/1‑148 (MQ=255)
                                                                                                       aCCTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCgggtggg                                            <  1:2925709/145‑1 (MQ=255)
                                                                                                        ccTGCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTg                                        <  2:1635058/148‑1 (MQ=255)
                                                                                                           gCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTcg                                     >  1:3386220/1‑148 (MQ=255)
                                                                                                           gCAGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTcg                                     >  1:375907/1‑148 (MQ=255)
                                                                                                             aGTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTcgcg                                   <  2:42069/148‑1 (MQ=255)
                                                                                                              gTTTCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCgcgc                                  >  1:2179633/1‑148 (MQ=255)
                                                                                                                ttCATCAATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTg                                >  2:2922912/1‑148 (MQ=255)
                                                                                                                    tcaATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGt                            >  2:3382580/1‑148 (MQ=255)
                                                                                                                    tcaATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGt                            <  2:2606019/148‑1 (MQ=255)
                                                                                                                     caATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGtt                           >  1:3218728/1‑148 (MQ=255)
                                                                                                                     caATAATTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGtt                           <  2:2882001/148‑1 (MQ=255)
                                                                                                                      aataatTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTc                          >  1:2741213/1‑148 (MQ=255)
                                                                                                                       ataatTACCTGCTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCTGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTcg                         <  1:1906630/148‑1 (MQ=255)
                                                                                                                        taatTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTcgc                        >  2:1840857/1‑148 (MQ=255)
                                                                                                                         aatTACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGACTAGCGGGTGGGTTTGTCGCGCTGAGGTTcgcg                       >  1:893212/1‑148 (MQ=255)
                                                                                                                           ttACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTACGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGGTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTc                     >  2:1275598/1‑148 (MQ=255)
                                                                                                                           ttACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTc                     <  2:3463642/148‑1 (MQ=255)
                                                                                                                            tACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTcc                    <  1:3184521/148‑1 (MQ=255)
                                                                                                                            tACCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTcc                    <  1:3285106/148‑1 (MQ=255)
                                                                                                                             aCCTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCa                   >  2:472029/1‑148 (MQ=255)
                                                                                                                              ccTGGTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAg                  >  2:2754440/1‑148 (MQ=255)
                                                                                                                              ccTGGTCTGCAATCTGAATTTTTTTCTATTTAACATTAAGCTTCAGATAGTATTTCAGATCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAg                  >  1:670818/1‑148 (MQ=255)
                                                                                                                                 ggTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTaa               >  2:1897874/1‑148 (MQ=255)
                                                                                                                                  gTCTGCAATCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTaaa              <  1:541422/148‑1 (MQ=255)
                                                                                                                                         aTCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCa       <  2:2440904/148‑1 (MQ=255)
                                                                                                                                         aTCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCa       >  1:139825/1‑148 (MQ=255)
                                                                                                                                         aTCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCa       <  2:368838/148‑1 (MQ=255)
                                                                                                                                          tCTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCaa      >  2:2094372/1‑148 (MQ=255)
                                                                                                                                           cTGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAAt     >  1:2626311/1‑148 (MQ=255)
                                                                                                                                            tGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATa    <  2:1423684/148‑1 (MQ=255)
                                                                                                                                            tGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATa    <  2:1424096/148‑1 (MQ=255)
                                                                                                                                            tGAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATa    <  1:472029/148‑1 (MQ=255)
                                                                                                                                             gAATTTTTTTCTCTTTAACATTAAGCTTCAGATAGTTTTTCAGTTCCTGTGACAGAAGGTTGGGAATAAGCGGTTGACGCTTCCACCAGTGCGGAATAGCGGGTGGGTTTGTCGCGCTGAGGTTCGCGTCCAGTAAAAAAGGCAATaa   >  1:1005801/1‑148 (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: