Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,626,589 T→A 8.1% pseudogene (37/414 nt) yhiL ← pseudogene, DUF4049 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,626,5890TA 8.1% 792.1 / 49.0 258pseudogene (37/414 nt)yhiLpseudogene, DUF4049 family protein
Reads supporting (aligned to +/- strand):  ref base T (136/101);  new base A (8/13);  total (144/114)
Fisher's exact test for biased strand distribution p-value = 1.09e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

Only 100 of 267 total aligned reads displayed.

TGTTATCAAAATAATTAGCATATACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTTATACGTATCACCGGCAGAGAGT  >  CP009273/3626442‑3626736
                                                                                                                                                   |                                                                                                                                                   
   tATCAAAATAATTAGCATATACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATaaaa                                                                                                                                                  <  1:1957774/148‑1 (MQ=255)
   tATCAAAATAATTAGCATATACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATaaaa                                                                                                                                                  >  2:3289065/1‑148 (MQ=255)
   tATCAAAATAATTAGCATATACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTaaa                                                                                                                                                  <  2:1372629/148‑1 (MQ=255)
     tCAAAATAATTAGCATATACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATg                                                                                                                                                >  1:3048517/1‑148 (MQ=255)
       aaaaTAATTAGCATATACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATAAAATGGa                                                                                                                                              <  1:1731302/148‑1 (MQ=255)
         aataatTAGCATATACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAAc                                                                                                                                            >  1:2176518/1‑148 (MQ=255)
         aataatTAGCATATACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGAATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAAc                                                                                                                                            >  1:1192097/1‑148 (MQ=255)
           taatTAGCATATACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCt                                                                                                                                          >  2:637086/1‑148 (MQ=255)
            aatTAGCATATACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTg                                                                                                                                         >  2:321491/1‑148 (MQ=255)
             atTAGCATATACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGc                                                                                                                                        >  1:598398/1‑148 (MQ=255)
             atTAGCATATACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTATTCATCTGATTAAATGGAACCTGc                                                                                                                                        >  2:517345/1‑148 (MQ=255)
              ttAGCATATACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCa                                                                                                                                       >  2:3153751/1‑148 (MQ=255)
              ttAGCATATACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCa                                                                                                                                       >  2:3151628/1‑148 (MQ=255)
              ttAGCATATACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATAAAATGGAACCTGCa                                                                                                                                       <  1:2267679/148‑1 (MQ=255)
               tAGCATATACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCaa                                                                                                                                      >  1:1652997/1‑148 (MQ=255)
                  cATATACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGc                                                                                                                                   >  1:2490258/1‑148 (MQ=255)
                   atatACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCa                                                                                                                                  <  1:386717/148‑1 (MQ=255)
                   atatACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCa                                                                                                                                  >  2:722391/1‑148 (MQ=255)
                    tataCAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCAATGTGTTCTTTATGTTTTTTATTAAAAATATTGTCTAATTCCTCAGGATCAGTAGGATTTTTCATCTTATTAAATGGAACATGCAAAGCAc                                                                                                                                 >  1:937049/1‑148 (MQ=255)
                     ataCAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCAcc                                                                                                                                >  2:3106769/1‑148 (MQ=255)
                      taCAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCg                                                                                                                               >  2:617632/1‑148 (MQ=255)
                        cAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCgag                                                                                                                             <  1:1034470/148‑1 (MQ=255)
                         aGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCgaga                                                                                                                            <  1:2329979/148‑1 (MQ=255)
                           gtgtGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAGAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACa                                                                                                                          <  2:2682735/148‑1 (MQ=255)
                           gtgtGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACa                                                                                                                          <  1:1485548/148‑1 (MQ=255)
                           gtgtGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACa                                                                                                                          <  2:2020044/148‑1 (MQ=255)
                                   ttCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAAGAtt                                                                                                                  >  2:901364/1‑148 (MQ=255)
                                   ttCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACAtt                                                                                                                  <  1:2759203/148‑1 (MQ=255)
                                   ttCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACAtt                                                                                                                  >  2:903268/1‑148 (MQ=255)
                                    tCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACAttt                                                                                                                 <  2:1643507/148‑1 (MQ=255)
                                      aaaGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAAACTGCAAAGCACCGAGACAATAACAttttt                                                                                                               >  2:2542381/1‑148 (MQ=255)
                                          tATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTttct                                                                                                           >  1:889344/1‑148 (MQ=255)
                                            ttGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTttcttc                                                                                                         >  1:1349479/1‑148 (MQ=255)
                                            ttGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTttcttc                                                                                                         >  1:407811/1‑148 (MQ=255)
                                            ttGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTttcttc                                                                                                         <  2:29059/148‑1 (MQ=255)
                                             tGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCa                                                                                                        >  1:1979156/1‑148 (MQ=255)
                                             tGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCa                                                                                                        >  1:1114763/1‑148 (MQ=255)
                                               cTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCAtc                                                                                                      <  2:1949293/148‑1 (MQ=255)
                                                ttCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCAtct                                                                                                     <  2:854161/148‑1 (MQ=255)
                                                 tCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCAtctc                                                                                                    <  1:56055/148‑1 (MQ=255)
                                                    gATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCtta                                                                                                 >  1:616994/1‑148 (MQ=255)
                                                       tAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCttatta                                                                                              <  2:1218432/148‑1 (MQ=255)
                                                       tAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATAAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTAATCTCttctta                                                                                              >  2:2238253/1‑148 (MQ=255)
                                                          aTGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTc                                                                                           >  1:3167839/1‑148 (MQ=255)
                                                            ggAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCg                                                                                         <  1:364738/148‑1 (MQ=255)
                                                                     gATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACACTAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTg                                                                                >  2:2299224/1‑148 (MQ=255)
                                                                      atcatcCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTgg                                                                               >  2:3081177/1‑148 (MQ=255)
                                                                      atcatcCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTAATCTCTTATTATTCCGTGATGACTgg                                                                               >  2:3081261/1‑148 (MQ=255)
                                                                          tcCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTaa                                                                           <  2:926132/148‑1 (MQ=255)
                                                                             aTGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAAc                                                                        >  2:1553866/1‑148 (MQ=255)
                                                                              tgtgTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACt                                                                       <  1:1624953/148‑1 (MQ=255)
                                                                                 gtTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGc                                                                    >  1:45534/1‑148 (MQ=255)
                                                                                   tCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCTAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCt                                                                  <  2:351081/148‑1 (MQ=255)
                                                                                   tCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCt                                                                  <  2:1098676/148‑1 (MQ=255)
                                                                                    cTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCtt                                                                 >  2:3131481/1‑148 (MQ=255)
                                                                                       tATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCt                                                              >  2:1308568/1‑148 (MQ=255)
                                                                                       tATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCt                                                              >  2:2194051/1‑148 (MQ=255)
                                                                                        aTGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTg                                                             <  1:2904030/148‑1 (MQ=255)
                                                                                        aTGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTg                                                             >  2:3310141/1‑148 (MQ=255)
                                                                                          gTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATCCCGTGATGACTGGTTAAAACTTGCCTTTCTGAg                                                           <  2:495628/148‑1 (MQ=255)
                                                                                            tttttATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTc                                                         >  2:1911376/1‑148 (MQ=255)
                                                                                            tttttATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCAGTGATTACTGGTTAAAACTTGCATTAATGAGTc                                                         >  1:3289918/1‑148 (MQ=255)
                                                                                                  ttaAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAAtt                                                   <  2:2490258/148‑1 (MQ=255)
                                                                                                   taAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTa                                                  >  2:714144/1‑148 (MQ=255)
                                                                                                      aaaTATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACaa                                               <  2:1988412/148‑1 (MQ=255)
                                                                                                      aaaTATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACaa                                               <  2:1239184/148‑1 (MQ=255)
                                                                                                        atatTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAAc                                             >  2:2376725/1‑148 (MQ=255)
                                                                                                         tatTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACa                                            >  1:2561266/1‑148 (MQ=255)
                                                                                                          atTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACAt                                           <  1:722391/148‑1 (MQ=255)
                                                                                                                tAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCg                                     <  1:3106769/148‑1 (MQ=255)
                                                                                                                tAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCg                                     <  2:2148105/148‑1 (MQ=255)
                                                                                                                 aaTTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGt                                    <  1:1302589/148‑1 (MQ=255)
                                                                                                                  aTTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGtt                                   >  1:2414594/1‑148 (MQ=255)
                                                                                                                   ttCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGttt                                  >  2:955136/1‑148 (MQ=255)
                                                                                                                     ccTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTtta                                <  1:2788757/148‑1 (MQ=255)
                                                                                                                     ccTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTATTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTtta                                <  2:427201/148‑1 (MQ=255)
                                                                                                                       tCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTttaat                              >  2:2258849/1‑148 (MQ=255)
                                                                                                                          ggATCCGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTttaattaa                           >  2:1538611/1‑148 (MQ=255)
                                                                                                                          ggATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTttaattaa                           <  2:1759133/148‑1 (MQ=255)
                                                                                                                          ggATCAGTAGGATTTTTCATCTGATAAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTttaattaa                           >  1:83963/1‑148 (MQ=255)
                                                                                                                           gATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAg                          <  2:296561/148‑1 (MQ=255)
                                                                                                                           gATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAg                          >  1:1918372/1‑148 (MQ=255)
                                                                                                                           gATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAg                          <  2:1803746/148‑1 (MQ=255)
                                                                                                                               aGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTtat                      <  1:1776568/148‑1 (MQ=255)
                                                                                                                               aGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTtat                      <  1:3144067/148‑1 (MQ=255)
                                                                                                                                     aTTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTTATACGTAt                <  2:2244287/148‑1 (MQ=255)
                                                                                                                                           cATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTTATACGTATCACCgg          >  1:3301388/1‑148 (MQ=255)
                                                                                                                                           cATCTGATAAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTTATACGTATCACCgg          <  1:3193380/148‑1 (MQ=255)
                                                                                                                                            aTCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTTATACGTATCACCGGc         <  2:3003304/148‑1 (MQ=255)
                                                                                                                                             tCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTTATACGTATCACCGGCa        >  2:664787/1‑148 (MQ=255)
                                                                                                                                              cTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTTATACGTATCACCGGCag       <  1:3151676/148‑1 (MQ=255)
                                                                                                                                              cTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTTATACGTATCACCGGCag       <  2:889344/148‑1 (MQ=255)
                                                                                                                                              cTGATAAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTTATACGTATCACCGGCag       >  2:3246035/1‑148 (MQ=255)
                                                                                                                                              cTGATAAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTTATACGTATCACCGGCag       <  1:2314562/148‑1 (MQ=255)
                                                                                                                                               tGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTTATACGTATCACCGGCaga      <  1:3020410/148‑1 (MQ=255)
                                                                                                                                               tGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTTATACGTATCACCGGCaga      <  2:3048517/148‑1 (MQ=255)
                                                                                                                                                 aTTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTTATACGTATCACCGGCagaga    >  1:1017956/1‑148 (MQ=255)
                                                                                                                                                   tAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTTATACGTATCACCGGCAGAGAGt  <  1:321491/148‑1 (MQ=255)
                                                                                                                                                   tAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTTATACGTATCACCGGCAGAGAGt  <  2:1584336/148‑1 (MQ=255)
                                                                                                                                                   tAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTTATACGTATCACCGGCAGAGAGt  >  2:1858411/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGTTATCAAAATAATTAGCATATACAGGTGTGGGTTTCAAAGTATTGCTTCGGATTAGATGGAATAGAGGATCATCCATGTGTTCTTTATGTTTTTTATTAAAAATATTGGCTAATTCCTCAGGATCAGTAGGATTTTTCATCTGATTAAATGGAACCTGCAAAGCACCGAGACAATAACATTTTTTTCTTCATCTCTTATTATTCCGTGATGACTGGTTAAAACTTGCCTTTCTGAGTCATAATTACAAACATCGAGCGTTTTAATTAAGTTATACGTATCACCGGCAGAGAGT  >  CP009273/3626442‑3626736

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: