Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,913,709 C→A 15.6% E133* (GAA→TAA)  gudX ← glucarate dehydratase‑related protein, substrate unknown

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,913,7090CA15.6% 726.0 / 76.6 237E133* (GAA→TAA) gudXglucarate dehydratase‑related protein, substrate unknown
Reads supporting (aligned to +/- strand):  ref base C (112/88);  new base A (22/15);  total (134/103)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

Only 100 of 240 total aligned reads displayed.

TTGCTTTCTGATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTG  >  CP009273/2913562‑2913853
                                                                                                                                                   |                                                                                                                                                
ttGCTTTCTGATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTc                                                                                                                                                  >  2:1941167/1‑148 (MQ=255)
   cTTTCTGATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCa                                                                                                                                               >  2:2217557/1‑148 (MQ=255)
    tttCTGATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAg                                                                                                                                              >  1:2816429/1‑148 (MQ=255)
     ttCTGATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGa                                                                                                                                             >  1:1030466/1‑148 (MQ=255)
      tCTGATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGAc                                                                                                                                            <  2:2245568/148‑1 (MQ=255)
         gATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCgg                                                                                                                                         >  1:1644126/1‑148 (MQ=255)
         gATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCgg                                                                                                                                         >  2:1475902/1‑148 (MQ=255)
         gATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCgg                                                                                                                                         <  1:346398/148‑1 (MQ=255)
                  aCTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAg                                                                                                                                >  2:3483968/1‑148 (MQ=255)
                     gATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAgcgc                                                                                                                             >  2:1340945/1‑148 (MQ=255)
                       tACCACTCATGGTTGCCCGGCGTATTTTCCACATTAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCt                                                                                                                           <  1:1211608/148‑1 (MQ=255)
                       tACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGACCATTCAGCGCCt                                                                                                                           >  2:491845/1‑148 (MQ=255)
                         ccACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTa                                                                                                                         >  1:116709/1‑148 (MQ=255)
                           aCTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTAcc                                                                                                                       <  1:3196699/148‑1 (MQ=255)
                            cTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCt                                                                                                                      >  2:1875389/1‑148 (MQ=255)
                            cTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCt                                                                                                                      <  1:2808100/148‑1 (MQ=255)
                             tCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTa                                                                                                                     >  1:1236959/1‑148 (MQ=255)
                              cATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAg                                                                                                                    >  2:2808223/1‑148 (MQ=255)
                              cATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAg                                                                                                                    >  2:283801/1‑148 (MQ=255)
                                tGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCa                                                                                                                  >  1:1772701/1‑148 (MQ=255)
                                tGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCa                                                                                                                  >  1:865147/1‑148 (MQ=255)
                                  gTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAgg                                                                                                                >  2:660600/1‑148 (MQ=255)
                                    tGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTc                                                                                                              >  2:266065/1‑148 (MQ=255)
                                     gCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCa                                                                                                             <  2:404228/148‑1 (MQ=255)
                                        cGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGc                                                                                                          <  1:1460580/148‑1 (MQ=255)
                                         ggCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCa                                                                                                         >  2:2721696/1‑148 (MQ=255)
                                          gCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCaa                                                                                                        >  2:234213/1‑148 (MQ=255)
                                           cGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCaaa                                                                                                       >  2:587800/1‑148 (MQ=255)
                                              aTTTTCCACATAAGGAAGGTCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCg                                                                                                    >  1:3280246/1‑148 (MQ=255)
                                                tttCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGc                                                                                                  <  2:289618/148‑1 (MQ=255)
                                                  tCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCtt                                                                                                <  1:907589/148‑1 (MQ=255)
                                                   ccACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTc                                                                                               <  2:1005268/148‑1 (MQ=255)
                                                     acaTAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCa                                                                                             >  1:3134764/1‑148 (MQ=255)
                                                         aaGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAgcgc                                                                                         >  2:1773574/1‑148 (MQ=255)
                                                         aaGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGCAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAgcgc                                                                                         >  1:636499/1‑148 (MQ=255)
                                                              aGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCa                                                                                    <  2:3514820/148‑1 (MQ=255)
                                                              aGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCa                                                                                    <  1:1630714/148‑1 (MQ=255)
                                                                 tCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACgg                                                                                 >  1:2797556/1‑148 (MQ=255)
                                                                 tCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACgg                                                                                 >  2:1153886/1‑148 (MQ=255)
                                                                     ttttGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGtt                                                                             >  1:1344979/1‑148 (MQ=255)
                                                                           tCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAAcgcg                                                                       >  1:441566/1‑148 (MQ=255)
                                                                           tCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAAcgcg                                                                       >  2:1984428/1‑148 (MQ=255)
                                                                           tCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAAcgcg                                                                       <  2:2719930/148‑1 (MQ=255)
                                                                              gATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTaa                                                                    >  1:552585/1‑148 (MQ=255)
                                                                                tCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAAtt                                                                  <  1:268402/148‑1 (MQ=255)
                                                                                   ccGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCaa                                                               <  2:1805628/148‑1 (MQ=255)
                                                                                    cGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCaaa                                                              <  2:1569393/148‑1 (MQ=255)
                                                                                    cGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCaaa                                                              >  1:1061628/1‑148 (MQ=255)
                                                                                     gATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAg                                                             <  2:1237563/148‑1 (MQ=255)
                                                                                     gATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAg                                                             >  1:246755/1‑148 (MQ=255)
                                                                                       tataAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGt                                                           >  2:1882655/1‑148 (MQ=255)
                                                                                        ataAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTc                                                          <  1:283790/148‑1 (MQ=255)
                                                                                        ataAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTc                                                          >  2:1497786/1‑148 (MQ=255)
                                                                                        ataAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTc                                                          <  1:2274758/148‑1 (MQ=255)
                                                                                          aaaaCAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCa                                                        >  2:2197487/1‑148 (MQ=255)
                                                                                          aaaaCAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCa                                                        >  2:927702/1‑148 (MQ=255)
                                                                                          aaaaCAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCa                                                        >  2:2705656/1‑148 (MQ=255)
                                                                                           aaaCAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAg                                                       >  1:2922214/1‑148 (MQ=255)
                                                                                           aaaCAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAg                                                       >  1:2881479/1‑148 (MQ=255)
                                                                                           aaaCAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAg                                                       >  1:1547083/1‑148 (MQ=255)
                                                                                              cAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCa                                                    >  2:2223419/1‑148 (MQ=255)
                                                                                              cAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCa                                                    >  2:2499747/1‑148 (MQ=255)
                                                                                              cAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCa                                                    <  2:962442/148‑1 (MQ=255)
                                                                                               aGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCAc                                                   <  2:3362639/148‑1 (MQ=255)
                                                                                                gATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCAcc                                                  <  1:1417226/148‑1 (MQ=255)
                                                                                                 aTAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCt                                                 <  2:1542576/148‑1 (MQ=255)
                                                                                                 aTAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCt                                                 >  2:2913977/1‑148 (MQ=255)
                                                                                                 aTAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCt                                                 <  1:1032564/148‑1 (MQ=255)
                                                                                                  tAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCtt                                                >  2:1425331/1‑148 (MQ=255)
                                                                                                   aaCCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCttt                                               >  1:3031227/1‑148 (MQ=255)
                                                                                                    aCCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTa                                              >  1:3341531/1‑148 (MQ=255)
                                                                                                    aCCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTa                                              >  1:1582986/1‑148 (MQ=255)
                                                                                                     ccGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTAc                                             >  2:3223718/1‑148 (MQ=255)
                                                                                                         ggACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGaa                                         >  2:2968285/1‑148 (MQ=255)
                                                                                                          gACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAg                                        >  2:1235309/1‑148 (MQ=255)
                                                                                                            cGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGt                                      <  2:2759659/148‑1 (MQ=255)
                                                                                                             ggTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTa                                     <  2:3197306/148‑1 (MQ=255)
                                                                                                             ggTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTa                                     <  1:1941167/148‑1 (MQ=255)
                                                                                                                aaTAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCa                                  >  1:2112871/1‑148 (MQ=255)
                                                                                                                   aGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCaaaa                               <  1:717211/148‑1 (MQ=255)
                                                                                                                       tCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATcgg                           <  2:2833420/148‑1 (MQ=255)
                                                                                                                       tCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATcgg                           <  2:1058265/148‑1 (MQ=255)
                                                                                                                        cgcgTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATcggc                          >  1:359324/1‑148 (MQ=255)
                                                                                                                          cgTTGCTTGCCTGGCCCTAACAGTTCTCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGACACGGCGTTAACGCGTAATACAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAAAcggcgg                        >  1:2180929/1‑148 (MQ=255)
                                                                                                                          cgTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATcggcgg                        >  1:2584170/1‑148 (MQ=255)
                                                                                                                            ttgcttgcCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGcc                      >  1:3014130/1‑148 (MQ=255)
                                                                                                                              gcttgcCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTg                    <  2:3019576/148‑1 (MQ=255)
                                                                                                                               cttgcCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGa                   <  1:2628410/148‑1 (MQ=255)
                                                                                                                                 tgcCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGAtt                 >  1:2742947/1‑148 (MQ=255)
                                                                                                                                  gcCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTa                <  1:587800/148‑1 (MQ=255)
                                                                                                                                   ccTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTAc               <  1:2423895/148‑1 (MQ=255)
                                                                                                                                   ccTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTAc               <  2:2440882/148‑1 (MQ=255)
                                                                                                                                   ccTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTAc               <  1:1207027/148‑1 (MQ=255)
                                                                                                                                    cTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTAcc              <  1:3081514/148‑1 (MQ=255)
                                                                                                                                      ggCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGNCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCtt            >  1:1802702/1‑148 (MQ=255)
                                                                                                                                      ggCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCtt            <  1:1875389/148‑1 (MQ=255)
                                                                                                                                       gCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCttt           >  1:115703/1‑148 (MQ=255)
                                                                                                                                       gCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCttt           >  2:1652510/1‑148 (MQ=255)
                                                                                                                                           tAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGc       >  1:554273/1‑148 (MQ=255)
                                                                                                                                             aCAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCAc     >  1:1771665/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TTGCTTTCTGATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTG  >  CP009273/2913562‑2913853

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: