Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,567,414 A→G F367S (TTC→TCC)  gadC ← glutamate:gamma‑aminobutyric acid antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,567,4140AG99.2% 1130.5 / ‑2.9 267F367S (TTC→TCC) gadCglutamate:gamma‑aminobutyric acid antiporter
Reads supporting (aligned to +/- strand):  ref base A (0/2);  new base G (123/142);  total (123/144)
Fisher's exact test for biased strand distribution p-value = 5.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

Only 100 of 280 total aligned reads displayed.

TTAAGTCAGGATGTTTAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGAAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATT  >  NC_000913/1567316‑1567511
                                                                                                  |                                                                                                  
ttAAGTCAGGATGTTTAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCagggag                                                                                                  >  2:5455635/1‑101 (MQ=255)
 tAAGTCAGGATGTTTAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCagggagg                                                                                                 <  1:6037579/101‑1 (MQ=255)
 tAAGTCAGGATGTTTAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCagggagg                                                                                                 >  1:265206/1‑101 (MQ=255)
  aaGTCAGGATGTTTAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCagggagga                                                                                                >  2:6314430/1‑101 (MQ=255)
   aGTCAGGATGTTTAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGAc                                                                                               <  1:621787/101‑1 (MQ=255)
    gTCAGGATGTTTAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACa                                                                                              >  1:4812094/1‑101 (MQ=255)
    gTCAGGATGTTTAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACa                                                                                              >  2:3059572/1‑101 (MQ=255)
     tCAGGATGTTTAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACAt                                                                                             <  1:5131525/101‑1 (MQ=255)
       aGGATGTTTAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACAtgt                                                                                           >  1:5343414/1‑101 (MQ=255)
       aGGATGTTTAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACAtgt                                                                                           >  2:5397476/1‑101 (MQ=255)
        ggATGTTTAAGAACCAACACAGTGTAGCCAATAAACAGCATGAAGTGGGCACACAGAAAAACCACCACCGTCAGCGCCAATGCGATCAGGGAGGACAtgtt                                                                                          <  1:3471716/101‑1 (MQ=255)
        ggATGTTTAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACAtgtt                                                                                          <  1:1003612/101‑1 (MQ=255)
        ggATGTTTAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACAtgtt                                                                                          <  1:5455635/101‑1 (MQ=255)
        ggATGTTTAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACAtgtt                                                                                          >  2:5375157/1‑101 (MQ=255)
        ggATGTTTAAGAACCAACACAATGTAGCCAATAAACAGCATCAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACAtgtt                                                                                          <  2:7017219/101‑1 (MQ=255)
          aTGTTTAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACAtgttgt                                                                                        <  2:1760422/101‑1 (MQ=255)
           tGTTTAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATgttgtt                                                                                       <  2:1327088/101‑1 (MQ=255)
             tttAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCCCCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTAc                                                                                    <  2:752315/101‑1 (MQ=255)
             tttAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTAc                                                                                    <  1:950164/101‑1 (MQ=255)
              ttAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTAcc                                                                                   >  1:4785356/1‑101 (MQ=255)
               tAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCg                                                                                  >  1:748496/1‑101 (MQ=255)
                 aGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGcc                                                                                <  1:3863855/101‑1 (MQ=255)
                  gAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCa                                                                               <  1:2188737/101‑1 (MQ=255)
                   aaCCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCAc                                                                              <  1:299185/101‑1 (MQ=255)
                    aCCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCAcc                                                                             <  1:1729418/101‑1 (MQ=255)
                      cAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCgg                                                                           >  2:6841590/1‑101 (MQ=255)
                      cAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCgg                                                                           <  1:6402833/101‑1 (MQ=255)
                       aaCACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGt                                                                          >  2:2291609/1‑101 (MQ=255)
                       aaCAAAAGGTAGCCAATAAACAGCATGAAATAAACACACAGATAAATCGCCACCGCCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGt                                                                          <  2:6647515/101‑1 (MQ=255)
                        acacAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTTCCGCCACCGGTa                                                                         <  2:1493447/101‑1 (MQ=255)
                         cacaATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTAt                                                                        <  1:6318796/101‑1 (MQ=255)
                         cacaATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTAt                                                                        >  2:2859773/1‑101 (MQ=255)
                         cacaATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTAt                                                                        <  1:1473725/101‑1 (MQ=255)
                          acaATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTAtt                                                                       >  2:4704689/1‑101 (MQ=255)
                          acaATGTAGCCAAGAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTAtt                                                                       <  2:4144261/101‑1 (MQ=255)
                           caATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTg                                                                      <  2:2263434/101‑1 (MQ=255)
                           caATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTg                                                                      >  1:1319134/1‑101 (MQ=255)
                            aaTGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTgg                                                                     >  2:554165/1‑101 (MQ=255)
                             aTGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCCCCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGt                                                                    <  2:2752650/101‑1 (MQ=255)
                               gTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGATGACATGTTGTTACCGCCACCGGTATTGGTGa                                                                  >  1:4865563/1‑101 (MQ=255)
                                tAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAg                                                                 <  2:3128958/101‑1 (MQ=255)
                                      aTAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGgatgat                                                           <  2:3958225/101‑1 (MQ=255)
                                      aTAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGgatgat                                                           <  2:5026067/101‑1 (MQ=255)
                                        aaaCAGCATGAAATAAGCACACAGATAAATCTCCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTTCCGCCACCGGTATTGGTGAGGATGATCa                                                         <  2:1931577/101‑1 (MQ=255)
                                         aaCAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCaa                                                        <  1:2500513/101‑1 (MQ=255)
                                           cAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAAcg                                                      >  2:862224/1‑101 (MQ=255)
                                           cAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAAcg                                                      <  1:4916395/101‑1 (MQ=255)
                                             gCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAAcgcg                                                    >  1:28361/1‑101 (MQ=255)
                                               aTGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGAt                                                  <  1:357630/101‑1 (MQ=255)
                                               aTGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGAt                                                  <  1:525561/101‑1 (MQ=255)
                                                tGAAATAAGCACCCAGATAAATCACCCCCGACAGCGCCAGTGCGATCAGGGAGGACATGTTGTTCAACTCAACCCTATTGGTGAGGATGATCAACCCGATa                                                 <  2:6943053/101‑1 (MQ=255)
                                                tGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATa                                                 >  2:5170595/1‑101 (MQ=255)
                                                tGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATa                                                 <  2:105992/101‑1 (MQ=255)
                                                 gAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCTACCGGTATTGGTGAGGATGATCAACGCGATAg                                                <  2:3556486/101‑1 (MQ=255)
                                                 gAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAg                                                <  1:5679780/101‑1 (MQ=255)
                                                 gAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAg                                                <  1:6645848/101‑1 (MQ=255)
                                                  aaaTAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGCGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGa                                               >  2:5286548/1‑101 (MQ=255)
                                                   aaTAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGAc                                              >  1:4422097/1‑101 (MQ=255)
                                                     tAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGt                                            <  2:3868063/101‑1 (MQ=255)
                                                      aaGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTa                                           >  1:1636659/1‑101 (MQ=255)
                                                          acacaGATAAATCTCCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCa                                       >  2:5647794/1‑101 (MQ=255)
                                                           cacaGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCAc                                      <  1:2295114/101‑1 (MQ=255)
                                                            acaGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGCATTGGTGAGGATGATCAACGCGATAGACGTAATCAcc                                     >  1:417404/1‑101 (MQ=255)
                                                             caGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCAcaa                                    >  2:5982017/1‑99 (MQ=255)
                                                             caGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCa                                    >  1:3954256/1‑101 (MQ=255)
                                                               gATAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGc                                  <  1:3122876/101‑1 (MQ=255)
                                                                aTAAATCTCCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCt                                 <  2:319782/101‑1 (MQ=255)
                                                                aTAAATCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCt                                 >  2:4593164/1‑101 (MQ=255)
                                                                  aaaTCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGc                               >  1:17137/1‑101 (MQ=255)
                                                                  aaaTCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGc                               >  2:6045698/1‑101 (MQ=255)
                                                                   aaTCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCg                              >  1:6480252/1‑101 (MQ=255)
                                                                   aaTCACCACCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCg                              >  1:4603327/1‑101 (MQ=255)
                                                                      caccacCGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGaaa                           >  1:5388667/1‑101 (MQ=255)
                                                                       accaccGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAAt                          >  2:6159319/1‑101 (MQ=255)
                                                                         caccGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGa                        >  1:2091048/1‑101 (MQ=255)
                                                                         caccGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGa                        >  2:5087361/1‑101 (MQ=255)
                                                                           ccGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGAcc                      <  2:2393481/101‑1 (MQ=255)
                                                                            cGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCGCCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCa                     <  1:2365798/101‑1 (MQ=255)
                                                                            cGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCa                     >  1:3417159/1‑101 (MQ=255)
                                                                            cGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCa                     >  1:3694196/1‑101 (MQ=255)
                                                                            cGTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCa                     <  1:5569646/101‑1 (MQ=255)
                                                                             gTCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAg                    >  2:4111973/1‑101 (MQ=255)
                                                                              tCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGc                   <  2:2655246/101‑1 (MQ=255)
                                                                              tCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGc                   <  1:6898204/101‑1 (MQ=255)
                                                                              tCAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGc                   <  2:1369234/101‑1 (MQ=255)
                                                                               cAGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCg                  <  2:4434178/101‑1 (MQ=255)
                                                                                aGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGt                 >  1:4243926/1‑101 (MQ=255)
                                                                                aGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGt                 <  2:3611174/101‑1 (MQ=255)
                                                                                aGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGt                 >  1:5922474/1‑101 (MQ=255)
                                                                                aGCGCCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGt                 <  2:2642475/101‑1 (MQ=255)
                                                                                  cgcCAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGAATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTa               >  1:4051761/1‑101 (MQ=255)
                                                                                    ccAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTAcc             >  1:2692656/1‑101 (MQ=255)
                                                                                    ccAGTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTAcc             >  1:833858/1‑101 (MQ=255)
                                                                                       gTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGt          >  1:4213413/1‑101 (MQ=255)
                                                                                       gTGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGt          <  2:6486222/101‑1 (MQ=255)
                                                                                        tGCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTa         <  1:1096968/101‑1 (MQ=255)
                                                                                         gCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTAc        <  1:5693802/101‑1 (MQ=255)
                                                                                         gCGATCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTAc        <  2:1233915/101‑1 (MQ=255)
                                                                                            aTCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGcc     <  1:2247381/101‑1 (MQ=255)
                                                                                            aTCAGGGAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGcc     >  2:3747659/1‑101 (MQ=255)
                                                                                                  |                                                                                                  
TTAAGTCAGGATGTTTAAGAACCAACACAATGTAGCCAATAAACAGCATGAAATAAGCACACAGATAAATCACCACCGTCAGCGCCAGTGCGATCAGGAAGGACATGTTGTTACCGCCACCGGTATTGGTGAGGATGATCAACGCGATAGACGTAATCACCAGCTGCGAAATGACCAGCGTTACCGGTACGCCATT  >  NC_000913/1567316‑1567511

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 17 ≤ ATCG/ATCG < 26 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: