Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,428,271 A→C 8.4% S152S (TCT→TCG insH1 ← IS5 transposase and trans‑activator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,428,2710AC 8.4% 154.9 / 7.0 60S152S (TCT→TCGinsH1IS5 transposase and trans‑activator
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base C (3/2);  total (31/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

CTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGA  >  NC_000913/1428127‑1428413
                                                                                                                                                |                                                                                                                                              
cTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTagag                                                                                                                                             <  5:39599/148‑1 (MQ=32)
    ttGGCCGCGGTGGTGACCAGGCTGTGGGTCATGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTgg                                                                                                                                         >  6:69520/1‑148 (MQ=17)
    tgggCCGCGGTGGTGACCTGGCTGTGGGTCAGGCCACTCTGGGCATCGAGACCAATGTGGGCCTTCATGCCTAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCTTCTCCGGCTAGCGTTGCGGCTCTTTGTTCTTGTTAGAGCTgg                                                                                                                                         <  2:88618/146‑1 (MQ=11)
                   aCCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGAtgg                                                                                                                          <  7:301156/148‑1 (MQ=32)
                     cAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGACTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATGGCGGTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGAtggtg                                                                                                                        <  2:191557/148‑1 (MQ=12)
                      aGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGAtggtgg                                                                                                                       >  6:205468/1‑148 (MQ=32)
                       ggCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGAtggtgg                                                                                                                       >  7:119260/1‑147 (MQ=32)
                          tgtgGGTCCGGCCACTCTTGGCATCGACACCAATGCGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTGGTTCTTGGTAGAGCTGGGTGCCTCAatga                                                                                                                             <  2:82209/138‑1 (MQ=12)
                           gtgGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCCTGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAatgat                                                                                                                            >  1:82209/1‑138 (MQ=17)
                                tCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCaa                                                                                                             <  5:205468/148‑1 (MQ=34)
                                    gCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTg                                                                                                         <  1:130606/148‑1 (MQ=37)
                                     ccACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCGGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGttc                                                                                                        >  5:313591/1‑146 (MQ=37)
                                       aCTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCCCCACAGTGCCt                                                                                                      >  2:249731/1‑148 (MQ=255)
                                          ctTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTg                                                                                                                                          <  2:293154/109‑1 (MQ=32)
                                          ctTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTg                                                                                                                                          >  1:293154/1‑109 (MQ=32)
                                           ttGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTg                                                                                                    <  2:392326/146‑1 (MQ=37)
                                            tGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGt                                                                                                 >  6:30384/1‑148 (MQ=39)
                                              gCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGtc                                                                                                <  2:202008/147‑1 (MQ=39)
                                               cATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGtcat                                                                                              <  5:69520/148‑1 (MQ=39)
                                                  cGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGTTTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGTTGGGTGCCTCAATGATGGTGGCATCCACCaa                                                                                                             <  7:132576/130‑1 (MQ=25)
                                                  cGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGTTTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGTTGGGTGCCTCAATGATGGTGGCATCCACCaa                                                                                                             >  8:132576/1‑130 (MQ=25)
                                                  cGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTCGAGCTGGGTGCCTCAATGATGGTGGCATCCACCaaa                                                                                                            >  6:96967/1‑131 (MQ=11)
                                                  cGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTCGAGCTGGGTGCCTCAATGATGGTGGCATCCACCaaa                                                                                                            <  5:96967/131‑1 (MQ=11)
                                                  cGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTcatcat                                                                                           >  8:10314/1‑148 (MQ=39)
                                                   gACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATg                                                                                          <  5:30384/148‑1 (MQ=39)
                                                     cacCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGAc                                                                                        >  2:291131/1‑148 (MQ=39)
                                                                 ccTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGgc                                                                            >  1:380868/1‑148 (MQ=39)
                                                                  cTTCATTCCAACCTGCCACTGATTGCCGTTCTTGGTCTGATTCATCTCCTGATCTCGTTGCTGCTCTTTGTTCTTGGTAGTGCTGGGTGCCTCATTGATGGTGGCATCCACCAAAGTGCCTTTGGTCATCATGACGTCTGCTTCGgcc                                                                           <  4:294935/148‑1 (MQ=14)
                                                                    tCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCaaa                                                                                                            <  6:210414/113‑1 (MQ=37)
                                                                    tCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCaaa                                                                                                            >  5:210414/1‑113 (MQ=37)
                                                                               tGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGa                                                              >  2:388115/1‑148 (MQ=39)
                                                                                   aCTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGt                                                          >  6:317052/1‑148 (MQ=39)
                                                                                     tGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGACTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCt                                                        >  3:177492/1‑148 (MQ=255)
                                                                                      gATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCtt                                                       >  2:291048/1‑148 (MQ=39)
                                                                                           ccTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAAc                                                   <  8:339454/147‑1 (MQ=39)
                                                                                                cTTTGTCTGATGATTCTCAGGATCGTGTTGCTGCTCTTTGTTCTTGGTTGTTCTGGGTGGCTCAATGATGTTGGCTGCGTCCAAAGTGGGTTGGGTCCTCATGACGCCTGCTGCGGCCAGCCGGCGATTGATGGTCTTGAACAATTg                                              <  6:313591/147‑1 (MQ=2)
                                                                                                cTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTg                                              <  7:10314/147‑1 (MQ=39)
                                                                                                   ggTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACg                                           <  8:119260/147‑1 (MQ=255)
                                                                                                    gTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACggg                                         >  5:273135/1‑148 (MQ=255)
                                                                                                        gATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTg                                                                                                    <  3:256263/85‑1 (MQ=255)
                                                                                                        gATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTg                                                                                                    >  4:256263/1‑85 (MQ=255)
                                                                                                        gATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTg                                              <  1:290989/139‑1 (MQ=40)
                                                                                                        gATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTg                                              >  2:290989/1‑139 (MQ=40)
                                                                                                        gATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGCCGGGCCAg                                     >  1:130305/1‑148 (MQ=255)
                                                                                                               ctcCGGATCGCGTTGCTGCTCTTTGTTCTTGGTCGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGAtgc                              >  4:107254/1‑148 (MQ=37)
                                                                                                               ctcCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGAtgc                              >  6:269352/1‑148 (MQ=255)
                                                                                                               ctcCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGAtgc                              >  4:210359/1‑148 (MQ=255)
                                                                                                                  cGGATCGCGTTGCTGCTCTTTGTTCTTGGTCGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGAtgctgc                           <  3:107254/148‑1 (MQ=37)
                                                                                                                  cGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGAtgct                             <  3:210359/146‑1 (MQ=255)
                                                                                                                  cGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGAtgc                              <  5:269352/145‑1 (MQ=255)
                                                                                                                         cgTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCa                    <  1:291131/148‑1 (MQ=255)
                                                                                                                           ttGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAgg                  >  5:415068/1‑148 (MQ=255)
                                                                                                                                  tctTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCgg            >  6:19860/1‑147 (MQ=255)
                                                                                                                                      tGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGaaa         <  7:294247/146‑1 (MQ=255)
                                                                                                                                      tGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAAtt       >  6:362080/1‑148 (MQ=255)
                                                                                                                                        ttCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTc      <  6:415068/147‑1 (MQ=255)
                                                                                                                                          cTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCat    >  5:179345/1‑147 (MQ=255)
                                                                                                                                          cTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCat    >  6:126727/1‑147 (MQ=255)
                                                                                                                                           ttGGTCGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATAGAAGGGCCAGTTGATGCTGCTCGCGCAGGTGGCGGAAATTCatga  >  1:301282/1‑148 (MQ=255)
                                                                                                                                             ggTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCatg   <  5:431533/145‑1 (MQ=255)
                                                                                                                                             ggTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCatg   <  5:362080/145‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                              
CTCGTTGGCCGCGGTGGTGACCAGGCTGTGGGTCAGGCCACTCTTGGCATCGACACCAATGTGGGCCTTCATGCCAAAGTGCCACTGATTGCCTTTCTTGGTCTGATGCATCTCCGGATCGCGTTGCTGCTCTTTGTTCTTGGTAGAGCTGGGTGCCTCAATGATGGTGGCATCCACCAAAGTGCCTTGGGTCATCATGACGCCTGCTTCGGCCAGCCAGCGATTGATGGTCTTGAACAATTGACGGGCCAGTTGATGCTGCTCGAGCAGGTGGCGGAAATTCATGA  >  NC_000913/1428127‑1428413

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