Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,520,083 A→C 16.1% 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,0830AC16.1% 461.6 / 78.1 188intergenic (+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 (79/78);  new base C (16/15);  total (95/93)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.88e-01

Only 100 of 191 total aligned reads displayed.

CGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATG  >  NC_000913/1519988‑1520182
                                                                                               |                                                                                                   
cgcCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccg                                                                                                 >  1:1328894/1‑100 (MQ=255)
  ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgcc                                                                                               <  1:2983809/100‑1 (MQ=255)
  ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgcc                                                                                               >  1:578421/1‑100 (MQ=255)
   cacaGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCa                                                                                              >  1:2754865/1‑100 (MQ=255)
    acaGACGTTTATCGGCCTGCTGAATGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgt                                                                                             <  2:2468106/100‑3 (MQ=255)
    acaGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgt                                                                                             <  1:3450221/100‑3 (MQ=255)
    acaGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgt                                                                                             <  1:1844550/100‑3 (MQ=255)
     cgggcGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCTAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                            <  2:2561429/96‑4 (MQ=255)
     caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                            <  2:586764/100‑4 (MQ=255)
     caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                            <  1:2084707/100‑4 (MQ=255)
     caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                            <  2:1359054/100‑4 (MQ=255)
     caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATa                                                                                            <  1:3344663/100‑1 (MQ=255)
     caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATa                                                                                            >  1:3070991/1‑100 (MQ=255)
     caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATa                                                                                            >  1:1676044/1‑100 (MQ=255)
      aGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgttg                                                                                           <  1:2283888/100‑5 (MQ=255)
       gACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATaaa                                                                                          <  1:3413864/100‑1 (MQ=255)
        aCGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATaaaa                                                                                         <  1:500815/100‑1 (MQ=255)
         cGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATaaaaa                                                                                        <  1:2855259/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)
             tATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCa                                                                                    <  2:2380586/100‑1 (MQ=255)
             tATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCa                                                                                    <  1:2553756/100‑1 (MQ=255)
              aTCGGCCTGCTGAAGGGTATAAACTTCGGCAAAGTGGTGAGCCGCGTGGCGTGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAg                                                                                   <  2:455305/100‑1 (MQ=255)
              aTCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAg                                                                                   >  1:3510421/1‑100 (MQ=255)
              aTCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAg                                                                                   <  2:1799499/100‑1 (MQ=255)
                cGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGc                                                                                 <  1:817203/100‑1 (MQ=255)
                cGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGc                                                                                 >  1:2367645/1‑100 (MQ=255)
                  gCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCaa                                                                               >  2:353420/1‑100 (MQ=255)
                              gTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAtttt                                                                   >  2:1006648/1‑100 (MQ=255)
                              gTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAtttt                                                                   >  1:11872/1‑100 (MQ=255)
                              gTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAtttt                                                                   >  2:3501927/1‑100 (MQ=255)
                              gTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAtttt                                                                   >  1:214379/1‑100 (MQ=255)
                               tAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTc                                                                  >  2:1400140/1‑100 (MQ=255)
                               tAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTc                                                                  >  2:294724/1‑100 (MQ=255)
                                aaaaaCTTCGGCAAAGTGGTGGTCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCa                                                                 <  2:152660/100‑1 (MQ=255)
                                aaaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAttgt                                                                   >  1:1443326/1‑96 (MQ=255)
                                aaaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCa                                                                 >  1:678368/1‑100 (MQ=255)
                                aaaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCa                                                                 >  1:1784582/1‑100 (MQ=255)
                                  aaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGc                                                               <  1:3416945/100‑1 (MQ=255)
                                  aaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGc                                                               <  2:1154985/100‑1 (MQ=255)
                                   aaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCt                                                              <  2:1292619/100‑1 (MQ=255)
                                    aCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCtt                                                             >  1:287405/1‑100 (MQ=255)
                                         ggCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTtta                                                           >  2:2183704/1‑97 (MQ=255)
                                          gCAAAGTGGTGATCCGCGTGGCGGGTGATTATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaatt                                                        <  2:3553341/99‑1 (MQ=255)
                                           cAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaattaa                                                      >  1:86365/1‑100 (MQ=255)
                                            aaaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaac                                                     <  1:3317461/100‑1 (MQ=255)
                                            aaaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaac                                                     <  1:2560440/100‑1 (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:799100/100‑1 (MQ=255)
                                                 ggtgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAAcc                                                >  1:1025822/1‑100 (MQ=255)
                                                   tgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcac                                              <  1:2430093/100‑1 (MQ=255)
                                                   tgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcac                                              >  2:1087488/1‑100 (MQ=255)
                                                      tCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCAcatc                                           <  2:1674565/100‑1 (MQ=255)
                                                       ccGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCAcatca                                          >  2:2073085/1‑100 (MQ=255)
                                                           gTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGt                                      >  1:778481/1‑100 (MQ=255)
                                                             ggCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTc                                    <  1:430036/100‑1 (MQ=255)
                                                               cGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctg                                  >  1:368781/1‑100 (MQ=255)
                                                                 ggTGATGATAAAAGATTCAACGGCGGTGTAAGCCGCCATAAAAAAACAAGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctgct                                <  2:2502842/100‑1 (MQ=255)
                                                                   tgatgaTTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCtgctgc                              <  1:3230887/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)
                                                                       gtctAAAGTTTCAACGGCGGCGTAAGCAGCCGTAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTTCTGCtaat                          >  1:2450263/4‑100 (MQ=255)
                                                                       gatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaat                          >  1:1998830/1‑100 (MQ=255)
                                                                       gatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaat                          <  2:481703/100‑1 (MQ=255)
                                                                        atTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaata                         <  1:2073098/100‑1 (MQ=255)
                                                                        atTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaata                         <  1:2613682/100‑1 (MQ=255)
                                                                        atTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaata                         <  2:3375130/100‑1 (MQ=255)
                                                                         ttAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaataa                        <  1:1695286/100‑1 (MQ=255)
                                                                            aaGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATaacaa                     <  1:106452/100‑1 (MQ=255)
                                                                              gTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAacaaca                   >  2:2969508/1‑100 (MQ=255)
                                                                              gTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAacaaca                   >  1:627615/1‑100 (MQ=255)
                                                                               tttttatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAAcaacaa                  >  1:6871/8‑100 (MQ=255)
                                                                               tttCAACGGCGGTGTAAGCCGCCATAAAAATCCAGGCTAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAAcaacaa                  <  2:2450246/100‑1 (MQ=255)
                                                                                ttCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaatagatcggaa                 >  1:2209192/1‑92 (MQ=255)
                                                                                 tttatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAg                >  1:107726/6‑100 (MQ=255)
                                                                                 tCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAg                >  1:3283357/1‑100 (MQ=255)
                                                                                  ttatggcggcTTACGCCTCCATAAAAATCCAGGCAACAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTTTTCTTCTGCCAATAACAACaaa                 >  2:676741/5‑98 (MQ=255)
                                                                                  cAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAgg               <  2:2891487/100‑1 (MQ=255)
                                                                                  cAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAgg               >  2:3306757/1‑100 (MQ=255)
                                                                                  cAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATAATTGTTCTGCTGCTAATAACAACAAAgg               >  2:3354187/1‑100 (MQ=255)
                                                                                  cAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCGTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAg                <  2:2890238/99‑1 (MQ=255)
                                                                                   tatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGt              >  1:1513756/4‑100 (MQ=255)
                                                                                   tatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGt              >  1:1428247/4‑100 (MQ=255)
                                                                                   tatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGt              >  1:2701851/4‑100 (MQ=255)
                                                                                   tatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGt              >  1:511889/4‑100 (MQ=255)
                                                                                   aaCGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGt              >  1:1756353/1‑100 (MQ=255)
                                                                                   aaCGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGt              >  2:2405939/1‑100 (MQ=255)
                                                                                    atggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTg             >  2:1625807/3‑100 (MQ=255)
                                                                                    aCGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTg             <  2:2675813/100‑1 (MQ=255)
                                                                                     cggcggCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGa            >  2:3100954/1‑100 (MQ=255)
                                                                                     cggcggCGTAAGCAGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGa            >  2:2019053/1‑100 (MQ=255)
                                                                                       gcggcgTAAGCCGCCATAAAAATCCAAACAAAAAGCGTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGt          <  2:1703511/100‑1 (MQ=255)
                                                                                         ggcgTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCTCATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGt          <  2:1946155/98‑1 (MQ=255)
                                                                                           cgTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGtt      >  1:1122127/1‑100 (MQ=255)
                                                                                               aGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATg  <  1:544387/100‑1 (MQ=255)
                                                                                               |                                                                                                   
CGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATG  >  NC_000913/1519988‑1520182

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: