Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 576,998 G→T 100% pseudogene (189/309 nt) aaaD → pseudogene; DLP12 prophage; tail fiber assembly protein family;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273576,9980GT98.4% 1029.5 / 4.3 311pseudogene (189/309 nt)aaaDpseudogene; DLP12 prophage; tail fiber assembly protein family;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base G (2/3);  new base T (149/157);  total (151/160)
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.92e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 311 total aligned reads displayed.

GTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTC  >  CP009273/576852‑577142
                                                                                                                                                  |                                                                                                                                                
gTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTAt                                                                                                                                                >  2:2649993/1‑149 (MQ=255)
 tCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATg                                                                                                                                               >  2:2247981/1‑149 (MQ=255)
  ccGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATgg                                                                                                                                              >  2:1565422/1‑149 (MQ=255)
     cgcCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCa                                                                                                                                           <  2:4513860/149‑1 (MQ=255)
         gggCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGa                                                                                                                                       >  1:220051/1‑149 (MQ=255)
           gCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACt                                                                                                                                     <  1:721213/149‑1 (MQ=255)
            cTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTa                                                                                                                                    <  2:2757038/149‑1 (MQ=255)
                  tCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCg                                                                                                                              >  1:2335747/1‑149 (MQ=255)
                   cACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGt                                                                                                                             <  2:4194974/149‑1 (MQ=255)
                   cACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGt                                                                                                                             <  2:1804768/149‑1 (MQ=255)
                       gTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGACAGAATCCGCATACCCGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCAACATGAATGCTATGGGCAGCGACTACATCCGTGAgg                                                                                                                         <  1:1646652/149‑1 (MQ=255)
                       gTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAgg                                                                                                                         <  1:2582582/149‑1 (MQ=255)
                          cAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGa                                                                                                                      <  2:140003/149‑1 (MQ=255)
                           aGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGaa                                                                                                                     <  2:87062/149‑1 (MQ=255)
                            ggCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAAt                                                                                                                    >  2:4155512/1‑149 (MQ=255)
                             gCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAAtg                                                                                                                   <  1:2174818/149‑1 (MQ=255)
                              ccGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAAtgt                                                                                                                  <  2:1050417/149‑1 (MQ=255)
                              ccGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAAtgt                                                                                                                  <  2:1287686/149‑1 (MQ=255)
                                     cacaGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAg                                                                                                           <  2:2170346/149‑1 (MQ=255)
                                       caGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACCTCCGTGAGGTGACTGTGGTGACGTc                                                                                                         >  1:1843696/1‑149 (MQ=255)
                                       caGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTc                                                                                                         >  2:2564784/1‑149 (MQ=255)
                                        aGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCt                                                                                                        <  2:4040127/149‑1 (MQ=255)
                                        aGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCt                                                                                                        <  2:3900405/149‑1 (MQ=255)
                                          aCCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGc                                                                                                      <  2:3098981/149‑1 (MQ=255)
                                                gTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgtgt                                                                                                <  1:4043215/149‑1 (MQ=255)
                                                gTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCgtgt                                                                                                >  2:490822/1‑149 (MQ=255)
                                                  tGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCg                                                                                              >  1:1315229/1‑149 (MQ=255)
                                                            gATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaa                                                                                    >  2:4130036/1‑149 (MQ=255)
                                                             aTGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaaa                                                                                   <  2:2372094/149‑1 (MQ=255)
                                                             aTGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaaa                                                                                   >  1:2513765/1‑149 (MQ=255)
                                                              tGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaaaa                                                                                  <  1:2632887/149‑1 (MQ=255)
                                                              tGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCaaaa                                                                                  <  1:1736328/149‑1 (MQ=255)
                                                               gCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAt                                                                                 <  2:2113778/149‑1 (MQ=255)
                                                                cTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtg                                                                                <  1:4694314/149‑1 (MQ=255)
                                                                 tAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgc                                                                               <  1:3260712/149‑1 (MQ=255)
                                                                   aTTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgctg                                                                             >  1:1033892/1‑149 (MQ=255)
                                                                    ttACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgctgc                                                                            >  2:1511693/1‑149 (MQ=255)
                                                                    ttACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAAtgctgc                                                                            >  2:3168283/1‑149 (MQ=255)
                                                                     tACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATgctgct                                                                           <  2:2066660/149‑1 (MQ=255)
                                                                      aCTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGctgctg                                                                          >  1:3887660/1‑149 (MQ=255)
                                                                      aCTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGctgctg                                                                          >  1:3095555/1‑149 (MQ=255)
                                                                        tATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGGGAGGTGAAGGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTggg                                                                        >  2:2352983/1‑149 (MQ=255)
                                                                           ctcCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGgtgt                                                                     >  1:1054766/1‑149 (MQ=255)
                                                                             cccGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGttt                                                                   <  1:1972537/149‑1 (MQ=255)
                                                                             cccGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGttt                                                                   <  1:4783651/149‑1 (MQ=255)
                                                                              ccGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTa                                                                  >  1:1420202/1‑149 (MQ=255)
                                                                               cGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTAt                                                                 <  1:1704258/149‑1 (MQ=255)
                                                                                 aaaGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGc                                                               <  1:4219199/149‑1 (MQ=255)
                                                                                  aaGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGACTGTGGTGAAGTCTGCCCGTGTCGGTGATTCCAAAATGCTGCTGGGTGTTTATGcc                                                              >  2:1943986/1‑149 (MQ=255)
                                                                                  aaGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGcc                                                              <  1:459656/149‑1 (MQ=255)
                                                                                     aaTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTAc                                                           <  2:4462622/149‑1 (MQ=255)
                                                                                     aaTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTAc                                                           >  2:4244766/1‑149 (MQ=255)
                                                                                     aaTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTAc                                                           >  1:244263/1‑149 (MQ=255)
                                                                                     aaTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTAc                                                           >  1:4284052/1‑149 (MQ=255)
                                                                                      aTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACt                                                          >  1:281597/1‑149 (MQ=255)
                                                                                      aTCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACt                                                          >  1:4021095/1‑149 (MQ=255)
                                                                                       tCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACtt                                                         <  1:1961736/149‑1 (MQ=255)
                                                                                       tCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACtt                                                         <  1:3213561/149‑1 (MQ=255)
                                                                                        ccGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACttt                                                        <  2:1900171/149‑1 (MQ=255)
                                                                                         cGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTta                                                       >  1:3575047/1‑149 (MQ=255)
                                                                                         cGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTta                                                       >  2:2382306/1‑149 (MQ=255)
                                                                                          gCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTtat                                                      <  2:594096/149‑1 (MQ=255)
                                                                                           cATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTtata                                                     >  1:798565/1‑149 (MQ=255)
                                                                                           cATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTtata                                                     >  1:3190144/1‑149 (MQ=255)
                                                                                              aCCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATagag                                                  <  1:2530627/149‑1 (MQ=255)
                                                                                                 aGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCAt                                               >  2:1022928/1‑149 (MQ=255)
                                                                                                    aaGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAag                                            >  1:4135213/1‑149 (MQ=255)
                                                                                                     aGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagc                                           <  1:4154912/149‑1 (MQ=255)
                                                                                                       ggCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagcag                                         >  2:2322942/1‑149 (MQ=255)
                                                                                                       ggCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagcag                                         <  1:4082067/149‑1 (MQ=255)
                                                                                                        gcgcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAagcagc                                        <  1:3330015/149‑1 (MQ=255)
                                                                                                          gcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATCAGCAgcgc                                      >  2:4287498/1‑149 (MQ=255)
                                                                                                          gcTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAgcgc                                      >  1:1555270/1‑149 (MQ=255)
                                                                                                            tGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCaa                                    >  1:4343952/1‑149 (MQ=255)
                                                                                                                  acacTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCt                              <  1:3590739/149‑1 (MQ=255)
                                                                                                                  acacTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCt                              <  2:3610367/149‑1 (MQ=255)
                                                                                                                  acacTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCt                              >  2:3509952/1‑149 (MQ=255)
                                                                                                                   cacTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCtt                             <  1:623101/149‑1 (MQ=255)
                                                                                                                   cacTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCtt                             >  1:2303901/1‑149 (MQ=255)
                                                                                                                    acTGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTa                            >  1:2928913/1‑149 (MQ=255)
                                                                                                                      tGCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTAtt                          >  2:2871310/1‑149 (MQ=255)
                                                                                                                       gCCCTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTc                         >  2:3706752/1‑149 (MQ=255)
                                                                                                                        cccTTTCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTcc                        <  1:1436215/149‑1 (MQ=255)
                                                                                                                           tttCAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGc                     <  1:950983/149‑1 (MQ=255)
                                                                                                                              cAGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTgg                  <  2:3240566/149‑1 (MQ=255)
                                                                                                                               aGCGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGc                 >  2:4754968/1‑149 (MQ=255)
                                                                                                                                 cGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCtt               <  2:3304425/149‑1 (MQ=255)
                                                                                                                                 cGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCtt               <  1:376978/149‑1 (MQ=255)
                                                                                                                                 cGGGCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCtt               <  1:3361353/149‑1 (MQ=255)
                                                                                                                                  gggCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAATCAGCGCAACACCCTTATTCCAGCTGGCTTc              >  2:4682096/1‑149 (MQ=255)
                                                                                                                                    gCCATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGt            >  1:1805216/1‑149 (MQ=255)
                                                                                                                                     ccATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTg           >  2:3016468/1‑149 (MQ=255)
                                                                                                                                     ccATCATGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTg           <  2:3429294/149‑1 (MQ=255)
                                                                                                                                      catcatGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTgg          >  1:1920412/1‑149 (MQ=255)
                                                                                                                                      catcatGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTgg          <  2:1744217/149‑1 (MQ=255)
                                                                                                                                         catGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTg       >  1:2731521/1‑149 (MQ=255)
                                                                                                                                           tGAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGtt     >  2:3933818/1‑149 (MQ=255)
                                                                                                                                            gAATGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGttt    <  2:1708963/149‑1 (MQ=255)
                                                                                                                                              aTGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTc  >  1:1061722/1‑149 (MQ=255)
                                                                                                                                              aTGCTATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTc  <  1:3955152/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GTCCGCGCCGGGCTTCGCTCACTGTTCAGGCCGGAGCCACAGACCGCCGTTGAATGGGCGGATGCTAATTACTATCTCCCGAAAGAATCCGCATACCAGGAAGGGCGCTGGGAAACACTGCCCTTTCAGCGGGCCATCATGAATGCGATGGGCAGCGACTACATCCGTGAGGTGAATGTGGTGAAGTCTGCCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTC  >  CP009273/576852‑577142

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: