Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,520,083 A→C 19.0% intergenic (+19/‑179) curA → / → mcbR curcumin/dihydrocurcumin reductase, NADPH‑dependent/colanic acid and biofilm gene transcriptional regulator, MqsR‑controlled

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,520,0830AC19.0% 491.1 / 119.2 221intergenic (+19/‑179)curA/mcbRcurcumin/dihydrocurcumin reductase, NADPH‑dependent/colanic acid and biofilm gene transcriptional regulator, MqsR‑controlled
Reads supporting (aligned to +/- strand):  ref base A (80/98);  new base C (19/23);  total (100/121)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.94e-01

Only 100 of 222 total aligned reads displayed.

ATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTAT  >  NC_000913/1519985‑1520181
                                                                                                  |                                                                                                  
 tGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgc                                                                                                  <  2:4189885/100‑1 (MQ=255)
   cgcCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccg                                                                                                >  2:1293318/1‑100 (MQ=255)
    gcCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgc                                                                                               <  1:4317624/100‑1 (MQ=255)
    gcCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgc                                                                                               <  1:2997821/100‑1 (MQ=255)
     ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgcc                                                                                              >  1:584361/1‑100 (MQ=255)
      cacaGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCa                                                                                             >  2:2999449/1‑100 (MQ=255)
       acaGACGTTTATCGGCCTGCTGAAGTGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCa                                                                                             >  1:1898284/1‑99 (MQ=255)
       acaGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgt                                                                                            <  1:4185423/100‑3 (MQ=255)
        caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                           <  2:4259170/100‑4 (MQ=255)
         aGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgttg                                                                                          <  2:3976190/100‑5 (MQ=255)
          gACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGTGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATaaa                                                                                         <  2:1399302/100‑1 (MQ=255)
          gACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATaaa                                                                                         <  2:356749/100‑1 (MQ=255)
          gACGTTTATCGACCTGGTGACTGGTAAAAACTGCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                           <  1:112988/98‑4 (MQ=255)
            cGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATaaaaa                                                                                       >  1:758317/1‑100 (MQ=255)
            cGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATaaaaa                                                                                       >  2:4396344/1‑100 (MQ=255)
            cGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATaaaaa                                                                                       >  2:708257/1‑100 (MQ=255)
            cGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATaaaaa                                                                                       <  2:2520445/100‑1 (MQ=255)
             gTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAATGGCGGCTTACgccgccgttgaaac                                                                                      <  1:2071584/100‑9 (MQ=38)
             gTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgttgaaac                                                                                      <  2:547371/100‑9 (MQ=255)
              tttATCGGCTTGCTGAAGGGTAAAAACGTCGGCAAAGTGGTGATCGGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATc                                                                                     <  1:780700/100‑1 (MQ=255)
              tttATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATc                                                                                     >  2:4121807/1‑100 (MQ=255)
                tATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCa                                                                                   >  1:3069607/1‑100 (MQ=255)
                 aTCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAg                                                                                  <  2:4230464/100‑1 (MQ=255)
                 aTCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAg                                                                                  <  1:3448441/100‑1 (MQ=255)
                 aTCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAg                                                                                  <  1:2297857/100‑1 (MQ=255)
                   cGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGc                                                                                <  2:3054317/100‑1 (MQ=255)
                    ggCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCa                                                                               <  2:343846/100‑1 (MQ=255)
                    ggCCTGCTGAAGGATAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTGAGACGCCATAAAAATCCAGGCa                                                                               >  1:3084422/1‑100 (MQ=255)
                     gCCTGCTGAAGGGTAAAAACTTCGGAAAAGTGGTGATCCGCGTGGCGGGTGATGCTTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCaa                                                                              >  1:2661729/1‑100 (MQ=255)
                         gctgAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATGAAAGATTCAACGGCGGCGTAAGCCGCCATAAAAATCGAGGCAACAAg                                                                          >  1:986137/1‑100 (MQ=255)
                          ctgAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGc                                                                         >  1:1386317/1‑100 (MQ=255)
                            gAAGGGTAAAAACTTCGGCAAAGGGGTGCTCTCCGTGGCGTGTGCTGATTCAAGCTTCAACGGCAGCGTAACCCGCCATAAAAATCCAGCCACAAAGCt                                                                        >  1:1255623/1‑99 (MQ=255)
                                ggTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGATTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAttt                                                                   >  1:919062/1‑100 (MQ=255)
                                 gTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAtttt                                                                  >  2:166646/1‑100 (MQ=255)
                                  tAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTc                                                                 <  2:4301466/100‑1 (MQ=255)
                                   aaaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCa                                                                <  2:1869607/100‑1 (MQ=255)
                                    aaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAg                                                               >  2:1564531/1‑100 (MQ=255)
                                     aaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGc                                                              <  2:1359726/100‑1 (MQ=255)
                                        cTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCttt                                                           >  2:3717400/1‑100 (MQ=255)
                                        cTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCttt                                                           <  2:453381/100‑1 (MQ=255)
                                          tCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaa                                                         >  2:3493079/1‑100 (MQ=255)
                                           cGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaat                                                        <  1:3578398/100‑1 (MQ=255)
                                             gCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaatta                                                      <  1:3302635/100‑1 (MQ=255)
                                             gCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTAAACGGCGGAGTAAGCCGCCATATAAATGCGGGAGAAAAGGTTATGTTCAgtg                                                             >  2:2868363/1‑91 (MQ=255)
                                              cAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaattaa                                                     >  1:2624700/1‑100 (MQ=255)
                                              cAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaattaa                                                     <  2:2773862/100‑1 (MQ=255)
                                              cAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaatta                                                      <  1:3314803/99‑1 (MQ=255)
                                               aaaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAAAAAAAAAAAAAAGCTTATTTTCAGCTttaatta                                                      <  1:131893/98‑1 (MQ=255)
                                                aaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaact                                                   >  1:1543022/1‑100 (MQ=255)
                                                  gtggtgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaactaa                                                 >  2:907153/1‑100 (MQ=255)
                                                   tggtgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaactaac                                                >  2:1735685/1‑100 (MQ=255)
                                                    ggtgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaactaac                                                <  1:4277862/99‑1 (MQ=255)
                                                    ggtgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAAcc                                               >  2:4018473/1‑100 (MQ=255)
                                                    ggtgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAAcc                                               >  1:2244977/1‑100 (MQ=255)
                                                    ggtgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAAcc                                               <  2:2741763/100‑1 (MQ=255)
                                                    ggtgATCCGCGTGGCGGGTGATGATTAAAGTTTAAACGGCGGCGTAAGCCGCCATAAAAATCCAAGCAAAAAGCTTATTTTCAGCTTTAATTAACTAAcc                                               <  1:360949/100‑1 (MQ=255)
                                                       gATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAACAAGCTTATTTTCAGCTTTAATTAACTAACcaca                                            >  1:3040460/1‑100 (MQ=255)
                                                         tCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCAcatc                                          >  2:2844664/1‑100 (MQ=255)
                                                         tCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCAcatc                                          >  2:2117898/1‑100 (MQ=255)
                                                          ccGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCAcatca                                         <  1:4066828/100‑1 (MQ=255)
                                                             cgTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTg                                      <  2:1144556/100‑1 (MQ=255)
                                                               tGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGtt                                    >  2:227013/1‑100 (MQ=255)
                                                               tGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGtt                                    >  2:3853147/1‑100 (MQ=255)
                                                                ggCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTc                                   <  1:722606/100‑1 (MQ=255)
                                                                 gCGGGTGATGACTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTct                                  <  2:2715244/100‑1 (MQ=255)
                                                                   gggTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctgc                                <  2:1200202/100‑1 (MQ=255)
                                                                      tgatgaTTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAAACAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctgct                               <  1:2528554/98‑1 (MQ=255)
                                                                       gatgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTgctgct                            <  1:2795722/100‑1 (MQ=255)
                                                                        atgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCta                           <  2:3323449/100‑1 (MQ=255)
                                                                         tgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTAAGCTTTAATTAACTAACCACATCATTGTTctgctg                              >  2:895539/1‑96 (MQ=255)
                                                                          gatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaat                         <  2:3923578/100‑1 (MQ=255)
                                                                          gatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaat                         <  2:588590/100‑1 (MQ=255)
                                                                           atTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaata                        >  1:1997707/1‑100 (MQ=255)
                                                                            ttAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaataa                       >  2:1124213/1‑100 (MQ=255)
                                                                             tAAAGTTTCAACGGCGGCGTAAGCCGCCATAAATATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATaac                      <  1:3728690/100‑1 (MQ=255)
                                                                             tAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCGTTAATTAACTAACCACATCATTGTTCTGCTGCtaataa                       >  1:628033/1‑99 (MQ=255)
                                                                               aaGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATaacaa                    <  2:3347068/100‑1 (MQ=255)
                                                                                gatttttatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGTTAATaacaac                   >  2:1957248/10‑100 (MQ=255)
                                                                                 gTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAacaaca                  >  2:827370/1‑100 (MQ=255)
                                                                                 gTTTCAACGGCGACGTAAGCCGCCATAAAAATCCAAACAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTGCTGCTGCTAATaacaac                   <  1:2117891/99‑1 (MQ=255)
                                                                                  tttCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAAcaacaa                 >  2:1984406/1‑100 (MQ=255)
                                                                                   ttttatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACaaa                >  1:1596567/7‑100 (MQ=255)
                                                                                    tttatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCCCATCATTGTTCTGCTGCTAATAacaaca                  >  1:3570387/6‑97 (MQ=255)
                                                                                     cAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAgg              <  1:1269992/100‑1 (MQ=255)
                                                                                     cAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAgg              >  1:4348871/1‑100 (MQ=255)
                                                                                      tatggcggcTTACGCCGCCATAAAAATCCTGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTTCTGCTAATAACAACAAAgg              >  1:1021494/4‑99 (MQ=255)
                                                                                      tatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTCTTTTCCGCCTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAgg              >  1:244218/4‑99 (MQ=255)
                                                                                      tatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAgg              >  1:2978286/4‑99 (MQ=255)
                                                                                      tatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAgg              >  1:1443145/4‑99 (MQ=255)
                                                                                      aaCGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGt             <  2:3201567/100‑1 (MQ=255)
                                                                                       atggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGt             >  1:3417070/3‑99 (MQ=255)
                                                                                       atggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTg            >  1:1786868/3‑100 (MQ=255)
                                                                                       aCGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACGAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTg            <  1:2779570/100‑1 (MQ=255)
                                                                                         ggcggcGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg          <  2:649275/100‑1 (MQ=255)
                                                                                            ggcgTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTgg       >  2:1542568/1‑100 (MQ=255)
                                                                                             gcgTAAGCCGCCATAAACATCCAGGAAAAAAGCTTATATTCAGCTATAATTAACTAACCACAACATTATTCTGCTGCTAATAACAAAAAAGGTGAGTgg       >  1:2357421/1‑99 (MQ=255)
                                                                                             gcgTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGt      >  2:3263712/1‑100 (MQ=255)
                                                                                             gcgTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGt      >  2:3131393/1‑100 (MQ=255)
                                                                                              cgTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGtt     <  2:1605964/100‑1 (MQ=255)
                                                                                              cgTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGtt     >  2:1312864/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
ATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTAT  >  NC_000913/1519985‑1520181

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: