Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,520,080 G→T 15.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,0800GT15.1% 436.6 / 62.4 160intergenic (+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 (72/64);  new base T (6/18);  total (78/82)
Fisher's exact test for biased strand distribution p-value = 1.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.89e-01

Only 100 of 171 total aligned reads displayed.

GAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTT  >  NC_000913/1519981‑1520179
                                                                                                   |                                                                                                   
gagaATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGgcggcg                                                                                                     >  1:1565849/1‑100 (MQ=255)
  gaATGGGCCACATACGTTTATCGGCCTGCGGAAGGGTAAGAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGACTAAAGTTTCAACCTCGGCGTa                                                                                                   <  2:1394163/100‑1 (MQ=255)
    aTGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAg                                                                                                 >  2:171205/1‑100 (MQ=255)
    aTGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAg                                                                                                 >  1:409976/1‑100 (MQ=255)
      gcgcCACAGACGTTTATCGGCCTGCTGAAAGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATTATTAAAGTTTCAACGGCGGCGTAAg                                                                                                 <  2:1585115/98‑1 (MQ=255)
        gcCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgc                                                                                             <  2:3100343/100‑1 (MQ=255)
        acaacaGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCCAAGAGGTGCTCCGCGTGGCGGGTGATGGTTAAACTTTCAACGGCGGCGTACgccg                                                                                              >  2:3838850/4‑99 (MQ=255)
         ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgcc                                                                                            <  1:704268/100‑1 (MQ=255)
         ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgcc                                                                                            <  2:1766767/100‑1 (MQ=255)
         ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgcc                                                                                            >  2:2883954/1‑100 (MQ=255)
         ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgcc                                                                                            <  2:3044131/100‑1 (MQ=255)
           acaGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgt                                                                                          <  1:2428467/100‑3 (MQ=255)
           acaGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGCTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgt                                                                                          <  2:2103571/100‑3 (MQ=255)
            caGAGGTTTGTCGGCGTGCTGGGGGGTAGAAACTTCGGCAAAGGGGGGACGCGCGTGGCGGGTGATGATTGTAGTTGCCACGGCGGCTTACgccggcgtt                                                                                         <  2:104327/100‑6 (MQ=255)
            caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                         <  1:1078812/100‑4 (MQ=255)
            caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                         <  1:2081913/100‑4 (MQ=255)
            caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                         <  2:3358763/100‑4 (MQ=255)
            caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATa                                                                                         <  1:3423949/100‑1 (MQ=255)
            caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGAGCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATa                                                                                         <  2:2266787/100‑1 (MQ=255)
             aGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgttg                                                                                        <  1:269553/100‑5 (MQ=255)
              gACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATaaa                                                                                       >  2:2955961/1‑100 (MQ=255)
              gACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCCTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATaaa                                                                                       <  2:2240369/100‑1 (MQ=255)
                cGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATaaaaa                                                                                     >  1:86516/1‑100 (MQ=255)
                 gTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGTGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAAt                                                                                    <  2:2195469/100‑1 (MQ=255)
                 gTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAATGGCGGCTTACgccgccgttgaaac                                                                                    <  1:2958331/100‑9 (MQ=38)
                 gTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgttgaaac                                                                                    <  2:1458402/100‑9 (MQ=255)
                 gTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAAt                                                                                    >  1:995264/1‑100 (MQ=255)
                   ttATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATcc                                                                                  <  1:1228627/100‑1 (MQ=255)
                    tATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCa                                                                                 >  1:1604733/1‑100 (MQ=255)
                     aTCGGCCTGCTGAAGGGTAAAAACCTCCCCAAACTGGTCCTCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAGAAAt                                                                                    <  2:3580603/96‑1 (MQ=255)
                       cGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAgg                                                                               <  2:2695904/99‑1 (MQ=255)
                            tgctgAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCaaaaa                                                                         >  2:3759497/1‑100 (MQ=255)
                              ctgAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAgn                                                                       <  1:805303/100‑2 (MQ=255)
                              ctgAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGc                                                                       <  1:2208474/100‑1 (MQ=255)
                              ctgAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGc                                                                       >  1:666993/1‑100 (MQ=255)
                                 aaGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCtt                                                                     <  1:592695/99‑1 (MQ=255)
                                 aaGGGTAAAAACCTCGGCAATGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTGAACGGCGGCGTGAGCCGGCTTAAAAACCCAGGCAAAAAGCTTa                                                                    >  2:32810/1‑100 (MQ=255)
                                  aGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCACAAAGCTTa                                                                    >  2:3355053/1‑99 (MQ=255)
                                   gggTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAtt                                                                  >  1:2947763/1‑100 (MQ=255)
                                    ggTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAttt                                                                 >  1:3698055/1‑100 (MQ=255)
                                    ggTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAttt                                                                 >  2:431113/1‑100 (MQ=255)
                                    ggTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAttt                                                                 >  2:2174539/1‑100 (MQ=255)
                                     gTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAtttt                                                                <  1:3759525/100‑1 (MQ=255)
                                     gTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAtttt                                                                >  1:3215059/1‑100 (MQ=255)
                                     gTAAAAACTTCGGCAAAGTGGTAATCCGCGTGGCGAGTAAAGATTAAAATTTCAACAGCGGCGTAAGCCGCCATAAAAAAAAAAAAAAAAAGCTTAttt                                                                 <  2:837871/99‑1 (MQ=255)
                                      tAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTc                                                               >  1:1993380/1‑100 (MQ=255)
                                      tAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTc                                                               <  1:565148/100‑1 (MQ=255)
                                      tAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTc                                                               >  1:2404852/1‑100 (MQ=255)
                                       aaaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCa                                                              <  1:619161/100‑1 (MQ=255)
                                       aaaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCa                                                              <  2:1731441/100‑1 (MQ=255)
                                        aaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAg                                                             >  1:1885968/1‑100 (MQ=255)
                                        aaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAg                                                             >  1:1521497/1‑100 (MQ=255)
                                        aaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAg                                                             >  1:2471536/1‑100 (MQ=255)
                                         aaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGc                                                            >  2:2449030/1‑100 (MQ=255)
                                         aaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGc                                                            <  1:2754184/100‑1 (MQ=255)
                                            cTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCttt                                                         <  1:2633929/100‑1 (MQ=255)
                                             ttCGGCAAAGTGGTGATCCGCGTGGCGGGTAATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAAAAAAAGAAAAAAGCTTATTTTCAGCTtta                                                        <  2:86516/100‑1 (MQ=255)
                                                ggCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaatt                                                     <  2:2442507/100‑1 (MQ=255)
                                                ggCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaa                                                       >  2:312158/1‑98 (MQ=255)
                                                 gCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaatta                                                    >  2:3537201/1‑100 (MQ=255)
                                                 gCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaatta                                                    >  1:3127196/1‑100 (MQ=255)
                                                 gCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaatta                                                    >  1:2834748/1‑100 (MQ=255)
                                                  cAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaattaa                                                   >  1:1472921/1‑100 (MQ=255)
                                                   aaaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAACAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaac                                                  <  1:2804431/100‑1 (MQ=255)
                                                   aaaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaac                                                  <  2:1007673/100‑1 (MQ=255)
                                                   aaaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaac                                                  <  2:3588968/100‑1 (MQ=255)
                                                   aaaGTGGTGATCCGCGTGGCGGGTGATGACTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAAAAAAAAAAAAAAGCTTATTTTCAGCTTTAATtaac                                                  <  2:3539602/100‑1 (MQ=255)
                                                       tggtgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaactaac                                              <  1:823730/100‑1 (MQ=255)
                                                           gATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcaca                                          >  2:1224896/1‑100 (MQ=255)
                                                           gATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcaca                                          >  2:2365733/1‑100 (MQ=255)
                                                           gATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcaca                                          >  2:591077/1‑100 (MQ=255)
                                                           gATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcaca                                          >  1:763402/1‑100 (MQ=255)
                                                                  gTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGt                                   >  1:1002842/1‑100 (MQ=255)
                                                                  gTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGt                                   >  2:2824024/1‑100 (MQ=255)
                                                                     gCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTct                                >  1:2039928/1‑100 (MQ=255)
                                                                      cGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctg                               >  1:1647744/1‑100 (MQ=255)
                                                                        ggTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctgct                             <  2:1206426/100‑1 (MQ=255)
                                                                        ggTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCATTAATTAACTAACCACATCATTGTTctgc                              >  1:1169976/1‑99 (MQ=255)
                                                                         gTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctgctg                            <  2:1993370/100‑1 (MQ=255)
                                                                          tgatgaTTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCtgctgc                           <  2:3127176/100‑1 (MQ=255)
                                                                           gatgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTgctgct                          <  1:2706429/100‑1 (MQ=255)
                                                                           gatgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTgctgct                          <  1:1224905/100‑1 (MQ=255)
                                                                            atgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAACGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCta                         <  2:3433773/100‑1 (MQ=255)
                                                                            atgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCta                         <  1:1843992/100‑1 (MQ=255)
                                                                            aagatTAAATTTTAAACGGCGGCGTAATCCGCCATAAAAAAAAAAAAAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTgctgct                          <  2:2552026/97‑1 (MQ=255)
                                                                             tgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaa                        <  1:1843616/100‑1 (MQ=255)
                                                                                ttAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaataa                     <  1:3471152/100‑1 (MQ=255)
                                                                                  aaaGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATaac                    <  1:3051335/99‑1 (MQ=255)
                                                                                     gTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAacaaca                >  1:483451/1‑100 (MQ=255)
                                                                                       ttttatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACaaa              >  1:3267414/7‑100 (MQ=255)
                                                                                       ttttatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACaaa              >  1:1349758/7‑100 (MQ=255)
                                                                                       ttCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACaaa              >  2:2739666/1‑100 (MQ=255)
                                                                                        tCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAg             >  1:2645360/1‑100 (MQ=255)
                                                                                         cAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAgg            >  1:1831751/1‑100 (MQ=255)
                                                                                          aaCGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACGAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGt           >  2:2142262/1‑100 (MQ=255)
                                                                                           atggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAggag          >  1:2321656/3‑98 (MQ=255)
                                                                                               cggcgTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTg      >  2:739582/1‑100 (MQ=255)
                                                                                                ggcgTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTgg     <  2:1472912/100‑1 (MQ=255)
                                                                                                 gcgTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGt    >  1:1273324/1‑100 (MQ=255)
                                                                                                 gcgTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAAATAACAAAATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTgg     >  2:3147051/1‑99 (MQ=255)
                                                                                                   |                                                                                                   
GAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTT  >  NC_000913/1519981‑1520179

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: