Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,913,709 C→A 7.9% 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,7090CA 7.9% 960.4 / 34.8 265E133* (GAA→TAA) gudXglucarate dehydratase‑related protein, substrate unknown
Reads supporting (aligned to +/- strand):  ref base C (135/109);  new base A (8/13);  total (143/122)
Fisher's exact test for biased strand distribution p-value = 1.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 266 total aligned reads displayed.

TTGCTTTCTGATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCT  >  CP009273/2913562‑2913855
                                                                                                                                                   |                                                                                                                                                  
 tGCTTTCTGATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCg                                                                                                                                                   >  2:1604009/1‑148 (MQ=255)
      tCTGATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGAc                                                                                                                                              <  2:275611/148‑1 (MQ=255)
          aTGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGa                                                                                                                                          <  2:2943585/148‑1 (MQ=255)
             ccGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACa                                                                                                                                       >  1:940027/1‑148 (MQ=255)
              cGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACAt                                                                                                                                      >  2:828095/1‑148 (MQ=255)
               gCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACAtt                                                                                                                                     <  2:395645/148‑1 (MQ=255)
                 aaCTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCa                                                                                                                                   <  1:1824614/148‑1 (MQ=255)
                 aaCTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCa                                                                                                                                   <  1:1739318/148‑1 (MQ=255)
                   cTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAgc                                                                                                                                 >  2:225282/1‑148 (MQ=255)
                   cTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAgc                                                                                                                                 >  2:452539/1‑148 (MQ=255)
                      aTACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGcc                                                                                                                              <  1:468292/148‑1 (MQ=255)
                       tACCACTCATGGTTGCCCGGAGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCt                                                                                                                             <  1:1222715/148‑1 (MQ=255)
                        aCCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCtt                                                                                                                            >  2:166080/1‑148 (MQ=255)
                            cTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCt                                                                                                                        <  1:2219910/148‑1 (MQ=255)
                              cATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAg                                                                                                                      >  2:1722967/1‑148 (MQ=255)
                               aTGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGc                                                                                                                     <  2:1334498/148‑1 (MQ=255)
                                tGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCa                                                                                                                    <  2:2162823/148‑1 (MQ=255)
                                  gTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAgg                                                                                                                  <  1:584626/148‑1 (MQ=255)
                                          gCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGACTGGCCCTAACAGTTCGCATACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCaa                                                                                                          >  1:830744/1‑148 (MQ=255)
                                           cGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCaaa                                                                                                         <  2:1535860/148‑1 (MQ=255)
                                             tATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGc                                                                                                       <  2:3390167/148‑1 (MQ=255)
                                               ttttCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCgg                                                                                                     <  1:1372912/148‑1 (MQ=255)
                                                        tAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAgcg                                                                                            >  1:761723/1‑148 (MQ=255)
                                                           ggAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAAGAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGcgccg                                                                                         <  2:2995215/148‑1 (MQ=255)
                                                            gAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGcgccgc                                                                                        <  1:285962/148‑1 (MQ=255)
                                                                aTCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACg                                                                                    >  1:1692430/1‑148 (MQ=255)
                                                                aTCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACg                                                                                    >  2:1746752/1‑148 (MQ=255)
                                                                 tCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCGCgg                                                                                   >  1:1476042/1‑148 (MQ=255)
                                                                   ggTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCg                                                                                 <  1:1317157/148‑1 (MQ=255)
                                                                   ggTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCg                                                                                 <  2:630099/148‑1 (MQ=255)
                                                                       ttGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTaa                                                                             >  2:1766646/1‑148 (MQ=255)
                                                                          gTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAAcgc                                                                          >  1:2337727/1‑148 (MQ=255)
                                                                          gTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAAcgc                                                                          >  1:2337674/1‑148 (MQ=255)
                                                                          gTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAAcgc                                                                          >  1:1840704/1‑148 (MQ=255)
                                                                          gTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAAcgc                                                                          >  1:772302/1‑148 (MQ=255)
                                                                           tCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAAcgcg                                                                         <  2:2141130/148‑1 (MQ=255)
                                                                            ccGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGt                                                                        <  1:386702/148‑1 (MQ=255)
                                                                            ccGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGt                                                                        <  2:1943865/148‑1 (MQ=255)
                                                                            ccGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGt                                                                        <  2:1944264/148‑1 (MQ=255)
                                                                            ccGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGt                                                                        >  2:1820201/1‑148 (MQ=255)
                                                                              gATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTAGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTaa                                                                      <  1:2992809/148‑1 (MQ=255)
                                                                                  aCCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCa                                                                  >  2:1359965/1‑148 (MQ=255)
                                                                                    cGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCaaa                                                                <  1:3460613/148‑1 (MQ=255)
                                                                                    cGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCaaa                                                                >  2:1237723/1‑148 (MQ=255)
                                                                                     gATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAg                                                               >  1:2614870/1‑148 (MQ=255)
                                                                                      atatAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAgg                                                              >  1:3313739/1‑148 (MQ=255)
                                                                                       tataAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGt                                                             >  2:1393472/1‑148 (MQ=255)
                                                                                       tataAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGt                                                             >  2:3556034/1‑148 (MQ=255)
                                                                                       tataAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGt                                                             >  1:2681032/1‑148 (MQ=255)
                                                                                        ataAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTc                                                            >  2:1003041/1‑148 (MQ=255)
                                                                                        ataAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGATCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTc                                                            >  2:2785926/1‑148 (MQ=255)
                                                                                            aaCAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAgg                                                        >  2:1058862/1‑148 (MQ=255)
                                                                                            aaCAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAgg                                                        >  2:3211876/1‑148 (MQ=255)
                                                                                                gATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCAcc                                                    <  2:940011/148‑1 (MQ=255)
                                                                                                   aaCCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCATAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCttt                                                 <  1:1287653/148‑1 (MQ=255)
                                                                                                   aaCCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCttt                                                 <  2:2426608/148‑1 (MQ=255)
                                                                                                   aaCCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCttt                                                 <  1:1107470/148‑1 (MQ=255)
                                                                                                    aCCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCCCGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTa                                                >  2:2662858/1‑148 (MQ=255)
                                                                                                     ccGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTAc                                               >  2:3212029/1‑148 (MQ=255)
                                                                                                      cGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCCCGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTAcc                                              <  2:783170/148‑1 (MQ=255)
                                                                                                        aGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGa                                            <  1:700335/148‑1 (MQ=255)
                                                                                                          gACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAg                                          <  1:2264375/148‑1 (MQ=255)
                                                                                                          gACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAg                                          <  2:534841/148‑1 (MQ=255)
                                                                                                            cGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGt                                        <  1:2805462/148‑1 (MQ=255)
                                                                                                            cGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGt                                        <  1:1560936/148‑1 (MQ=255)
                                                                                                               tAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATc                                     >  1:970770/1‑148 (MQ=255)
                                                                                                                 aTAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCaa                                   <  2:733613/148‑1 (MQ=255)
                                                                                                                  tAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCaaa                                  >  2:1106343/1‑148 (MQ=255)
                                                                                                                  tAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCaaa                                  >  2:2884950/1‑148 (MQ=255)
                                                                                                                   aGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAATGTCCAGTCACATTTACCGAAGGTATCaaaa                                 >  1:312502/1‑148 (MQ=255)
                                                                                                                   aGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCaaaa                                 >  1:1162286/1‑148 (MQ=255)
                                                                                                                     ccTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATc                               >  1:3517234/1‑148 (MQ=255)
                                                                                                                      cTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAACGGTCCAGGCACCTTTACCGAAGGTATCAAAATcg                              <  1:828095/148‑1 (MQ=255)
                                                                                                                      cTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATcg                              >  2:2046537/1‑148 (MQ=255)
                                                                                                                      cTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATcg                              <  2:2529661/148‑1 (MQ=255)
                                                                                                                      cTCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATcg                              <  2:2129359/148‑1 (MQ=255)
                                                                                                                       tCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATcgg                             <  1:1059891/148‑1 (MQ=255)
                                                                                                                       tCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATcgg                             <  1:1058562/148‑1 (MQ=255)
                                                                                                                       tCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATcgg                             >  2:1441698/1‑148 (MQ=255)
                                                                                                                       tCGCGTTGCTTGCCTGGCCCTAACAGTTAGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATcgg                             >  1:1853105/1‑148 (MQ=255)
                                                                                                                        cgcgTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATcggc                            >  1:237650/1‑148 (MQ=255)
                                                                                                                              gcttgcCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTg                      >  2:2452007/1‑148 (MQ=255)
                                                                                                                              gcttgcCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTg                      >  2:663755/1‑148 (MQ=255)
                                                                                                                              gcttgcCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTg                      <  1:1909535/148‑1 (MQ=255)
                                                                                                                               cttgcCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGa                     >  2:1783801/1‑148 (MQ=255)
                                                                                                                               cttgcCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGa                     <  1:166080/148‑1 (MQ=255)
                                                                                                                                  gcCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTa                  <  1:3372504/148‑1 (MQ=255)
                                                                                                                                    cTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTAcc                >  1:2499036/1‑148 (MQ=255)
                                                                                                                                     tGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCt               >  2:964703/1‑148 (MQ=255)
                                                                                                                                     tGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCt               <  2:3150144/148‑1 (MQ=255)
                                                                                                                                      ggCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCtt              >  2:3076731/1‑148 (MQ=255)
                                                                                                                                       gCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTCATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCttt             >  2:1121133/1‑148 (MQ=255)
                                                                                                                                       gCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCttt             >  2:2294648/1‑148 (MQ=255)
                                                                                                                                        cccTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTa            <  1:2087499/148‑1 (MQ=255)
                                                                                                                                            aaCAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCa        >  1:1457559/1‑148 (MQ=255)
                                                                                                                                             aCAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCAc       >  1:3556767/1‑148 (MQ=255)
                                                                                                                                             aCAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCAc       >  1:1993731/1‑148 (MQ=255)
                                                                                                                                             aCAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCAc       >  1:193294/1‑148 (MQ=255)
                                                                                                                                                 ttCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACctgc   >  2:1283224/1‑148 (MQ=255)
                                                                                                                                                  tCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACctgct  >  2:1879027/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TTGCTTTCTGATGCCGCAACTGATACCACTCATGGTTGCCCGGCGTATTTTCCACATAAGGAAGATCGGTTTTGGTCCGATCACCGATATAAAACAGATAACCGAGGACGGTAATAGCCTCGCGTTGCTTGCCTGGCCCTAACAGTTCGCAGACCGGAACATTCAGCGCCTTACCTAGCAGGTCAAGCAAAGCGGCTTCCAGCGCCGCCACGGCGTTAACGCGTAATTCAAAGGTCCAGGCACCTTTACCGAAGGTATCAAAATCGGCGGCCTGATTACCTTTATGCACCTGCT  >  CP009273/2913562‑2913855

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: