Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 524,496 T→G 100% pseudogene (401/491 nt) ybbD → pseudogene

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273524,4960TG98.6% 1225.3 / ‑0.9 305pseudogene (401/491 nt)ybbDpseudogene
Reads supporting (aligned to +/- strand):  ref base T (0/6);  new base G (144/155);  total (144/161)
Fisher's exact test for biased strand distribution p-value = 3.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 309 total aligned reads displayed.

CCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGAGGAAAA  >  CP009273/524351‑524643
                                                                                                                                                 |                                                                                                                                                   
  cTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTaa                                                                                                                                                 <  1:2144651/148‑1 (MQ=255)
     aCAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAAcc                                                                                                                                              >  1:898168/1‑148 (MQ=255)
     aCAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAAcc                                                                                                                                              >  2:442703/1‑148 (MQ=255)
       aaaTGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCtt                                                                                                                                            >  2:1915151/1‑148 (MQ=255)
        aaTGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCttt                                                                                                                                           <  2:367393/148‑1 (MQ=255)
         atgatgatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCtttt                                                                                                                                          >  2:302165/1‑148 (MQ=255)
         atgatgatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCtttt                                                                                                                                          >  2:1103080/1‑148 (MQ=255)
         atgatgatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCtttt                                                                                                                                          >  2:674078/1‑148 (MQ=255)
          tgatgatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCttttt                                                                                                                                         >  2:2557392/1‑148 (MQ=255)
           gatgatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCtttttt                                                                                                                                        <  2:3410557/148‑1 (MQ=255)
           gatgatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCtttttt                                                                                                                                        <  2:1513938/148‑1 (MQ=255)
             tgatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCtttttttt                                                                                                                                      >  1:1606982/1‑148 (MQ=255)
             tgatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCtttttttt                                                                                                                                      >  2:1752693/1‑148 (MQ=255)
              gatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAAcc                                                                                                                                              >  1:2823252/1‑139 (MQ=255)
              gatgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTttg                                                                                                                                     >  2:1294958/1‑148 (MQ=255)
               atgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTttgt                                                                                                                                    <  2:817388/148‑1 (MQ=255)
               atgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTttgt                                                                                                                                    <  2:3535147/148‑1 (MQ=255)
               atgGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTttgt                                                                                                                                    <  2:2734667/148‑1 (MQ=255)
                tgGGGGTAGTTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTttgtt                                                                                                                                   <  2:2381606/148‑1 (MQ=255)
                  ggggTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTgg                                                                                                                                 <  1:3116876/148‑1 (MQ=255)
                  ggggTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTgg                                                                                                                                 <  1:1208608/148‑1 (MQ=255)
                    ggTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTg                                                                                                                               >  2:2516779/1‑148 (MQ=255)
                       aGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGaaa                                                                                                                            <  2:426633/148‑1 (MQ=255)
                          tttAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTa                                                                                                                         <  2:2133903/148‑1 (MQ=255)
                           ttAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTaa                                                                                                                        <  2:3585556/148‑1 (MQ=255)
                               gAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGtt                                                                                                                    >  1:1404020/1‑148 (MQ=255)
                               gAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGtt                                                                                                                    >  1:1452001/1‑148 (MQ=255)
                                aaaaCATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTa                                                                                                                   <  2:313696/148‑1 (MQ=255)
                                aaaaCATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTa                                                                                                                   <  2:323339/148‑1 (MQ=255)
                                aaaaCATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTa                                                                                                                   <  2:2903390/148‑1 (MQ=255)
                                     aTTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCat                                                                                                              <  2:700801/148‑1 (MQ=255)
                                      ttGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCatt                                                                                                             <  1:1989922/148‑1 (MQ=255)
                                       tGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCatta                                                                                                            <  1:3279973/148‑1 (MQ=255)
                                             tGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTtat                                                                                                      >  2:2474249/1‑148 (MQ=255)
                                              gAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTtata                                                                                                     <  1:1495256/148‑1 (MQ=255)
                                                gggCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAAc                                                                                                   >  1:3636077/1‑148 (MQ=255)
                                                    cAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAat                                                                                               <  1:747441/148‑1 (MQ=255)
                                                     aGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAATGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAata                                                                                              <  2:1847292/148‑1 (MQ=255)
                                                      gggAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAatat                                                                                             >  1:345857/1‑148 (MQ=255)
                                                      gggAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAatat                                                                                             >  2:2837055/1‑148 (MQ=255)
                                                        gAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTg                                                                                           <  1:1060943/148‑1 (MQ=255)
                                                         aaCAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTgg                                                                                          <  1:875351/148‑1 (MQ=255)
                                                          aCAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGa                                                                                         <  2:570735/148‑1 (MQ=255)
                                                           cAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAg                                                                                        >  1:3350082/1‑148 (MQ=255)
                                                             aGTGGCGAAACTCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCt                                                                                      <  2:3266393/148‑1 (MQ=255)
                                                             aGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCt                                                                                      <  1:1360004/148‑1 (MQ=255)
                                                             aGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCt                                                                                      >  1:1887655/1‑148 (MQ=255)
                                                             aGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCt                                                                                      <  2:2888072/148‑1 (MQ=255)
                                                               tGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTgg                                                                                    <  2:1136288/148‑1 (MQ=255)
                                                               tGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTgg                                                                                    >  1:2311041/1‑148 (MQ=255)
                                                                 tcGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGAt                                                                                  <  2:2407802/147‑1 (MQ=255)
                                                                     aaaTCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGcc                                                                              >  1:376542/1‑148 (MQ=255)
                                                                     aaaTCGTAACAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACATTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGcc                                                                              <  2:1989314/148‑1 (MQ=255)
                                                                      aaTCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCt                                                                             >  1:596258/1‑148 (MQ=255)
                                                                       aTCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTg                                                                            <  1:1004778/148‑1 (MQ=255)
                                                                       aTCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTg                                                                            >  2:3762174/1‑148 (MQ=255)
                                                                         cGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGAt                                                                          >  1:73797/1‑148 (MQ=255)
                                                                           tATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGc                                                                        <  2:3120245/148‑1 (MQ=255)
                                                                            aTCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCt                                                                       >  1:3041868/1‑148 (MQ=255)
                                                                             tCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGcgt                                                                      >  2:2647893/1‑146 (MQ=255)
                                                                              cAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTtg                                                                     >  1:2965091/1‑148 (MQ=255)
                                                                              aaaaGAATGATCCAGATACAAAAAAAGTGTGGTTTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAAATAATTTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTtg                                                                     <  1:127588/147‑1 (MQ=255)
                                                                                aaGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTtgtg                                                                   <  2:1440146/148‑1 (MQ=255)
                                                                                 aGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGc                                                                  >  1:245996/1‑148 (MQ=255)
                                                                                   aaTGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCag                                                                >  1:3537213/1‑148 (MQ=255)
                                                                                   aaTGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCag                                                                >  1:435406/1‑148 (MQ=255)
                                                                                    aTGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCaga                                                               >  1:3806042/1‑148 (MQ=255)
                                                                                     tGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCagag                                                              >  1:421670/1‑148 (MQ=255)
                                                                                        tCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTaa                                                           <  2:1404020/148‑1 (MQ=255)
                                                                                         ccAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAAt                                                          >  2:1882797/1‑148 (MQ=255)
                                                                                            gATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGc                                                       <  1:1823484/148‑1 (MQ=255)
                                                                                             aTACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCt                                                      >  1:2685362/1‑148 (MQ=255)
                                                                                              tACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCtt                                                     >  1:1702030/1‑148 (MQ=255)
                                                                                              tACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCtt                                                     <  1:561227/148‑1 (MQ=255)
                                                                                                 aaaaaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGaa                                                  >  1:2042690/1‑148 (MQ=255)
                                                                                                    aaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTa                                               <  1:1750591/148‑1 (MQ=255)
                                                                                                    aaaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTa                                               >  1:3781554/1‑148 (MQ=255)
                                                                                                     aaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTaa                                              >  2:2278789/1‑148 (MQ=255)
                                                                                                     aaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTaa                                              >  1:1299311/1‑148 (MQ=255)
                                                                                                     aaaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTaa                                              >  2:1277881/1‑148 (MQ=255)
                                                                                                      aaGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATTTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAag                                             >  1:2400721/1‑148 (MQ=255)
                                                                                                       aGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAaga                                            <  1:3363219/148‑1 (MQ=255)
                                                                                                       aGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAaga                                            <  2:2436132/148‑1 (MQ=255)
                                                                                                         tgtgGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGt                                          >  1:1183799/1‑148 (MQ=255)
                                                                                                                       ggTAAGATACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGa                            >  2:3327283/1‑148 (MQ=255)
                                                                                                                              tACCAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAgtcgtc                     >  1:2087276/1‑148 (MQ=255)
                                                                                                                                 cAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAAc                  >  2:963556/1‑148 (MQ=255)
                                                                                                                                 cAGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAAc                  >  1:1848027/1‑148 (MQ=255)
                                                                                                                                  aGGGCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACt                 <  2:584751/148‑1 (MQ=255)
                                                                                                                                   gggCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACtt                <  2:1568139/148‑1 (MQ=255)
                                                                                                                                   gggCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACtt                <  2:3429380/148‑1 (MQ=255)
                                                                                                                                    ggCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTc               >  1:3719318/1‑148 (MQ=255)
                                                                                                                                    ggCAAGGGGATGGGGTAAACCTATTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAAATTc               <  2:39715/148‑1 (MQ=255)
                                                                                                                                     gCAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGTAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCg              <  2:3771930/148‑1 (MQ=255)
                                                                                                                                      cAAGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGt             >  2:1863321/1‑148 (MQ=255)
                                                                                                                                       aaGGGGATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTa            <  2:3349995/148‑1 (MQ=255)
                                                                                                                                         ggggATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGt          >  2:164127/1‑148 (MQ=255)
                                                                                                                                         ggggATGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGt          >  2:2649913/1‑148 (MQ=255)
                                                                                                                                             aTGGGGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGAgg      >  2:2201444/1‑148 (MQ=255)
                                                                                                                                                 tgtAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGAGGaaaa  >  2:162004/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CCCTAACAAATGATGATGGGGGTAGGTTTAAGAAAACATTGACAGTGAGGGCCAGGGAACAAGTGGCGAAAATCGTATCAAAGAATGATCCAGATACAAAAAAAGTGTGGTGTAAATATGGTAAGATACCAGGGCAAGGGGATGGTGTAAACCTTTTTTTTGTTGGTGAAATTAATGTTACGCATTATTTTATAACAAATATTGGAGCTGGATTGCCTGATGCTTGTGCAGAGTAATTGCTTGAATTAAGAGTCTATCCCATATCGAAGTCGTCAACTTCGTAGTGAGGAAAA  >  CP009273/524351‑524643

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