Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,861,908 T→A 5.1% V62V (GTA→GTT)  gllC ← 4‑carboxy‑4‑hydroxy‑2‑oxoadipic acid aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,861,9080TA 5.1% 250.9 / 3.1 78V62V (GTA→GTT) gllC4‑carboxy‑4‑hydroxy‑2‑oxoadipic acid aldolase
Reads supporting (aligned to +/- strand):  ref base T (38/36);  new base A (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

TGGAGCGAGGTCGCCAGCAGGTCGCCGAAATAGCCATCGCTGCACGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCGATCATGTGTGGCTCAACGCGTGGGATATTA  >  NC_002947/2861768‑2862055
                                                                                                                                            |                                                                                                                                                   
tGGAGCGAGGTCGCCAGCAGGTCGCCGAAATAGCCATCGCTGCACGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTg                                                                                                                                            <  2:1820697/150‑1 (MQ=255)
    ccgatctCGCCAGCAGGTCGCCGAAATAGCCATCGCTGCACGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGNANCATCCANTTGTCGCCTGGGNCCACCAGCGCGGTAACCGTGCTGGNgg                                                                                                                                        <  1:1460147/144‑1 (MQ=255)
             ccAGCAGGTCGCCGAAATAGCCATCGCTGCACGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCca                                                                                                                               <  1:607593/150‑1 (MQ=255)
                    gTCGCCGAAATAGCCATCGCTGCACGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTg                                                                                                                        >  1:394430/1‑150 (MQ=255)
                      cGCCGAAATAGCCATCGCTGCACGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCt                                                                                                                      <  2:394430/150‑1 (MQ=255)
                      cGCCGAAATAGCCATCGCTGCACGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCt                                                                                                                      <  2:1160328/150‑1 (MQ=255)
                       gCCGAAATAGCCATCGCTGCACGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTg                                                                                                                     >  1:1221943/1‑150 (MQ=255)
                          gAAATAGCCATCGCTGCACGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGcc                                                                                                                                          >  1:1146916/1‑126 (MQ=255)
                          gAAATAGCCATCGCTGCACGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGcc                                                                                                                                          <  2:1146916/126‑1 (MQ=255)
                          gAAATAGCCATCGCTGCACGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGAt                                                                                                                  <  1:597299/150‑1 (MQ=255)
                             aTAGCCATCGCTGCACGGCGAGCTGGGCGCCACCACGAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAgg                                                                                                               >  2:1096820/1‑150 (MQ=255)
                                 ccATCGCTGCACGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGgcgc                                                                                                           >  2:317869/1‑150 (MQ=255)
                                   aTCGCTGCACGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCAc                                                                                                         >  1:1350602/1‑150 (MQ=255)
                                   aTCGCTGCACGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCAc                                                                                                         >  1:230334/1‑150 (MQ=255)
                                     cGCTGCACGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACgg                                                                                                       <  2:1298272/150‑1 (MQ=255)
                                      gcggCACGGCGACCTGGGCGGCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCATGGc                                                                                                      <  2:262067/147‑1 (MQ=255)
                                      gCTGCACGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGc                                                                                                      <  2:261930/150‑1 (MQ=255)
                                       cTGCACGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCa                                                                                                     <  2:1350602/150‑1 (MQ=255)
                                           aCGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTAt                                                                                                 >  2:2099178/1‑150 (MQ=255)
                                            cGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTAtt                                                                                                <  1:1476329/150‑1 (MQ=255)
                                                                ccaGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACCTCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGAACtt                                                                            <  1:872269/150‑1 (MQ=255)
                                                                       gTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGa                                                                                              <  1:605819/125‑1 (MQ=255)
                                                                       gTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGa                                                                                              >  2:605819/1‑125 (MQ=255)
                                                                              gggCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCAcc                                                              <  2:1949136/150‑1 (MQ=255)
                                                                                gcggcACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGt                                                            >  1:1593128/1‑150 (MQ=255)
                                                                                  ggcACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTgg                                                          >  1:2355267/1‑150 (MQ=255)
                                                                                     aCTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCt                                                                                                                      <  2:895651/87‑1 (MQ=255)
                                                                                     aCTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCt                                                                                                                      >  1:895651/1‑87 (MQ=255)
                                                                                     aCTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTg                                                           <  2:1916860/146‑1 (MQ=255)
                                                                                     aCTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTg                                                           >  1:1916860/1‑146 (MQ=255)
                                                                                      cTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGac                                                      >  1:1237342/1‑150 (MQ=255)
                                                                                       tGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGaca                                                     <  2:1593128/150‑1 (MQ=255)
                                                                                       tGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGaca                                                     <  1:2091708/150‑1 (MQ=255)
                                                                                        gCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGacac                                                    >  2:116747/1‑150 (MQ=255)
                                                                                          tCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCCAATGCATCGTGGCGacactc                                                  >  2:674208/1‑148 (MQ=255)
                                                                                           ccACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCa                                                 <  1:799881/150‑1 (MQ=255)
                                                                                            cACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAg                                                >  2:534268/1‑150 (MQ=255)
                                                                                               ggCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCg                                             <  2:1221943/150‑1 (MQ=255)
                                                                                                gCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCgg                                            >  2:1416814/1‑150 (MQ=255)
                                                                                                    cagggAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCAACCTGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGCGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCa                                        >  1:442151/4‑150 (MQ=255)
                                                                                                     aTGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATa                                                                                                                 >  2:1709351/1‑76 (MQ=255)
                                                                                                     aTGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATa                                                                                                                 <  1:1709351/76‑1 (MQ=255)
                                                                                                     aTGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAg                                       >  2:2039443/1‑150 (MQ=255)
                                                                                                      tGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCGACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGTACCGTGGCGACACCCAGCCGGCCCAgg                                      >  1:243643/1‑150 (MQ=255)
                                                                                                      tGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAgg                                      <  2:2355267/150‑1 (MQ=255)
                                                                                                        gAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGc                                    >  2:576347/1‑150 (MQ=255)
                                                                                                         aaCATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCa                                   >  1:1005392/1‑150 (MQ=255)
                                                                                                            aTCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCTACACCCAGCCAGGCCAGGCCAGCg                                >  2:911645/1‑150 (MQ=255)
                                                                                                            aTCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCg                                >  1:1134690/1‑150 (MQ=255)
                                                                                                             tCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGagcagc                                                                                         <  1:1068451/92‑1 (MQ=255)
                                                                                                             tCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGagcagc                                                                                         >  2:1068451/1‑92 (MQ=255)
                                                                                                             tCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCGa                               >  2:1800691/1‑150 (MQ=255)
                                                                                                               cAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTCGCGACACCCACCCGCCCCACGCCCTCGATc                             >  2:1340715/1‑150 (MQ=255)
                                                                                                                    gTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCtt                                                                                    <  1:2076120/90‑1 (MQ=255)
                                                                                                                    gTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCtt                                                                                    >  2:2076120/1‑90 (MQ=255)
                                                                                                                    gTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCGATCAtgtg                        <  2:1237342/150‑1 (MQ=255)
                                                                                                                     gcGCCTGGTGCCACCAATACGGTTACGGCGGTGCTGGTCCTCGCCCGACCTCGCTGGATAGGGCGCACGGACGTATGGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCGATCAtgtgt                       <  1:2340694/149‑1 (MQ=255)
                                                                                                                      cGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCGCGCGGCCCAGGGCATCGATCAtgtgtg                      <  1:1271314/150‑1 (MQ=255)
                                                                                                                             ggCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCCGCCGGCCCAGGGCATCGATCATGTGTGGCTCAAc               <  2:540691/150‑1 (MQ=255)
                                                                                                                             ggCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGTCGGCCCAGGGCATCGATCATGTGTGGCTCAAc               >  1:2214385/1‑150 (MQ=255)
                                                                                                                              gCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGCCACCCAGCCGGCCCAGGGCATCGATCATGTGTGGCTCAAcg              >  1:996985/1‑150 (MQ=255)
                                                                                                                               ccggcaGCGCGGTAACCGTGCTGGTGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTgg                                                          >  2:369448/5‑105 (MQ=37)
                                                                                                                               ccggcaGCGCGGTAACCGTGCTGGTGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTgg                                                          <  1:369448/101‑1 (MQ=37)
                                                                                                                                caccaGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCGATCATGTGTGGCTCAAcgcg            >  2:1512301/1‑150 (MQ=255)
                                                                                                                                 accaGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCGATCATGTGTGGCTCAACGCGt           >  1:1129825/1‑150 (MQ=255)
                                                                                                                                 accaGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCGATCATGTGTGGCTCAACGCGt           <  2:1005392/150‑1 (MQ=255)
                                                                                                                                  gcaGCGCGGTAACCGTGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCGATCATGTGTGGCTCAACGCGTg          >  1:561882/2‑150 (MQ=255)
                                                                                                                                  ccaGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCa                                                                    <  2:2059983/92‑1 (MQ=255)
                                                                                                                                  ccaGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCa                                                                    >  1:2059983/1‑92 (MQ=255)
                                                                                                                                  ccaGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCGATCATGTGTGGCTCAACGCGTg          <  1:534268/150‑1 (MQ=255)
                                                                                                                                  ccaGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCGATCATGTGTGGCTCAACGCGTg          <  1:168204/150‑1 (MQ=255)
                                                                                                                                   caGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCGATCATGTGTGGCTCAACGCGTgg         <  2:437185/150‑1 (MQ=255)
                                                                                                                                   caGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCGATCATGTGTGGCTCAACGCGTgg         >  2:492054/1‑150 (MQ=255)
                                                                                                                                   caGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCGATCATGTGTGGCTCAACGCGTgg         <  1:1416814/150‑1 (MQ=255)
                                                                                                                                    aGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCGATCATGTGTGGCt                   <  1:2312776/139‑1 (MQ=255)
                                                                                                                                    aGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCGATCATGTGTGGCt                   >  2:2312776/1‑139 (MQ=255)
                                                                                                                                         ggTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCGATCATGTGTGGCTCAACGCGTGGGATAtt   >  1:396209/1‑150 (MQ=255)
                                                                                                                                          gTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCGATCATGTGTGGCTCAACGCGTGGGATATTa  <  1:249602/150‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
TGGAGCGAGGTCGCCAGCAGGTCGCCGAAATAGCCATCGCTGCACGGCGAGCTGGGCGCCACCACCAGTACGTCGCCTGGGCGGCACTGCTCCACGGCAACATGGAACATCCAGTTGTCGCCTGGGGCCACCAGCACGGTTACCGCGCTGCCGGCCACCGCCACACCTTGCTGGATAGGGCGCACGGCAGTATTGAGCAGCCCCTTGCGACCTTGGGCCTCATGCACCGTGGCGACACCCAGCCGGCCCAGGGCATCGATCATGTGTGGCTCAACGCGTGGGATATTA  >  NC_002947/2861768‑2862055

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