Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,851,695 C→A 8.7% G24C (GGT→TGT)  yidJ ← sulfatase/phosphatase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,851,6950CA 8.7% 765.6 / 29.9 216G24C (GGT→TGT) yidJsulfatase/phosphatase superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (85/112);  new base A (12/7);  total (97/119)
Fisher's exact test for biased strand distribution p-value = 1.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 216 total aligned reads displayed.

GGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGA  >  CP009273/3851549‑3851840
                                                                                                                                                  |                                                                                                                                                 
ggTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCg                                                                                                                                                  <  1:3403668/148‑1 (MQ=255)
 gTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGc                                                                                                                                                 <  2:1310065/148‑1 (MQ=255)
    cACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCtat                                                                                                                                              <  2:2692428/148‑1 (MQ=255)
     aCGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCtata                                                                                                                                             <  2:655395/148‑1 (MQ=255)
     aCGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCtata                                                                                                                                             >  2:3444295/1‑148 (MQ=255)
      cGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGACTTAAAGCGAATACCTTCCGCCGCCAGACTATCAATCTTTTGCGTATTCAGCGGTTTACCGCTATAg                                                                                                                                            >  2:583625/1‑148 (MQ=255)
        ggCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACAGCTATAGCa                                                                                                                                          <  1:553827/148‑1 (MQ=255)
        ggCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACAGCTATAGCa                                                                                                                                          >  2:3757590/1‑148 (MQ=255)
         gCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCaa                                                                                                                                         >  2:2156247/1‑148 (MQ=255)
             gACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGa                                                                                                                                     >  1:3245935/1‑148 (MQ=255)
             gACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGa                                                                                                                                     <  1:2184789/148‑1 (MQ=255)
             gACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGAACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACAGCTATAGCAACCGa                                                                                                                                     <  2:813051/148‑1 (MQ=255)
                 ggTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCat                                                                                                                                 <  2:2117429/148‑1 (MQ=255)
                   ttAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCata                                                                                                                                >  1:3359707/1‑147 (MQ=255)
                   ttAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCata                                                                                                                                <  2:3359707/147‑1 (MQ=255)
                          aGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGc                                                                                                                        <  2:1673070/148‑1 (MQ=255)
                          aGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGc                                                                                                                        <  1:3668431/148‑1 (MQ=255)
                           gATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGcc                                                                                                                       >  2:3434560/1‑148 (MQ=255)
                           gATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGcc                                                                                                                       >  1:2831254/1‑148 (MQ=255)
                            aTACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACGGGTGAACTGGTGTAGGCGTAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCCGCGGTTTACCGATATAGCAACCGACCATATTGGTGGCCt                                                                                                                      <  1:3559087/148‑1 (MQ=255)
                            aTACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCt                                                                                                                      >  1:1363826/1‑148 (MQ=255)
                            aTACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCt                                                                                                                      <  1:2947076/148‑1 (MQ=255)
                            aTACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCt                                                                                                                      >  2:275461/1‑148 (MQ=255)
                            aTACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGATATAGCAACCGACCATATTGGTGGCCt                                                                                                                      <  1:24762/148‑1 (MQ=255)
                               ccGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTggg                                                                                                                   <  2:3406745/148‑1 (MQ=255)
                                cGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGt                                                                                                                  >  1:881957/1‑148 (MQ=255)
                                   tAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATc                                                                                                               <  1:3086851/148‑1 (MQ=255)
                                      aTAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGt                                                                                                            >  1:1771403/1‑148 (MQ=255)
                                         gTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCAt                                                                                                         <  2:3082816/148‑1 (MQ=255)
                                           ccGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGa                                                                                                       <  2:2047183/148‑1 (MQ=255)
                                              gCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGa                                                                                                    <  2:2962114/148‑1 (MQ=255)
                                              gCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGa                                                                                                    <  1:2534853/148‑1 (MQ=255)
                                               cccGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGaa                                                                                                   <  1:2602119/148‑1 (MQ=255)
                                               cccGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGaa                                                                                                   <  2:2402529/148‑1 (MQ=255)
                                                 cGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACa                                                                                                 >  2:685088/1‑148 (MQ=255)
                                                    gCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGaa                                                                                              <  1:1669364/148‑1 (MQ=255)
                                                       ggCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCAGATTGGTGGCCGGGGTCTCGGTCATGACGACCAGACAAt                                                                                           >  2:2646097/1‑148 (MQ=255)
                                                        gCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAAtt                                                                                          <  1:3649084/148‑1 (MQ=255)
                                                        gCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAAtt                                                                                          >  1:728723/1‑148 (MQ=255)
                                                        gCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAAtt                                                                                          <  2:410726/148‑1 (MQ=255)
                                                        gCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAAtc                                                                                          >  2:3784741/1‑147 (MQ=255)
                                                              aaaCCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCg                                                                                    <  1:719796/148‑1 (MQ=255)
                                                               aaCCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCACCATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGt                                                                                   <  2:2244558/148‑1 (MQ=255)
                                                               aaCCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGt                                                                                   >  2:3008414/1‑148 (MQ=255)
                                                                aCCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGtt                                                                                  >  2:2475327/1‑148 (MQ=255)
                                                                  cGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTc                                                                                <  1:3606765/148‑1 (MQ=255)
                                                                  cGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTc                                                                                <  1:1711719/148‑1 (MQ=255)
                                                                        aCAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGc                                                                          >  2:2871719/1‑148 (MQ=255)
                                                                        aCAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGc                                                                          <  1:1155931/148‑1 (MQ=255)
                                                                            gtgtAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATc                                                                      <  1:2280485/148‑1 (MQ=255)
                                                                             tgtAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATcc                                                                     <  2:3067101/148‑1 (MQ=255)
                                                                                aGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTc                                                                  <  2:2831228/148‑1 (MQ=255)
                                                                                     gAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTc                                                             >  1:1081816/1‑148 (MQ=255)
                                                                                        ttAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCggg                                                          <  1:941128/148‑1 (MQ=255)
                                                                                           aaGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCCTATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCTTTTCATAGGCTATCCTTCCCGTCGGGAGt                                                       <  2:1689150/148‑1 (MQ=255)
                                                                                           aaGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCTTTTCATAGGCTATCCTTCCCGTCGGGAGt                                                       <  2:2636686/148‑1 (MQ=255)
                                                                                           aaGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCTTTTCATAGGCTATCCTTCCCGTCGGGAGt                                                       <  2:2636675/148‑1 (MQ=255)
                                                                                             gCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTaa                                                     <  1:1525161/148‑1 (MQ=255)
                                                                                             gCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTaa                                                     <  1:429605/148‑1 (MQ=255)
                                                                                             gCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCAGATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTaa                                                     <  1:2236565/148‑1 (MQ=255)
                                                                                              cGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTaaa                                                    >  2:3634090/1‑148 (MQ=255)
                                                                                               gAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAg                                                   <  2:2206240/148‑1 (MQ=255)
                                                                                               gAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAg                                                   <  2:295462/148‑1 (MQ=255)
                                                                                                aaTACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATCTTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAgg                                                  <  1:2646097/148‑1 (MQ=255)
                                                                                                aaTACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAgg                                                  >  2:786983/1‑148 (MQ=255)
                                                                                                aaTACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAgg                                                  >  2:787034/1‑148 (MQ=255)
                                                                                                aaTACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAgg                                                  >  2:3302558/1‑148 (MQ=255)
                                                                                                aaTACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAgg                                                  <  1:592582/148‑1 (MQ=255)
                                                                                                aaTACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAgg                                                  >  1:579780/1‑148 (MQ=255)
                                                                                                 aTACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGt                                                 <  2:1311216/148‑1 (MQ=255)
                                                                                                 aTACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGt                                                 >  1:1947308/1‑148 (MQ=255)
                                                                                                 aTACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGt                                                 >  1:2430450/1‑148 (MQ=255)
                                                                                                   aCCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTa                                               >  2:1531296/1‑148 (MQ=255)
                                                                                                   aCCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTa                                               <  1:405122/148‑1 (MQ=255)
                                                                                                       tCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGc                                           <  2:1291353/148‑1 (MQ=255)
                                                                                                                aGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTg                                  <  2:2940767/148‑1 (MQ=255)
                                                                                                                 gACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGa                                 >  1:2816435/1‑148 (MQ=255)
                                                                                                                   cTATCAATATTTTGCGTATTCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACg                               <  2:1039870/148‑1 (MQ=255)
                                                                                                                    tATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGa                              >  2:2925484/1‑148 (MQ=255)
                                                                                                                      tCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATg                            >  1:1012867/1‑148 (MQ=255)
                                                                                                                      tCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATg                            >  2:1537412/1‑148 (MQ=255)
                                                                                                                      tCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATg                            >  1:1999301/1‑148 (MQ=255)
                                                                                                                       cAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATgc                           >  1:3544878/1‑148 (MQ=255)
                                                                                                                       cAATATTTTGCGTATTCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGACCAGAAAATTGGGGCGTTTCATAGACTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATgc                           >  2:1516119/1‑148 (MQ=255)
                                                                                                                       cAATATTTTGCGTATTCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGACTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATgc                           >  2:2273605/1‑148 (MQ=255)
                                                                                                                         atatTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATgccg                         >  1:1197259/1‑148 (MQ=255)
                                                                                                                         atatTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATgccg                         >  2:1731976/1‑148 (MQ=255)
                                                                                                                             tttGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAg                     >  1:2309768/1‑148 (MQ=255)
                                                                                                                               tGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCt                   <  2:1494292/148‑1 (MQ=255)
                                                                                                                                     ttCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCaaa             >  2:3223875/1‑148 (MQ=255)
                                                                                                                                      tCAGCGGTTTACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCaaaa            <  1:1800733/148‑1 (MQ=255)
                                                                                                                                       cAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCaaaaa           <  2:3786761/148‑1 (MQ=255)
                                                                                                                                        aGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAAt          <  2:1159319/148‑1 (MQ=255)
                                                                                                                                          cGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTAAGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCa        >  2:3217967/1‑148 (MQ=255)
                                                                                                                                           ggTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAg       <  2:3030920/148‑1 (MQ=255)
                                                                                                                                            gTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAgg      <  2:1363826/148‑1 (MQ=255)
                                                                                                                                             tttACAGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGt     >  2:725563/1‑148 (MQ=255)
                                                                                                                                               tACCGCTATAGCAACCGACCATAGTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAg   <  2:3649970/148‑1 (MQ=255)
                                                                                                                                                aCCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGa  <  2:368897/148‑1 (MQ=255)
                                                                                                                                                aCCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGa  <  2:2016551/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GGTCCACGGGCCGGACTGGTTAGCGTAGATACCGGTAAATAGTCCGGCCCGTGCAGGCGTACAAACCGGTGAACAGGTGTAGGCGGAATTAAAGCGAATACCTTCCGCCGCCAGACTATCAATATTTTGCGTATTCAGCGGTTTACCGCTATAGCAACCGACCATATTGGTGGCCTGGGTATCGGTCATGACGAACAGAAAATTGGGGCGTTTCATAGGCTATCCTTCCCGTCGGGAGTAAAGGTTACGGCTGGGTCGTGACGATGCCGCCAGCTGTCAAAAATCAGGTAGA  >  CP009273/3851549‑3851840

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: