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.A96.1% 335.5 / 6.8 102coding (575/750 nt)appYglobal transcriptional activator; DLP12 prophage
Reads supporting (aligned to +/- strand):  ref base . (0/4);  new base A (46/52);  total (46/56)
Fisher's exact test for biased strand distribution p-value = 1.25e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.50e-01

Only 100 of 112 total aligned reads displayed.

CGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGC‑‑AAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATT  >  NC_000913/584156‑584345
                                                                                                    |                                                                                             
cGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGc                                                                                                <  1:1775037/100‑1 (MQ=255)
 ggAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGC‑‑a                                                                                             >  2:271435/1‑100 (MQ=255)
  gAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGC‑‑aa                                                                                            >  2:1812489/1‑100 (MQ=255)
       gATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGC‑‑aaaaaaa                                                                                      >  1:2786229/1‑100 (MQ=255)
           tataTGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGGCAGTTACTGAAATATTGAGGGATACTAGGATGGTGTATGa‑‑aaaaaaaactc                                                                                  <  2:1490817/100‑5 (MQ=255)
           tataCGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGC‑‑aaaaaaaactc                                                                                  <  1:3041838/100‑5 (MQ=255)
           tataCGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGC‑‑aaaaaaaaact                                                                                  >  2:3025755/1‑96 (MQ=255)
           tataCGAGTAAAAGTATAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGC‑‑aaaaaaaactc                                                                                  <  2:2295796/100‑5 (MQ=255)
            ataCGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGC‑‑aaaaaaaactca                                                                                 <  1:2402717/100‑6 (MQ=255)
             taCGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTAAGAGATACTAGGATGAGGTATGCAAAAAAAACTCAt                                                                                >  1:770420/1‑100 (MQ=255)
               cGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATaa                                                                              >  1:2970319/1‑100 (MQ=255)
                 aGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAAACTCATAAc                                                                             <  1:3025788/100‑1 (MQ=255)
                 aGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACt                                                                            >  2:1726571/1‑100 (MQ=255)
                 aGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAAGTATGCAAAAAAAACTCATAAc                                                                             <  2:389315/99‑1 (MQ=255)
                   tGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTc                                                                          >  2:3044259/1‑100 (MQ=255)
                    gAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCa                                                                         >  1:2933567/1‑100 (MQ=255)
                    gAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAAATCATAACTTCa                                                                         >  2:813064/1‑100 (MQ=255)
                     aaaGTGTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCa                                                                         <  2:1286744/99‑1 (MQ=255)
                       tgTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAAc                                                                      <  1:107298/99‑1 (MQ=255)
                       aGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTAAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAAc                                                                      >  1:833982/1‑100 (MQ=255)
                        gTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCaaa                                                                       >  2:107297/1‑98 (MQ=255)
                        gTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAAct                                                                     <  1:401070/100‑1 (MQ=255)
                           taataaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCtt                                                                  <  2:1103301/100‑1 (MQ=255)
                             ataaAAAAAAGATTCAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGC‑‑AAAAAATACTCATAACTTCAAACTCTTAt                                                                <  2:1808388/100‑1 (MQ=255)
                             ataaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTAt                                                                <  1:2514239/100‑1 (MQ=255)
                             ataaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTAt                                                                >  1:1170809/1‑100 (MQ=255)
                              taaAAAAAAGATTAAGGGATGAAGGAACGTTTTTTATTGAAATATTGAGTGATTCTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTAtt                                                               <  1:892130/100‑1 (MQ=255)
                              taaAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTAtt                                                               <  1:2432202/100‑1 (MQ=255)
                               aaaaaaaaGATTAATGGATGAATGAACGTCATTTTCTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTc                                                              <  1:312022/100‑1 (MQ=255)
                               aaaaaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTc                                                              <  1:1726591/100‑1 (MQ=255)
                               aaaaaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTc                                                              >  2:2259641/1‑100 (MQ=255)
                               aaaaaaaaGATTAAGGGATAAAGGAACGTCATTTACTGAAATATTGAGAGATATTAGGATGATGTATGCAAAAAAAACTCATAACTTCAAACTCTTAtt                                                               >  2:1806082/1‑99 (MQ=255)
                                 aaaaaaGATTAAGGGATGAAGGCACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCt                                                             >  1:2483323/1‑99 (MQ=255)
                                 aaaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTAATAACTTCAAACTCTTATTCTa                                                            >  1:2090658/1‑100 (MQ=255)
                                  aaaaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGTATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTAt                                                           <  1:18437/100‑1 (MQ=255)
                                   aaaaaaCTAATGGATTAAGGAACGTCATTTACTGAAATATTGAGTGTTATTATTAAAAAATAAAAAAAAAAAACTCATAACTTCAAACTCTTATTCTATc                                                          <  2:844962/100‑1 (MQ=255)
                                    aaaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCa                                                         <  2:1472733/100‑1 (MQ=255)
                                     aaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCaa                                                        <  2:803819/100‑1 (MQ=255)
                                     aaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCaa                                                        <  2:803912/100‑1 (MQ=255)
                                     aaGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCaa                                                        <  1:2876934/100‑1 (MQ=255)
                                      aGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAAt                                                       <  2:2152169/100‑1 (MQ=255)
                                       gATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATg                                                      >  2:794991/1‑100 (MQ=255)
                                        aTTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGt                                                     >  2:449921/1‑100 (MQ=255)
                                         ttAAGGGATGAATGAACGTTATTTACTGAAATATTGAGAGATAATAGAATGAGGTATGCAAAAAAAACTCATAACTTCATACTCTTATTCTATCAATGTc                                                    <  2:1495254/100‑1 (MQ=255)
                                           aaGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGt                                                  >  2:593760/1‑100 (MQ=255)
                                           aaGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGt                                                  <  2:2480723/100‑1 (MQ=255)
                                            aGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGTATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGt                                                  <  2:89081/99‑1 (MQ=255)
                                            aGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTa                                                 <  2:56971/100‑1 (MQ=255)
                                             gggATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAg                                                >  2:492488/1‑100 (MQ=255)
                                              ggATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGc                                               <  2:2914652/100‑1 (MQ=255)
                                              ggATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGc                                               >  2:482668/1‑100 (MQ=255)
                                                 tGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGt                                                  >  2:1333446/1‑94 (MQ=255)
                                                     ggAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGaaa                                        >  1:2436836/1‑100 (MQ=255)
                                                       aaCGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtg                                      >  2:786199/1‑100 (MQ=255)
                                                        aCGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtg                                      >  1:1983428/1‑99 (MQ=255)
                                                         cGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtgtg                                    <  2:2506363/100‑1 (MQ=255)
                                                         cGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAAtgtg                                    >  2:1675928/1‑100 (MQ=255)
                                                          gTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGC‑‑AAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGc                                  <  2:2188867/100‑1 (MQ=255)
                                                           tCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTgg                                   <  2:2786197/99‑1 (MQ=255)
                                                           tCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGc                                  <  1:215583/100‑1 (MQ=255)
                                                           tCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGc                                  <  1:1203642/100‑1 (MQ=255)
                                                            cATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCt                                 >  2:509434/1‑100 (MQ=255)
                                                             aTTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCta                                <  1:2178306/100‑1 (MQ=255)
                                                             aTTTACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCta                                <  1:361600/100‑1 (MQ=255)
                                                              tttACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCtat                               <  2:1679526/100‑1 (MQ=255)
                                                              tttACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCtat                               <  2:1106858/100‑1 (MQ=255)
                                                                tACTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATaa                             <  2:225758/100‑1 (MQ=255)
                                                                 aCTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAAc                            >  2:1360287/1‑100 (MQ=255)
                                                                  cTGAAATGTTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACa                           >  1:1587466/1‑100 (MQ=255)
                                                                  cTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACa                           >  2:2237021/1‑100 (MQ=255)
                                                                  cTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGACTATAACa                           >  1:3037472/1‑100 (MQ=255)
                                                                  cTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCAAGAAATGTGGCTATAACa                           >  2:436504/1‑100 (MQ=255)
                                                                  cTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCATAACCAAGAAATTTGACTATAACa                           >  1:2530922/1‑100 (MQ=255)
                                                                  cTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTAATTCTATCAATGTCGTAGCCCAGAAATATGGCTATaata                           >  1:1007068/1‑98 (MQ=255)
                                                                  cTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAAATCTTATTATATCAATATCGTAAAACAGAAATGTGGCTATAACa                           >  1:810919/1‑100 (MQ=255)
                                                                  cTGAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTAATAACTTCAAACTTTTATTCTATCAATGTAGTAGCCAAGAAATGTGGTTATAACa                           >  2:1408365/1‑100 (MQ=255)
                                                                    gAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCATAAATGAGGCTATAACAGt                         >  1:1867629/1‑100 (MQ=255)
                                                                    gAAATATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAAAGTCGTAGCCAAGAAATGTGGCTATAACAGt                         >  2:2592557/1‑100 (MQ=255)
                                                                     aaaTATTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAACCCAGAAATGTGGCTATAACAGTa                        >  1:1725693/1‑100 (MQ=255)
                                                                        tatTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACtt                     <  1:2195404/100‑1 (MQ=255)
                                                                         atTGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttc                    <  2:3015258/100‑1 (MQ=255)
                                                                           tGAGAGATACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcat                  >  1:2580769/1‑100 (MQ=255)
                                                                            gagagaTACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcata                 <  2:770416/100‑1 (MQ=255)
                                                                            gagagaTACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcata                 >  1:2126111/1‑100 (MQ=255)
                                                                               agaTACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatattt              <  1:2650447/100‑1 (MQ=255)
                                                                                 aTACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttca            <  1:940243/100‑1 (MQ=255)
                                                                                  tACTAGGATGAGGTATGCAAAAAAAACTCATATCTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcat           <  1:1478900/100‑1 (MQ=255)
                                                                                  tACTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcat           <  2:1278232/100‑1 (MQ=255)
                                                                                  tACTAGGATGAGGTATGC‑‑AAAAAAATCTCATAACTTCAATCTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACttcatatttcat           <  1:34632/100‑1 (MQ=255)
                                                                                   aCTAGGATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGACTATAACAGTACttcatatttcata          >  1:1205669/1‑100 (MQ=255)
                                                                                       ggATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATAtgtg      >  2:2211921/1‑100 (MQ=255)
                                                                                       ggATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATAtgtg      >  1:1135159/1‑100 (MQ=255)
                                                                                       ggATGAAGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATAtg        <  2:1867607/98‑1 (MQ=255)
                                                                                        gATGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGc     <  2:1468867/100‑1 (MQ=255)
                                                                                         atggagTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCa    <  2:2980646/95‑1 (MQ=255)
                                                                                         aTGAGGTATGCAAAAAAAACTCATAACTTCAATCTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCa    <  1:122814/100‑1 (MQ=255)
                                                                                         aTGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCa    <  2:1711791/100‑1 (MQ=255)
                                                                                         aTGAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATATCAAAACCCAAAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCa    <  1:2529868/100‑1 (MQ=255)
                                                                                          tGAGGTATGC‑‑AAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCAtt  <  1:2378731/100‑1 (MQ=255)
                                                                                           gAGGTATGCAAAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCAtt  <  1:3032250/100‑1 (MQ=255)
                                                                                                    |                                                                                             
CGGAATTGATTTATACGAGTGAAAGTTTAATAAAAAAAAGATTAAGGGATGAAGGAACGTCATTTACTGAAATATTGAGAGATACTAGGATGAGGTATGC‑‑AAAAAAACTCATAACTTCAAACTCTTATTCTATCAATGTCGTAGCCCAGAAATGTGGCTATAACAGTACTTCATATTTCATATGTGCATT  >  NC_000913/584156‑584345

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