Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 523,413 G→T 31.9% pseudogene (6/708 nt) rhsH → pseudogene, Rhs familly protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273523,4130GT31.9% 451.8 / 235.5 257pseudogene (6/708 nt)rhsHpseudogene, Rhs familly protein
Reads supporting (aligned to +/- strand):  ref base G (116/59);  new base T (48/34);  total (164/93)
Fisher's exact test for biased strand distribution p-value = 2.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.13e-01

Only 100 of 261 total aligned reads displayed.

AAACCATATTTAATTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCG  >  CP009273/523265‑523562
                                                                                                                                                    |                                                                                                                                                     
 aaCCATATTTAATTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCt                                                                                                                                                     >  2:1518695/1‑150 (MQ=255)
    cATATTTAATTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTgat                                                                                                                                                  >  2:2567842/1‑150 (MQ=255)
      tatTTAATTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTgatgg                                                                                                                                                >  1:1128410/1‑150 (MQ=255)
          tAATTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGc                                                                                                                                            >  2:4763141/1‑150 (MQ=255)
            aTTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCgg                                                                                                                                          >  1:3020173/1‑150 (MQ=255)
            aTTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCgg                                                                                                                                          >  2:516381/1‑150 (MQ=255)
                   aaCTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGaaa                                                                                                                                   >  2:3884493/1‑150 (MQ=255)
                     cTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACa                                                                                                                                 >  1:446562/1‑150 (MQ=255)
                          aGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCg                                                                                                                            >  2:1199682/1‑150 (MQ=255)
                            aTAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTg                                                                                                                          >  1:1777724/1‑150 (MQ=255)
                            aTAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTg                                                                                                                          >  2:2755535/1‑150 (MQ=255)
                             tAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTgg                                                                                                                         >  1:2868361/1‑150 (MQ=255)
                              aGTTGTGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGa                                                                                                                        <  1:4337705/150‑1 (MQ=255)
                              aGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGa                                                                                                                        >  1:1979296/1‑150 (MQ=255)
                              aGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGa                                                                                                                        <  1:4302408/150‑1 (MQ=255)
                              aGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGa                                                                                                                        <  2:2910572/150‑1 (MQ=255)
                                   ggTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCggg                                                                                                                   >  1:921991/1‑150 (MQ=255)
                                     ttgttgAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGa                                                                                                                 <  2:1603357/150‑1 (MQ=255)
                                        ttgAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTa                                                                                                              <  2:3038520/150‑1 (MQ=255)
                                          gAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTatg                                                                                                            <  2:4498975/150‑1 (MQ=255)
                                               aGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAg                                                                                                       >  1:4144039/1‑150 (MQ=255)
                                               aGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAg                                                                                                       >  2:2287089/1‑150 (MQ=255)
                                                 cAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATGAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTg                                                                                                     >  2:1067276/1‑150 (MQ=255)
                                                    aGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTgggg                                                                                                  >  1:2444619/1‑150 (MQ=255)
                                                    aGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTgggg                                                                                                  >  2:3395348/1‑150 (MQ=255)
                                                    aGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTgggg                                                                                                  >  2:2014787/1‑150 (MQ=255)
                                                      ttACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCa                                                                                                <  2:3654258/150‑1 (MQ=255)
                                                      ttACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCa                                                                                                <  2:4145521/150‑1 (MQ=255)
                                                       tACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCaa                                                                                               >  1:2954219/1‑150 (MQ=255)
                                                          cTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCa                                                                                            <  1:2379955/150‑1 (MQ=255)
                                                           ttATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAg                                                                                           >  2:1443034/1‑150 (MQ=255)
                                                            tATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGc                                                                                          >  2:1424035/1‑150 (MQ=255)
                                                             aTGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCt                                                                                         <  2:4694480/150‑1 (MQ=255)
                                                             aTGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCt                                                                                         <  1:4508226/150‑1 (MQ=255)
                                                                aaGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCtgaa                                                                                      <  1:2036998/150‑1 (MQ=255)
                                                                 aGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCtgaat                                                                                     >  1:4358868/1‑150 (MQ=255)
                                                                  gTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCtgaatg                                                                                    >  1:4426467/1‑150 (MQ=255)
                                                                    tGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCtgaatgaa                                                                                  <  1:3997718/150‑1 (MQ=255)
                                                                      cGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAaga                                                                                >  2:4328329/1‑150 (MQ=255)
                                                                        ggatgtttatttattATATAAATTAATAATAAAAATGGATGTGTTTTGAATTTACTACATAGTAAATAATGCTGTCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAgaga                                                                              <  2:4526149/143‑1 (MQ=255)
                                                                         gtggtgTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCACCCAGCTGAATGAAGAGaa                                                                             >  2:2888662/1‑150 (MQ=255)
                                                                         gtggtgTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGaa                                                                             >  2:508977/1‑150 (MQ=255)
                                                                           ggtgTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAAcc                                                                           <  1:4696302/150‑1 (MQ=255)
                                                                           ggtgTAGTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAAcc                                                                           <  1:2723858/150‑1 (MQ=255)
                                                                             tgtATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCg                                                                         >  1:1127403/1‑150 (MQ=255)
                                                                             tgtATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCg                                                                         <  2:3717058/150‑1 (MQ=255)
                                                                                 tttattATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATc                                                                     >  1:1280817/1‑150 (MQ=255)
                                                                                   tattatAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCAc                                                                   >  2:3274543/1‑150 (MQ=255)
                                                                                   tattatAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCAc                                                                   >  1:3881351/1‑150 (MQ=255)
                                                                                    attatAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACAGTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCAcc                                                                  <  2:736006/150‑1 (MQ=255)
                                                                                      tataGAAATTAATAAGAAAAATTGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGCTGGTAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTg                                                                <  1:1565090/150‑1 (MQ=255)
                                                                                      tataGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTg                                                                <  1:2614611/150‑1 (MQ=255)
                                                                                       ataGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGc                                                               <  2:3737620/150‑1 (MQ=255)
                                                                                       ataGAAATTAATAAGAAAAATGGATGGGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGc                                                               >  2:3120403/1‑150 (MQ=255)
                                                                                           aaaTTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCa                                                           >  1:2453338/1‑150 (MQ=255)
                                                                                             aTTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATCATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGTAGAGAACCCGCATCACCTGCACCAGc                                                         >  2:4661524/1‑150 (MQ=255)
                                                                                             aTTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGc                                                         >  1:775290/1‑150 (MQ=255)
                                                                                              ttAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGcc                                                        <  1:3821154/150‑1 (MQ=255)
                                                                                                  taaGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTAc                                                    <  2:1136423/150‑1 (MQ=255)
                                                                                                    aGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCg                                                  >  2:3601114/1‑150 (MQ=255)
                                                                                                    aGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCg                                                  >  2:4761166/1‑150 (MQ=255)
                                                                                                      aaaaaTGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGc                                                >  1:3800651/1‑150 (MQ=255)
                                                                                                        aaaTGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGCGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTg                                              >  1:4263506/1‑150 (MQ=255)
                                                                                                          aTGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGcc                                            >  1:4557507/1‑150 (MQ=255)
                                                                                                          aTGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGcc                                            >  2:300659/1‑150 (MQ=255)
                                                                                                          aTGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGcc                                            >  2:2525595/1‑150 (MQ=255)
                                                                                                             gATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCggg                                         <  1:1035606/150‑1 (MQ=255)
                                                                                                             gATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCggg                                         <  1:4171213/150‑1 (MQ=255)
                                                                                                              aTGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCgggg                                        <  2:407125/150‑1 (MQ=255)
                                                                                                               tgtgTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGgc                                       >  2:2460561/1‑150 (MQ=255)
                                                                                                                 tgtTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGgcag                                     >  1:773711/1‑150 (MQ=255)
                                                                                                                 tgtTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGgcag                                     >  1:1607308/1‑150 (MQ=255)
                                                                                                                  gtTTTGAATTTCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGgcagc                                    >  2:788003/1‑150 (MQ=255)
                                                                                                                  gtTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGgcagc                                    >  2:1891228/1‑150 (MQ=255)
                                                                                                                   ttttGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGgccgca                                   >  1:4779690/1‑150 (MQ=255)
                                                                                                                   ttttGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGgcagca                                   <  1:3666651/150‑1 (MQ=255)
                                                                                                                     ttGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGACTGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGt                                 <  2:316187/150‑1 (MQ=255)
                                                                                                                     ttGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACGCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGt                                 <  1:504486/150‑1 (MQ=255)
                                                                                                                     ttGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGt                                 <  2:1715036/150‑1 (MQ=255)
                                                                                                                          tttCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCCGTATGAt                            >  2:4418734/1‑150 (MQ=255)
                                                                                                                          tttCCTACATAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGAt                            >  2:4418757/1‑150 (MQ=255)
                                                                                                                               tACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGa                       >  1:1413180/1‑150 (MQ=255)
                                                                                                                                 cATAGTAAATAATGCTGGTGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGt                     <  1:2135598/150‑1 (MQ=255)
                                                                                                                                 cATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGt                     <  2:3561753/150‑1 (MQ=255)
                                                                                                                                 cATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGt                     <  2:2957044/150‑1 (MQ=255)
                                                                                                                                  aTAGTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTc                    >  2:2509808/1‑150 (MQ=255)
                                                                                                                                   tAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCg                   >  1:3730398/1‑150 (MQ=255)
                                                                                                                                   tAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCg                   <  2:1046114/150‑1 (MQ=255)
                                                                                                                                     gTAAATAATGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCggg                 >  2:886097/1‑150 (MQ=255)
                                                                                                                                     gTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCggg                 <  2:3500895/150‑1 (MQ=255)
                                                                                                                                      tAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAAGGAAGAGACCCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCGGCAGTATGATAAGGAGTCgggg                <  2:1254368/150‑1 (MQ=255)
                                                                                                                                        aataatGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTAGGGGCt              <  1:3467643/150‑1 (MQ=255)
                                                                                                                                        aataatGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTAGGGGCt              <  1:1061055/150‑1 (MQ=255)
                                                                                                                                        aataatGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCt              <  1:2373396/150‑1 (MQ=255)
                                                                                                                                         ataatGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTg             >  2:2772949/1‑150 (MQ=255)
                                                                                                                                          taatGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGt            >  2:3928477/1‑150 (MQ=255)
                                                                                                                                          taatGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGt            >  1:3363716/1‑150 (MQ=255)
                                                                                                                                            atGCTGGCTCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGtac          <  2:3839725/150‑1 (MQ=255)
                                                                                                                                             tGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGtact         <  1:1510852/150‑1 (MQ=255)
                                                                                                                                                   cgcTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCACATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACaacc   <  2:2310941/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AAACCATATTTAATTACGGAACTACCAGATAGTTGGGTTGTTGAGGGAGCAAAGTTACCTTATGAAGTTGCGGGTGGTGTATTTATTATAGAAATTAATAAGAAAAATGGATGTGTTTTGAATTTCCTACATAGTAAATAATGCTGGCGCTGATGGATGCGGATGGAAACATTGCGTGGAGCGGGGAGTATGATGAGTGGGGCAACCAGCTGAATGAAGAGAACCCGCATCACCTGCACCAGCCGTACCGGCTGCCGGGGCAGCAGTATGATAAGGAGTCGGGGCTGTACTACAACCG  >  CP009273/523265‑523562

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: