Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,520,080 G→T 16.1% intergenic (+16/‑182) 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,0800GT16.1% 499.4 / 75.7 186intergenic (+16/‑182)curA/mcbRcurcumin/dihydrocurcumin reductase, NADPH‑dependent/colanic acid and biofilm gene transcriptional regulator, MqsR‑controlled
Reads supporting (aligned to +/- strand):  ref base G (81/74);  new base T (15/15);  total (96/90)
Fisher's exact test for biased strand distribution p-value = 8.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 194 total aligned reads displayed.

GAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTT  >  NC_000913/1519981‑1520178
                                                                                                   |                                                                                                  
gagaATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGgcggcg                                                                                                    >  1:708053/1‑100 (MQ=255)
  gaATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTa                                                                                                  <  1:156734/100‑1 (MQ=255)
  gaAAGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTa                                                                                                  <  1:1962715/100‑1 (MQ=255)
       cgcCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccg                                                                                             >  1:1328894/1‑100 (MQ=255)
         ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgcc                                                                                           >  2:3330566/1‑100 (MQ=255)
         ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgcc                                                                                           >  1:3056972/1‑100 (MQ=255)
         ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgcc                                                                                           >  1:3234320/1‑100 (MQ=255)
         ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgcc                                                                                           >  1:578421/1‑100 (MQ=255)
           acaGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgt                                                                                         <  2:2461204/100‑3 (MQ=255)
           acaGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgt                                                                                         <  1:1844550/100‑3 (MQ=255)
            cgggcGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCTAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                        <  2:2561429/96‑4 (MQ=255)
            caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                        <  1:637411/100‑4 (MQ=255)
            caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                        <  1:2084707/100‑4 (MQ=255)
            caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                        <  2:1359054/100‑4 (MQ=255)
            caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                        <  2:586764/100‑4 (MQ=255)
            caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                        <  1:2648232/100‑4 (MQ=255)
            caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATa                                                                                        >  1:3070991/1‑100 (MQ=255)
             aGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgttg                                                                                       <  1:2283888/100‑5 (MQ=255)
             aGACGTTTATCGGCCGTCGGAATGGTGAAGACTTCGGCAAAGTGGGGAGCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgt                                                                                         <  2:2952155/98‑3 (MQ=255)
                cGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATaaaaa                                                                                    <  1:3395255/100‑1 (MQ=255)
                cGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGGGATCCGCGTGGCGGTTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATaaa                                                                                      <  2:2832369/98‑1 (MQ=255)
                 gTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAAt                                                                                   <  1:84539/100‑1 (MQ=255)
                   ttATCGGCCTGCTGAAGGGTAAAAACTTCTGCAAAGTGGTGATCCTCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCCTAAGCCGCCATAAAAATcc                                                                                 <  1:386464/100‑1 (MQ=255)
                   ttATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATcc                                                                                 >  2:3363108/1‑100 (MQ=255)
                     aTCGGCCTGCTGAAGGGTATAAACTTCGGCAAAGTGGTGAGCCGCGTGGCGTGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAg                                                                               <  2:455305/100‑1 (MQ=255)
                     aTCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAg                                                                               <  2:1799499/100‑1 (MQ=255)
                       cGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGc                                                                             >  1:2367645/1‑100 (MQ=255)
                       cGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGc                                                                             <  1:817203/100‑1 (MQ=255)
                         gCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCaa                                                                           >  1:3295589/1‑100 (MQ=255)
                         gCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCaa                                                                           >  2:353420/1‑100 (MQ=255)
                             gctgAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAg                                                                       >  1:2579469/1‑100 (MQ=255)
                                 aaGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTa                                                                   <  1:1631179/100‑1 (MQ=255)
                                  aGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAt                                                                  <  1:3331319/100‑1 (MQ=255)
                                     gTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAtttt                                                               >  1:11872/1‑100 (MQ=255)
                                     gTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAtttt                                                               >  1:214379/1‑100 (MQ=255)
                                     gTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAtttt                                                               >  2:3501927/1‑100 (MQ=255)
                                      tAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTc                                                              >  2:1400140/1‑100 (MQ=255)
                                       aaaaaCTTCGGCAAAGTGGTGGTCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCa                                                             <  2:152660/100‑1 (MQ=255)
                                       aaaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAttgt                                                               >  1:1443326/1‑96 (MQ=255)
                                       aaaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCa                                                             >  1:1784582/1‑100 (MQ=255)
                                         aaaTTTCGGCAAAGTTGTGATCCGCGTGGCGTGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAgn                                                           <  2:2425531/100‑2 (MQ=255)
                                         aaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGc                                                           <  2:1154985/100‑1 (MQ=255)
                                          aaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCt                                                          <  2:1292619/100‑1 (MQ=255)
                                          aaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCt                                                          >  1:3564179/1‑100 (MQ=255)
                                            cTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCttt                                                        >  2:1574201/1‑100 (MQ=255)
                                                ggCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaatt                                                    >  1:3418470/1‑100 (MQ=255)
                                                 gCAAAGTGGTGATCCGCGTGGCGGGTGATTATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaatt                                                    <  2:3553341/99‑1 (MQ=255)
                                                   aaaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaac                                                 >  1:2770316/1‑100 (MQ=255)
                                                    aaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaact                                                >  2:121803/1‑100 (MQ=255)
                                                       tggtgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaactaac                                             >  1:3224026/1‑100 (MQ=255)
                                                        ggtgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAAcc                                            >  2:144268/1‑100 (MQ=255)
                                                        ggtgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAAcc                                            >  1:1025822/1‑100 (MQ=255)
                                                        ggtgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAAcc                                            <  1:294726/100‑1 (MQ=255)
                                                        ggtgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAAcc                                            <  1:799100/100‑1 (MQ=255)
                                                          tgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcac                                          <  1:2430093/100‑1 (MQ=255)
                                                          tgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcac                                          >  2:1087488/1‑100 (MQ=255)
                                                           gATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcaca                                         >  1:541695/1‑100 (MQ=255)
                                                           gATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcaca                                         <  1:2708450/100‑1 (MQ=255)
                                                           gATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcaca                                         >  1:2559642/1‑100 (MQ=255)
                                                             tCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCAcatc                                       >  1:1686983/1‑100 (MQ=255)
                                                             tCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCAcatc                                       <  2:1674565/100‑1 (MQ=255)
                                                              ccGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCAcatca                                      >  2:2073085/1‑100 (MQ=255)
                                                               cgcgTGGCGGGTGATTATTAAAGTATCAACGGCGGCGTAACCCGCCATAAAAATCCAGTCAAACCGATTATTTTCAGCTTTAATTAACTAACCAcctca                                      >  2:2992493/1‑99 (MQ=255)
                                                                 cgTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTg                                   <  1:3416178/100‑1 (MQ=255)
                                                                  gTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGt                                  >  1:778481/1‑100 (MQ=255)
                                                                    ggCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTc                                <  1:430036/100‑1 (MQ=255)
                                                                      cGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctg                              >  1:368781/1‑100 (MQ=255)
                                                                        ggTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctgct                            <  1:510251/100‑1 (MQ=255)
                                                                        ggTGATGATAAAAGATTCAACGGCGGTGTAAGCCGCCATAAAAAAACAAGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctgct                            <  2:2502842/100‑1 (MQ=255)
                                                                           gatgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTgctgct                         <  1:179550/100‑1 (MQ=255)
                                                                            atgaaaaaaGTTTCAACGGCGGCGTAAGCCGCCATAAAAAAAAAAAAAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCta                        <  2:3047440/94‑1 (MQ=255)
                                                                             tgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaa                       <  2:3292493/100‑1 (MQ=255)
                                                                              tgccgatctttCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaat                      <  2:2209178/92‑1 (MQ=255)
                                                                              gatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaat                      >  1:1998830/1‑100 (MQ=255)
                                                                              gatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaat                      <  2:481703/100‑1 (MQ=255)
                                                                               atTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaata                     <  1:2073098/100‑1 (MQ=255)
                                                                               atTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaata                     <  2:3375130/100‑1 (MQ=255)
                                                                                   aaGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATaacaa                 <  1:106452/100‑1 (MQ=255)
                                                                                    aGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATaacaac                <  1:992503/100‑1 (MQ=255)
                                                                                    aGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATaacaac                <  1:3501952/100‑1 (MQ=255)
                                                                                     gattCAACGGCGGAGTAAGCCGCCATAAAAAGCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAacaaca               >  1:2453204/3‑100 (MQ=255)
                                                                                     gTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAacaaca               >  1:627615/1‑100 (MQ=255)
                                                                                     gTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAacaaca               >  2:2380711/1‑100 (MQ=255)
                                                                                      tttCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAAcaacaa              >  1:2873952/1‑100 (MQ=255)
                                                                                       ttCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaatagatcggaa             >  1:2209192/1‑92 (MQ=255)
                                                                                        tttatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAg            >  1:107726/6‑100 (MQ=255)
                                                                                         ttatggcggcTTACGCCTCCATAAAAATCCAGGCAACAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTTTTCTTCTGCCAATAACAACaaa             >  2:676741/5‑98 (MQ=255)
                                                                                         cAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAgg           <  1:679313/100‑1 (MQ=255)
                                                                                         cAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAgg           <  2:2891487/100‑1 (MQ=255)
                                                                                         cAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCGTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAg            <  2:2890238/99‑1 (MQ=255)
                                                                                          tatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGt          >  1:1513756/4‑100 (MQ=255)
                                                                                          tatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGt          >  1:1428247/4‑100 (MQ=255)
                                                                                           aCGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTg         <  2:2675813/100‑1 (MQ=255)
                                                                                            cggcggCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGa        >  2:3100954/1‑100 (MQ=255)
                                                                                            cggcggCGTAAGCAGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGa        >  2:2019053/1‑100 (MQ=255)
                                                                                             ggcggcGTAAGCCGCCATAAAAATCCAGGCAAAAAGCGTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg       >  1:1703519/1‑100 (MQ=255)
                                                                                              gcggcgTAAGCCGCCATAAAAATCCAAACAAAAAGCGTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGt      <  2:1703511/100‑1 (MQ=255)
                                                                                               cggcgTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTg     <  2:2940646/100‑1 (MQ=255)
                                                                                                ggcgTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCTCATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGt      <  2:1946155/98‑1 (MQ=255)
                                                                                                  cgTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGtt  >  1:1122127/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
GAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTT  >  NC_000913/1519981‑1520178

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