Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,047,705 C→T 5.8% A71V (GCC→GTC)  quiA → quinate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,047,7050CT 5.8% 238.3 / 2.7 79A71V (GCC→GTC) quiAquinate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (39/35);  new base T (2/3);  total (41/38)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

CATCGCCCTGCTCGGCCTCCTCCTCGCTGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTACCCGCTGCTGGCGCGTGCCAACGGCGCCAGCG  >  NC_002947/4047556‑4047854
                                                                                                                                                     |                                                                                                                                                     
catcGCCCTGCTCGGCCTCCTCCTCGCTGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTATGGTCTGATCGGGGGCGTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGc                                                                                                                                                       >  1:529648/1‑150 (MQ=255)
 atcGCCCTGCTCGGCCTCCTCCTCGCTGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGcc                                                                                                                                                      >  1:187712/1‑150 (MQ=255)
  tcGCCCTGCTCGGCCTCCTCCTCGCTGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCGCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCg                                                                                                                                                     >  1:260450/1‑150 (MQ=255)
   cGCACTGCTCGTCCTCCTCCTCGCTGCCGGCGGCGTCAAGCTCTCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTGCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCgt                                                                                                                                                    <  2:260450/150‑1 (MQ=255)
     cccTGCTCGGCCTCCTCCTCGCTGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCgtgt                                                                                                                                                  >  2:1160237/1‑150 (MQ=255)
      ccTGCTCGGCCTCCTCCTCGCTGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGtt                                                                                                                                                 >  1:1402865/1‑150 (MQ=255)
          cTCGGCCTCCTCCTCGCTGCCGCGGGCCTCACACGGGCCGCGCGTGGCTGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTg                                                                                                                                             <  2:178091/150‑1 (MQ=255)
            cGGCCTCCTCCTCGCTGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGt                                                                                                                                           <  1:733264/150‑1 (MQ=255)
             ggCCTCCTCCTCGCTGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCTGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAAGCGGTCGCCTGGCTGTATGCCGTGTTCCTGGTc                                                                                                                                          >  2:1565/1‑150 (MQ=255)
             ggCCTCCTCCTCGCTGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTACTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTAGTc                                                                                                                                          >  2:668743/1‑150 (MQ=255)
                ctcctcctcGCTGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGt                                                                                                                                       >  2:1284450/1‑150 (MQ=255)
                 tcctcctcGCTGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCgg                                                                                                                                        <  2:1453893/148‑1 (MQ=255)
                 tcctcctcGCTGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCgg                                                                                                                                        >  1:1453893/1‑148 (MQ=255)
                     cctcGCTGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGc                                                                                                                                  >  2:1468259/1‑150 (MQ=255)
                       tcGCTGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCa                                                                                                                                >  1:1444287/1‑150 (MQ=255)
                         gCTGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATc                                                                                                                              >  2:1032049/1‑150 (MQ=255)
                           tGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGTCGTGTTCCTGGTc                                                                                                                                          <  1:1067554/136‑1 (MQ=255)
                           tGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGTCGTGTTCCTGGTc                                                                                                                                          <  1:1067386/136‑1 (MQ=255)
                           tGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGTCGTGTTCCTGGTc                                                                                                                                          >  2:1067554/1‑136 (MQ=255)
                           tGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGTCGTGTTCCTGGTc                                                                                                                                          >  2:1067386/1‑136 (MQ=255)
                           tGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTg                                                                                                                                             <  1:206431/133‑1 (MQ=255)
                           tGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTg                                                                                                                                             >  2:206431/1‑133 (MQ=255)
                            gCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCGTCTgg                                                                                                                           <  1:1183045/150‑1 (MQ=255)
                              cggcggCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGc                                                                                                                         <  2:529648/150‑1 (MQ=255)
                                           cTCGCCGGTCTGGGCGGTTCCTGGTCCTCGCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCg                                                                                                            >  1:43260/1‑150 (MQ=255)
                                                  gTCTGGGCGGTTCCTGGGACTGCCTGTTCGGCGGCTTGGCAATGGCCACCGCCGGTGTCCTCATTGCTCGCCGCACGAAGGCGGGCACCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCACCTGGCCTCTCATCGATGCGGGCCTgg                                                                                                     <  2:559393/150‑1 (MQ=255)
                                                    cTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGgtg                                                                                                   <  1:1284450/150‑1 (MQ=255)
                                                     tGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGAGTCCTCATTGCTCGCCGAAATTATGTGGGCGCCTGGATGTATGCCGTGTTCGTTGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGgtgt                                                                                                  >  1:366609/1‑150 (MQ=255)
                                                      gggCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGtt                                                                                                 >  1:544877/1‑150 (MQ=255)
                                                      gggCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGAGGGCGCATGGCTGTATATCGAGGAGCTGCTCGTTACTGGTATCTGGGCACTGATCGATGCGGGCCTGGTGtt                                                                                                 <  1:446103/150‑1 (MQ=255)
                                                        gCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCt                                                                                               <  1:461346/150‑1 (MQ=255)
                                                         cGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTg                                                                                              <  1:1032049/150‑1 (MQ=255)
                                                            ttCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGcc                                                                                           >  1:863996/1‑150 (MQ=255)
                                                                   tACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTtc                                                                                    <  2:400793/150‑1 (MQ=255)
                                                                      ttCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTAATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTcg                                                                                 >  1:135461/1‑150 (MQ=255)
                                                                       tCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTcgc                                                                                >  1:1393534/1‑150 (MQ=255)
                                                                        ccTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTcgcg                                                                               >  1:296241/1‑150 (MQ=255)
                                                                         cTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACt                                                                                                                      <  1:566421/110‑1 (MQ=255)
                                                                         cTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACt                                                                                                                      >  2:566421/1‑110 (MQ=255)
                                                                             tCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGt                                                                          >  1:479984/1‑150 (MQ=255)
                                                                             tCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGt                                                                          >  1:678408/1‑150 (MQ=255)
                                                                                 cggcTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCAt                                                                      <  1:775664/150‑1 (MQ=255)
                                                                                        gccatggccatCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATACCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGgc                                                               >  1:198689/1‑150 (MQ=255)
                                                                                         ccatggccatCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGgcg                                                              >  2:556182/1‑150 (MQ=255)
                                                                                            tggccatCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTccgctcca                                                           <  2:366609/150‑8 (MQ=255)
                                                                                             ggccatCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCg                                                                                                                                         >  1:391737/1‑71 (MQ=255)
                                                                                             ggccatCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCg                                                                                                                                         <  2:391737/71‑1 (MQ=255)
                                                                                                 atCGCCGGTGTCCTCATTGCTCGCCGTAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGc                                                      >  2:688444/1‑150 (MQ=255)
                                                                                                  tCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCa                                                     <  2:678408/150‑1 (MQ=255)
                                                                                                         tgtCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTc                                                                        <  2:845453/124‑1 (MQ=255)
                                                                                                         tgtCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTc                                                                        >  1:845453/1‑124 (MQ=255)
                                                                                                          gtCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTg                                             >  2:604708/1‑150 (MQ=255)
                                                                                                            ccTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGgc                                           <  1:1160237/150‑1 (MQ=255)
                                                                                                             cTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTg                                                                                                                                   <  1:1321570/61‑1 (MQ=255)
                                                                                                             cTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTg                                                                                                                                   >  2:1321570/1‑61 (MQ=255)
                                                                                                             cTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCACCTGGGCACTGATCGATGCGAGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGgcg                                          >  2:550676/1‑150 (MQ=255)
                                                                                                                aTTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCCGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTg                                       <  1:289254/150‑1 (MQ=255)
                                                                                                                 ttGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTgg                                      <  2:863996/150‑1 (MQ=255)
                                                                                                                     tCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGc                                                                                            >  2:1057464/1‑92 (MQ=255)
                                                                                                                     tCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGc                                                                                            <  1:1057464/92‑1 (MQ=255)
                                                                                                                     tCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTa                                  >  1:1131195/1‑150 (MQ=255)
                                                                                                                     tCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTa                                  >  1:1037285/1‑150 (MQ=255)
                                                                                                                             agaagGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTACCCGCTGC                          >  1:62754‑M1/1‑150 (MQ=255)
                                                                                                                               aagGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTACCCGCTGCTG                        <  1:604708‑M1/150‑1 (MQ=255)
                                                                                                                                  gcgggcgCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTAATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTACCCGCTGCTGGCG                     <  2:135461‑M1/150‑1 (MQ=255)
                                                                                                                                   cgggcgCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTACCCGCTGCTGGCGC                    <  2:479984‑M1/150‑1 (MQ=255)
                                                                                                                                      gcgcCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTACCCGCTGCTGGCGCGTg                 <  2:1402865/150‑1 (MQ=255)
                                                                                                                                        gcCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGc                                                                                            >  1:194288/1‑73 (MQ=255)
                                                                                                                                        gcCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGc                                                                                            <  2:194288/73‑1 (MQ=255)
                                                                                                                                          cTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTACCCGCTGCTGGCGCGTGCCaa             >  2:742736/1‑150 (MQ=255)
                                                                                                                                            ggCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTACCCGCTGCTGGCGCGTGCCAACg           >  1:185459/1‑150 (MQ=255)
                                                                                                                                            ggCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTACCCGCTGCTGGCGCGTGCCAACg           <  2:187712/150‑1 (MQ=255)
                                                                                                                                             gCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTAGAGCCTGTTCATGTACGACGGAAGCTGGGTGGTCGTGGCGCTGGTCTACCCGCTGCTGGCGCGTGCCAACGg          <  2:387583‑M2/150‑2 (MQ=255)
                                                                                                                                              cTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTACCCGCTGCTGGCGCGTGCCAACGgc         >  2:1077712‑M2/1‑148 (MQ=255)
                                                                                                                                               tGTATGCCGTGTTCCTGGTCGGTACTGCAATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTACCCGCTGCTGGCGCGTGCCAACGgcg        <  1:1563639‑M2/150‑4 (MQ=255)
                                                                                                                                                gTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGTGCCGGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTACCCGCCGTTGGCACGCGCCagcagcgg       <  2:43260/150‑8 (MQ=255)
                                                                                                                                                  aTGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTACCCGCTGCTGGCGCGTGCCAACGgcgcca     <  1:1260551‑M2/150‑7 (MQ=255)
                                                                                                                                                    gCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTc                                                <  1:947709/105‑1 (MQ=255)
                                                                                                                                                    gCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTc                                                >  2:947709/1‑105 (MQ=255)
                                                                                                                                                    gCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTACCCGCTGCTGGCGCGTGCCAACGgcgccagc   <  1:1014797‑M2/150‑9 (MQ=255)
                                                                                                                                                    gCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTACCCGCTGCTGGCGCGTGCCAACGgcgccagc   >  2:1444287‑M2/1‑142 (MQ=255)
                                                                                                                                                     ccGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTACCCGCTGCTGGCGCcgttggcacgcgccagca  <  1:69953‑M1/150‑19 (MQ=255)
                                                                                                                                                     ccGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTACCCGCTGCTGGCGCcgttggcacgcgccagca  >  2:69953‑M1/1‑132 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CATCGCCCTGCTCGGCCTCCTCCTCGCTGCCGGCGGCGTCAAGCTCGCCGGTCTGGGCGGTTCCTGGTACTTCCTGATCGGCGGCTTGGCCATGGCCATCGCCGGTGTCCTCATTGCTCGCCGCAAGAAGGCGGGCGCCTGGCTGTATGCCGTGTTCCTGGTCGGTACTGCCATCTGGGCACTGATCGATGCGGGCCTGGTGTTCTGGCCGCTGTTCTCGCGCCTGTTCATGTTCGGCGCAATCGGCATGGTCGTGGCGCTGGTCTACCCGCTGCTGGCGCGTGCCAACGGCGCCAGCG  >  NC_002947/4047556‑4047854

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: