Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 584,262:1 (A)7→8 100% coding (582/750 nt) appY → global transcriptional activator; DLP12 prophage

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913584,2551.A93.8% 313.2 / 12.6 96coding (575/750 nt)appYglobal transcriptional activator; DLP12 prophage
Reads supporting (aligned to +/- strand):  ref base . (3/3);  new base A (47/43);  total (50/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 110 total aligned reads displayed.

GGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAA  >  NC_000913/584157‑584348
                                                                                                    |                                                                                               
ggAATTAATTTATACAAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGc                                                                                                  <  2:2133425/99‑1 (MQ=255)
  aaTTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCaaa                                                                                              >  1:2239388/1‑100 (MQ=255)
   aTTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCaaaa                                                                                             <  1:2229181/100‑1 (MQ=255)
      naTTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCaaaaaaa                                                                                         >  1:944202/2‑100 (MQ=255)
      gATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTACGGGATGAAGGAACGTCATTTAATGAAATATTGAGAGATACTAGGATGAGGTATGCaaaaaaa                                                                                         >  2:2557048/1‑100 (MQ=255)
      gATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCaaaaaaa                                                                                         >  1:676110/1‑100 (MQ=255)
      gATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCaaaaaaa                                                                                         >  2:682484/1‑100 (MQ=255)
      gATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAAAGTCATTTAATGAAATATTGAGAGATAATAGGATGAGGTATGCaaaaaaa                                                                                         >  2:2170809/1‑100 (MQ=255)
        tttATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCaaaaaaaac                                                                                       <  1:1478685/100‑3 (MQ=255)
        tttAAACGAGTGAAAATTTAATAAAAAAAAGATTAAGCGATGAATGAACTTCATTTACTGAAATATTGAGAAATACTAGGATGAGGTATGCaaaaa                                                                                            <  2:2070205/96‑1 (MQ=255)
         ttATACGATTGAAAAAAAAAAAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCaaaaa                                                                                            <  2:2005612/95‑1 (MQ=255)
         ttATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCaaaaaaaact                                                                                      <  1:489472/100‑4 (MQ=255)
         ttATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCaaaaaaaact                                                                                      <  1:1282237/100‑4 (MQ=255)
          tataCGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCaaaaaaaactc                                                                                     >  1:1569146/1‑96 (MQ=255)
            taCGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCAt                                                                                   >  1:2312526/1‑100 (MQ=255)
              cGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATaa                                                                                 >  2:764625/1‑100 (MQ=255)
              cGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATaa                                                                                 >  1:915185/1‑100 (MQ=255)
              cGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATaa                                                                                 >  1:1355657/1‑100 (MQ=255)
              cGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATaa                                                                                 >  1:133809/1‑100 (MQ=255)
                aGTAAAAGTTTAATAAAAAAAAGATTAAGGGTTGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAAc                                                                                <  1:1502319/99‑1 (MQ=255)
                 gTGAAAGTTTAATAAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGAAACAAGGATGAGGTATGCAAAAAAAACTCATAACt                                                                               >  1:1389284/1‑100 (MQ=255)
                 gTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACtt                                                                              >  2:1909048/1‑100 (MQ=255)
                 gTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTc                                                                             >  1:2042545/1‑100 (MQ=255)
                  tGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTc                                                                             <  1:463348/100‑1 (MQ=255)
                     aaGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCaaa                                                                          >  2:2531823/1‑99 (MQ=255)
                     aaGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAAc                                                                         <  1:2531838/100‑1 (MQ=255)
                     aaGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGAAACAAGGATGAGGTATGCAAAAAAAACTCATAACTTCaaa                                                                          >  2:31122/1‑100 (MQ=255)
                      aGTATAATAAAAAAAAGATTAAGGAATGAATGAACGTCATTTACTAAAATATTGAGAGATACTAGGATGAGGAATGCAAAAAAAACTCATAACTTCAAAc                                                                         <  2:352098/100‑1 (MQ=255)
                           aataaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTa                                                                    <  2:1458861/100‑1 (MQ=255)
                             taaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTa                                                                    >  1:944928/1‑98 (MQ=255)
                             taaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTAtt                                                                  <  2:2409197/100‑1 (MQ=255)
                             taaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTAtt                                                                  >  2:1718182/1‑100 (MQ=255)
                             taaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTAt                                                                   >  2:2152644/1‑99 (MQ=255)
                               aaaaaaaaGATTAAGGGATGAAGGAACGTCATTTCCTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCGAACACTTATTc                                                                 <  2:1343018/100‑1 (MQ=255)
                               aaaaaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTc                                                                 <  1:854788/100‑1 (MQ=255)
                               aaaaaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTc                                                                 <  1:1452554/100‑1 (MQ=255)
                               aaaaaaaaGATTAAGGGAAGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAGCTCATAACTTCAAACTCTTATTc                                                                 <  1:2297030/100‑1 (MQ=255)
                                  aaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAAGTATAAAAAAAAAACTCATAACTTCAAACTCTTATTCTAt                                                              <  2:944920/100‑1 (MQ=255)
                                   aaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATc                                                             >  1:149912/1‑100 (MQ=255)
                                       gATTAAGGGATGAAGGAACGTCATTTACTGCAATATTGAGAGATACTAGGATGAGGTATGCAAAAAACAATCATAACTTAAAACTCTTATTCTATCAAt                                                          >  2:1417448/1‑99 (MQ=255)
                                        aTTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGt                                                        <  1:1413505/100‑1 (MQ=255)
                                         ttAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTc                                                       >  1:494240/1‑100 (MQ=255)
                                          tAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCg                                                      <  1:98651/100‑1 (MQ=255)
                                           aaGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGt                                                     >  2:792150/1‑100 (MQ=255)
                                            aGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTa                                                    >  1:1496705/1‑100 (MQ=255)
                                              ggATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGc                                                  >  1:2299068/1‑100 (MQ=255)
                                              ggATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGc                                                  >  1:2012257/1‑100 (MQ=255)
                                               gATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGcc                                                 <  2:37847/100‑1 (MQ=255)
                                               gATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGcc                                                 >  2:609089/1‑100 (MQ=255)
                                                 tGAAGGAACGTCATTTACTGAAATATTCAGAGATACTAGGATTCGGTTTGCAAAAAAAACTCATAACTTCAAACTCTTCTTCTCTTAATGtatcatacta                                               >  2:1116416/1‑91 (MQ=255)
                                                  gAATGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTCCAAACTCTTATTCTATAAATGTCGTCGACCa                                               >  2:98376/1‑99 (MQ=255)
                                                    aGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCa                                               >  2:1226805/1‑97 (MQ=255)
                                                     ggAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGaaa                                           >  1:1536329/1‑100 (MQ=255)
                                                      gAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAt                                          >  2:1092071/1‑100 (MQ=255)
                                                      gAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtg                                         <  2:800486/100‑1 (MQ=255)
                                                       aaCGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtn                                         <  2:72219/100‑2 (MQ=255)
                                                         cGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtgtg                                       <  2:820227/100‑1 (MQ=255)
                                                          gTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTgg                                      <  2:11746/100‑1 (MQ=255)
                                                           tCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGc                                     <  1:2555926/100‑1 (MQ=255)
                                                           tCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGc                                     <  1:80502/100‑1 (MQ=255)
                                                           tCATTTACTGAAATATTGAGAGATACTAGGATGACGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGc                                     <  2:1886864/100‑1 (MQ=255)
                                                            cATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCt                                    >  1:173114/1‑100 (MQ=255)
                                                              tttACTGAAATATTGAGAGATACTAGGATAAAAAAAAAAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCtat                                  <  2:1861371/100‑1 (MQ=255)
                                                                tACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGCCTATaa                                >  1:944839/1‑100 (MQ=255)
                                                                 aCTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAAc                               <  2:1569552/100‑1 (MQ=255)
                                                                  cTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACa                              >  1:495261/1‑100 (MQ=255)
                                                                  cTGAAATATTGAGAGATACTAGGATGAGGTACGCAAAAAAAAATCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACa                              >  1:1731842/1‑100 (MQ=255)
                                                                   tGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAg                             <  2:598261/100‑1 (MQ=255)
                                                                   tGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAg                             <  1:279653/100‑1 (MQ=255)
                                                                    gAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGt                            >  2:52062/1‑100 (MQ=255)
                                                                    gAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGt                            >  1:333397/1‑100 (MQ=255)
                                                                     aaaTATTGAGAGATACTAGGATGAGGTATAAAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAAc                               <  2:1613171/96‑1 (MQ=255)
                                                                       atatTGAGAGATACTAGGATGAAGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACt                         <  2:2276136/100‑1 (MQ=255)
                                                                            gagagaTACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcata                    >  1:1052609/1‑100 (MQ=255)
                                                                            gagagaTACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcata                    >  1:477016/1‑100 (MQ=255)
                                                                            gagagaTACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcata                    >  2:1692084/1‑100 (MQ=255)
                                                                             agagaTACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatat                   <  2:823745/100‑1 (MQ=255)
                                                                              gagaTACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatt                  <  2:323408/100‑1 (MQ=255)
                                                                               agaTACTAGGAAGAAATAAGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatattt                 <  1:174501/100‑1 (MQ=255)
                                                                                gaTACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttc                >  1:541663/1‑100 (MQ=255)
                                                                                 aTACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttca               <  1:764632/100‑1 (MQ=255)
                                                                                 aTACTAGAATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGTCTATAACAGTACttcatatttca               <  1:2170824/100‑1 (MQ=255)
                                                                                  tACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcat              >  2:1162453/1‑100 (MQ=255)
                                                                                  tACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcat              <  2:431225/100‑1 (MQ=255)
                                                                                    cTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcatat            >  1:1439142/1‑100 (MQ=255)
                                                                                     tAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATAtg           >  1:842128/1‑100 (MQ=255)
                                                                                     tAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATAtg           <  1:47688/100‑1 (MQ=255)
                                                                                      aGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATAtgt          >  1:1350440/1‑100 (MQ=255)
                                                                                       ggATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATAtgtg         >  2:1652507/1‑100 (MQ=255)
                                                                                       ggATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATAATTCATAtgtg         >  2:1135870/1‑100 (MQ=255)
                                                                                         aTGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCa       <  1:2424545/100‑1 (MQ=255)
                                                                                         aTGAGGTATGCAAAAAAAACTCATAACTTCAAACGCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCa       <  2:2168895/100‑1 (MQ=255)
                                                                                         aTGAGATATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATAtgtg         <  2:1034217/98‑1 (MQ=255)
                                                                                          tGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCAt      <  1:1247859/100‑1 (MQ=255)
                                                                                           gAGGTATGCAAAAAAAACTCATAACTTCATACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCAt      <  2:2179413/99‑1 (MQ=255)
                                                                                           gAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCAtt     <  1:817541/100‑1 (MQ=255)
                                                                                             ggTATGCAAAAAAAACTCATAACTTCAAACTCATATTCTATCAATGTTGTAGACCAGAAATGTGGATAAAACAATACTTCATATTTCATAGGTGcct      >  2:2103553/1‑95 (MQ=255)
                                                                                              gTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTaa  <  2:1316461/100‑1 (MQ=255)
                                                                                              gTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTaa  <  2:661115/100‑1 (MQ=255)
                                                                                              gTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTaa  >  1:1687392/1‑100 (MQ=255)
                                                                                                    |                                                                                               
GGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATTTAA  >  NC_000913/584157‑584348

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