Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 524,983 G→T 12.2% pseudogene (28/93 nt) ylbI → pseudogene, Rhs family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273524,9830GT12.2% 936.0 / 67.7 287pseudogene (28/93 nt)ylbIpseudogene, Rhs family protein
Reads supporting (aligned to +/- strand):  ref base G (122/130);  new base T (18/17);  total (140/147)
Fisher's exact test for biased strand distribution p-value = 8.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.70e-01

Only 100 of 289 total aligned reads displayed.

ATATTCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTGAGATTCTCAAGGT  >  CP009273/524836‑525129
                                                                                                                                                   |                                                                                                                                                  
 tatTCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGa                                                                                                                                                   >  1:1156984/1‑148 (MQ=255)
   ttCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACa                                                                                                                                                 <  2:1609458/148‑1 (MQ=255)
   ttCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACa                                                                                                                                                 >  2:29169/1‑148 (MQ=255)
    tCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACAt                                                                                                                                                <  1:2790748/148‑1 (MQ=255)
    tCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACAt                                                                                                                                                >  1:1722739/1‑148 (MQ=255)
     cAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATg                                                                                                                                               <  1:2054867/148‑1 (MQ=255)
        ccATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGccc                                                                                                                                            >  1:2357342/1‑148 (MQ=255)
        ccATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGccc                                                                                                                                            <  1:1375224/148‑1 (MQ=255)
             cATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTa                                                                                                                                       <  1:2182515/148‑1 (MQ=255)
             cATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTa                                                                                                                                       <  1:3024694/148‑1 (MQ=255)
             cATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTa                                                                                                                                       <  1:2463190/148‑1 (MQ=255)
                   ttaaAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTAcggcgg                                                                                                                                 >  1:2105820/1‑148 (MQ=255)
                    taaAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACggcggc                                                                                                                                >  1:1202724/1‑148 (MQ=255)
                         tGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGAt                                                                                                                           <  2:3029351/148‑1 (MQ=255)
                           ttACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAg                                                                                                                         >  1:2509991/1‑148 (MQ=255)
                           ttACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAg                                                                                                                         <  2:3508957/148‑1 (MQ=255)
                           ttACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAg                                                                                                                         <  2:1156470/148‑1 (MQ=255)
                           ttACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAg                                                                                                                         >  2:1679596/1‑148 (MQ=255)
                            tACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGt                                                                                                                        >  1:1836213/1‑148 (MQ=255)
                            tACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGt                                                                                                                        >  2:2363112/1‑148 (MQ=255)
                                 gtTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAgg                                                                                                                   <  2:2285738/148‑1 (MQ=255)
                                          gAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCgg                                                                                                          >  1:1921510/1‑148 (MQ=255)
                                               ccTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTa                                                                                                     >  1:385188/1‑148 (MQ=255)
                                               ccTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTa                                                                                                     <  2:3207543/148‑1 (MQ=255)
                                               ccTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTa                                                                                                     <  2:1425051/148‑1 (MQ=255)
                                               ccTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCATGCCGGGAAAGAAAGAGGAAATAAACCAATAGCGGAAAACCAGCGGCACGCAAGCGGGACAAGCCCCAGAACGGCGGCCCCATAGTCCCGCGATCAGCGGGGGTa                                                                                                     >  2:2828013/1‑148 (MQ=255)
                                                cTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTAc                                                                                                    <  1:489938/148‑1 (MQ=255)
                                                cTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTAc                                                                                                    >  1:782317/1‑148 (MQ=255)
                                                  aTCTGAAAATAAATCATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATACACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTg                                                                                                  <  1:3493207/148‑1 (MQ=255)
                                                  aTCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTg                                                                                                  >  2:814410/1‑148 (MQ=255)
                                                    cTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGAt                                                                                                >  2:2030619/1‑148 (MQ=255)
                                                    cTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGAt                                                                                                >  2:2578841/1‑148 (MQ=255)
                                                      gAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAg                                                                                              <  2:2812439/148‑1 (MQ=255)
                                                      gAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAg                                                                                              >  2:3053235/1‑148 (MQ=255)
                                                       aaaaTAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGt                                                                                             <  2:2509991/148‑1 (MQ=255)
                                                           taaataaTCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTcgcg                                                                                         <  1:151325/148‑1 (MQ=255)
                                                           taaataaTCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTcgcg                                                                                         >  2:1501866/1‑148 (MQ=255)
                                                           taaataaTCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTcgcg                                                                                         >  1:2908633/1‑148 (MQ=255)
                                                           taaataaTCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTcgcg                                                                                         >  2:1795731/1‑148 (MQ=255)
                                                            aaataaTCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTCCTGATAGTCgcgc                                                                                        >  1:2533698/1‑148 (MQ=255)
                                                               taatCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGa                                                                                     <  2:1866659/148‑1 (MQ=255)
                                                               taatCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGa                                                                                     >  1:1047309/1‑148 (MQ=255)
                                                               taatCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGa                                                                                     <  1:1580637/148‑1 (MQ=255)
                                                                aatCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGAc                                                                                    >  2:2901922/1‑148 (MQ=255)
                                                                  tCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTg                                                                                  <  1:3027886/148‑1 (MQ=255)
                                                                   ccTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGc                                                                                 <  1:346546/148‑1 (MQ=255)
                                                                   ccTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGc                                                                                 <  2:2075497/148‑1 (MQ=255)
                                                                   ccTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGc                                                                                 <  1:2409061/148‑1 (MQ=255)
                                                                   ccTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGc                                                                                 >  2:2986689/1‑148 (MQ=255)
                                                                    cTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCa                                                                                <  2:2179581/148‑1 (MQ=255)
                                                                    cTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGCTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCa                                                                                <  1:3276582/148‑1 (MQ=255)
                                                                        cGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCga                                                                            <  2:694252/148‑1 (MQ=255)
                                                                          gATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCgaga                                                                          >  1:1893094/1‑148 (MQ=255)
                                                                            ttAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAAc                                                                        <  2:731042/148‑1 (MQ=255)
                                                                                  taaatTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATa                                                                  <  1:2324654/148‑1 (MQ=255)
                                                                                   aaatTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATAc                                                                 >  1:138566/1‑148 (MQ=255)
                                                                                     atTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACaa                                                               <  2:442184/148‑1 (MQ=255)
                                                                                      ttCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAAc                                                              >  2:327759/1‑148 (MQ=255)
                                                                                         ttGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACttt                                                           <  1:477080/148‑1 (MQ=255)
                                                                                           gCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACttttt                                                         <  2:1883555/148‑1 (MQ=255)
                                                                                                gaaagaaaGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGa                                                    <  2:632695/148‑1 (MQ=255)
                                                                                                gaaagaaaGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGa                                                    <  2:3324382/148‑1 (MQ=255)
                                                                                                    gaaagaGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGt                                                <  2:3148847/148‑1 (MQ=255)
                                                                                                     aaagaGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTATTTACGGAACGTc                                               >  1:775900/1‑148 (MQ=255)
                                                                                                       agagGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCAc                                             <  1:1788342/148‑1 (MQ=255)
                                                                                                        gagGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCAcc                                            >  1:355797/1‑148 (MQ=255)
                                                                                                          ggAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTg                                          >  1:1824132/1‑148 (MQ=255)
                                                                                                          ggAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTg                                          <  1:1995304/148‑1 (MQ=255)
                                                                                                           gAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGa                                         >  1:1037767/1‑148 (MQ=255)
                                                                                                               tAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATaaa                                     >  1:57322/1‑148 (MQ=255)
                                                                                                               tAAACCATGAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATaaa                                     <  2:2317337/148‑1 (MQ=255)
                                                                                                                aaaCCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATaaaa                                    <  2:1217372/148‑1 (MQ=255)
                                                                                                                  aCCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAAtt                                  >  1:2740925/1‑148 (MQ=255)
                                                                                                                   ccATTAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAAttt                                 >  1:2526952/1‑148 (MQ=255)
                                                                                                                     aTTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGt                               >  1:1116113/1‑148 (MQ=255)
                                                                                                                      ttAGCGGAAAACCAGCGGCACGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGtt                              >  2:549952/1‑148 (MQ=255)
                                                                                                                       tAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTc                             >  1:1560120/1‑148 (MQ=255)
                                                                                                                         gCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAg                           >  2:1704694/1‑148 (MQ=255)
                                                                                                                           ggAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTa                         <  1:664231/148‑1 (MQ=255)
                                                                                                                             aaaaCCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATc                       >  1:1759985/1‑148 (MQ=255)
                                                                                                                              aaaCCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCg                      >  1:3117865/1‑148 (MQ=255)
                                                                                                                                aCCAGCGGCACGCCAGGGTGACATGCCCCAGTAAGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTc                    >  2:2564118/1‑148 (MQ=255)
                                                                                                                                 ccAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCGCCTGATAAAATTTGTTCAGTATCGTCt                   >  2:3099069/1‑148 (MQ=255)
                                                                                                                                  cAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCtt                  <  2:384895/148‑1 (MQ=255)
                                                                                                                                     cGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTg               <  1:2571233/148‑1 (MQ=255)
                                                                                                                                      ggCACGCCAGGGTGACCTGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACATTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTga              <  1:2245561/148‑1 (MQ=255)
                                                                                                                                      ggCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTga              <  1:1002068/148‑1 (MQ=255)
                                                                                                                                        cACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTgaga            <  2:1353143/148‑1 (MQ=255)
                                                                                                                                         aCGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTGAGAt           <  2:2714364/148‑1 (MQ=255)
                                                                                                                                          cGCCAGGGTTACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTGAGAtt          >  2:942802/1‑148 (MQ=255)
                                                                                                                                          cGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGAATGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTGAGAtt          <  2:570091/148‑1 (MQ=255)
                                                                                                                                            ccAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTGAGATtct        <  1:3233286/148‑1 (MQ=255)
                                                                                                                                             cAGGGTTACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTGAGATtctc       >  1:1280285/1‑148 (MQ=255)
                                                                                                                                             cAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCTCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTGAGATtctc       <  2:387895/148‑1 (MQ=255)
                                                                                                                                             cAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTGAGATtctc       >  2:1648904/1‑148 (MQ=255)
                                                                                                                                             cAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTGAGATtctc       <  2:3035212/148‑1 (MQ=255)
                                                                                                                                               gggTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTGAGATTCTCaa     >  2:2107647/1‑148 (MQ=255)
                                                                                                                                                ggTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTGAGATTCTCAAg    >  1:1489205/1‑148 (MQ=255)
                                                                                                                                                 gTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTGAGATTCTCAAgg   <  2:1908297/148‑1 (MQ=255)
                                                                                                                                                 gTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTGAGATTCTCAAgg   >  2:103582/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATATTCAGCCATGCATTAATTAAAATGTTACGTGTTTAATGTGAGTCCCTATCTGAAAATAAATAATCCTTCCGGATTAAAATAAATTCTTGCCGGGAAAGAAAGAGGAAATAAACCATTAGCGGAAAACCAGCGGCACGCCAGGGTGACATGCCCCAGTACGGCGGCCCGATAGTCCAGGGTTCAGCGGGGGTACTGATAGTCGCGCCGACTGCATCGAGAACTCCATACAACTTTTTACGGAACGTCACCTGATAAAATTTGTTCAGTATCGTCTTTTGAGATTCTCAAGGT  >  CP009273/524836‑525129

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: