Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,520,080 G→T 19.3% 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,0800GT19.3% 525.5 / 107.1 217intergenic (+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 (76/99);  new base T (19/23);  total (95/122)
Fisher's exact test for biased strand distribution p-value = 8.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.63e-01

Only 100 of 225 total aligned reads displayed.

GAGAATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTT  >  NC_000913/1519981‑1520178
                                                                                                   |                                                                                                  
gagaATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGgcggcg                                                                                                    >  2:416253/1‑100 (MQ=255)
    aTGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAg                                                                                                <  2:2263734/100‑1 (MQ=255)
     tGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgc                                                                                               >  2:1588047/1‑100 (MQ=255)
     tGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgc                                                                                               <  2:3991851/100‑1 (MQ=255)
       cgcCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccg                                                                                             <  1:3244812/100‑1 (MQ=255)
        gcCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgc                                                                                            <  2:1974042/100‑1 (MQ=255)
         ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgcc                                                                                           >  1:584361/1‑100 (MQ=255)
         ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgcc                                                                                           <  2:3436862/100‑1 (MQ=255)
         ccACAGACGTTTATCGGCCTGCAGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGTGAGCTTAATCAAGTTTAAACGGCGGCGTAAgcag                                                                                             >  1:3039761/1‑96 (MQ=255)
          cacaGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTTAACGGCGGCGTAAGCCGCCa                                                                                          <  1:2688030/100‑1 (MQ=255)
          cacaGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCa                                                                                          <  2:1300981/100‑1 (MQ=255)
          cacaGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCa                                                                                          <  1:1375445/100‑1 (MQ=255)
          cacaGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCa                                                                                          <  1:126214/100‑1 (MQ=255)
           acaGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgt                                                                                         <  1:752409/100‑3 (MQ=255)
           acaGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgt                                                                                         <  1:4185423/100‑3 (MQ=255)
           acaGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgt                                                                                         <  2:2807709/100‑3 (MQ=255)
            caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                        <  2:1617628/100‑4 (MQ=255)
            caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                        <  2:1048235/100‑4 (MQ=255)
            caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                        <  2:1126965/100‑4 (MQ=255)
            caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                        <  2:4259170/100‑4 (MQ=255)
              gACGTTTATCGACCTGGTGACTGGTAAAAACTGCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                        <  1:112988/98‑4 (MQ=255)
               aCGGTTATCGCCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAATGTTTCAACGGCGGCTTACgccgccgttgaa                                                                                     <  1:2899707/100‑7 (MQ=255)
                cGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATaaaaa                                                                                    >  2:708257/1‑100 (MQ=255)
                 gTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAATGGCGGCTTACgccgccgttgaaac                                                                                   <  1:2071584/100‑9 (MQ=38)
                 gTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgttgaaac                                                                                   <  1:3447530/100‑9 (MQ=255)
                 gTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgttgaaac                                                                                   <  2:547371/100‑9 (MQ=255)
                   ttATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATcc                                                                                 <  1:1931606/100‑1 (MQ=255)
                    tATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCa                                                                                >  2:4122988/1‑100 (MQ=255)
                    tATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGGGGTGATCCGCGTGGCGGGTGATGATTAAAGCTTCAACGGCGGCGTAAGCCGCCATAAAAATCCa                                                                                <  1:2639705/100‑1 (MQ=255)
                     aTCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAg                                                                               <  2:4230464/100‑1 (MQ=255)
                     aTCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAg                                                                               >  2:1483189/1‑100 (MQ=255)
                       cGGCCTGCTGAAGGGTAAAAACTTCGGCCAAGTGGTGATCAGCGTGGCGGGTGATGATTAAAGTTTCACCGGCGGCGTAAGCCGCCATAAAAATCCAGGc                                                                             <  1:364434/100‑1 (MQ=255)
                       cGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGc                                                                             <  2:3054317/100‑1 (MQ=255)
                        ggCCTGCTGCAGGGTAAAAACTTCGGCAAAGTGGTGCTCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGTCGCCATAAAAATCCAGGCa                                                                            <  1:4336237/100‑1 (MQ=255)
                              ctgAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGc                                                                      <  1:3063866/100‑1 (MQ=255)
                                gAAGGGTAAAATCTTCGGCAAAGTGGTGATCCGCGTGGCGGGAGATGATTAGAGTTGCAACAGCGGCGTAGGCCGGCAGAGAAATTCATGCAAAAAGCtt                                                                    >  1:3529060/1‑100 (MQ=255)
                                gAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCtt                                                                    >  2:3743885/1‑100 (MQ=255)
                                    ggTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAttt                                                                >  2:1381926/1‑100 (MQ=255)
                                    ggTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGATTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAttt                                                                >  1:919062/1‑100 (MQ=255)
                                      tAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTc                                                              <  2:4301466/100‑1 (MQ=255)
                                          aaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCt                                                          >  1:2795523/1‑100 (MQ=255)
                                           actgccgCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCtt                                                         <  1:3070306/94‑1 (MQ=255)
                                           aCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCtt                                                         >  2:4101879/1‑100 (MQ=255)
                                            cTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCttt                                                        <  2:453381/100‑1 (MQ=255)
                                               cGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaat                                                     >  2:3609435/1‑100 (MQ=255)
                                                ggCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaatt                                                    <  2:3089566/100‑1 (MQ=255)
                                                ggCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaatt                                                    >  2:1910264/1‑100 (MQ=255)
                                                 gCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaatta                                                   >  2:4390587/1‑100 (MQ=255)
                                                 gCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAAACAAGACAAAAAGCTTATTTTCAGCTttaatta                                                   <  1:3319135/100‑1 (MQ=255)
                                                  cAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaattaa                                                  <  2:2773862/100‑1 (MQ=255)
                                                  cAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaatta                                                   <  1:3314588/99‑1 (MQ=255)
                                                   aaaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaac                                                 >  2:3156690/1‑100 (MQ=255)
                                                   aaaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAAAAAAAAAAAAAAGCTTATTTTCAGCTttaatta                                                   <  1:131893/98‑1 (MQ=255)
                                                      gtggtgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaactaa                                              >  2:907153/1‑100 (MQ=255)
                                                        ggtgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaactaac                                             <  1:4277862/99‑1 (MQ=255)
                                                         gtgATCCGCGTGGCGGGTGATGATTAAAGTTTCATCGGCGGCGTAAGCCGCCATAAAAATCCAAGAAAAAAGCTTATTTTCAGCTTTAATTAACTAAcc                                            <  1:3588230/99‑1 (MQ=255)
                                                          tgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcac                                          >  1:859180/1‑100 (MQ=255)
                                                          tgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcac                                          >  2:272858/1‑100 (MQ=255)
                                                          tgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcac                                          >  2:2706556/1‑100 (MQ=255)
                                                             tCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCAcatc                                       <  1:2620748/100‑1 (MQ=255)
                                                              ccacgTGGCGAGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGTTATAAAAAACCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcac                                          <  1:484200/93‑1 (MQ=255)
                                                              ccGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCAcatca                                      <  1:4066828/100‑1 (MQ=255)
                                                              ccGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCAcatca                                      <  2:1845252/100‑1 (MQ=255)
                                                                   tGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGtt                                 >  2:3853147/1‑100 (MQ=255)
                                                                    ggCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTc                                <  1:4069233/100‑1 (MQ=255)
                                                                       gggTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctgc                             <  2:1200202/100‑1 (MQ=255)
                                                                          tgatgaTTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAAACAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctgct                            <  1:2528554/98‑1 (MQ=255)
                                                                           gatgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTgctgct                         <  1:2280125/100‑1 (MQ=255)
                                                                            atgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCta                        <  2:3323449/100‑1 (MQ=255)
                                                                             tgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaa                       >  1:3575593/1‑100 (MQ=255)
                                                                             tgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTAAGCTTTAATTAACTAACCACATCATTGTTctgctg                           >  2:895539/1‑96 (MQ=255)
                                                                              gatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaat                      <  2:3923578/100‑1 (MQ=255)
                                                                              gatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaat                      <  2:588590/100‑1 (MQ=255)
                                                                              gatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaat                      <  2:1727041/100‑1 (MQ=255)
                                                                               gatctaGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaata                     <  2:1339472/95‑1 (MQ=255)
                                                                               atTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaata                     <  2:3466661/100‑1 (MQ=255)
                                                                                 tAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATaac                   <  2:1145763/100‑1 (MQ=255)
                                                                                 tAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCGTTAATTAACTAACCACATCATTGTTCTGCTGCtaataa                    >  1:628033/1‑99 (MQ=255)
                                                                                   aaGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATaacaa                 <  2:3347068/100‑1 (MQ=255)
                                                                                    ggtttttatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATaacaac                >  2:2987365/10‑100 (MQ=255)
                                                                                     gTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAacaaca               >  2:827370/1‑100 (MQ=255)
                                                                                     gTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAacaaca               >  2:2028371/1‑100 (MQ=255)
                                                                                     gTTTCAACGGCGACGTAAGCCGCCATAAAAATCCAAACAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTGCTGCTGCTAATaacaac                <  1:2117891/99‑1 (MQ=255)
                                                                                     atttttatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAacaaca               >  1:779964/9‑100 (MQ=255)
                                                                                       ttttatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACaaa             >  2:1562635/7‑100 (MQ=255)
                                                                                       ttttatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACaaa             >  1:1596567/7‑100 (MQ=255)
                                                                                        tCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAg            <  2:2420737/100‑1 (MQ=255)
                                                                                         ttatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAgg           >  2:3263270/5‑100 (MQ=255)
                                                                                         cAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAgg           <  1:1456359/100‑1 (MQ=255)
                                                                                          tatggcggcTTACGCCGCCATAAAAATCCTGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTTCTGCTAATAACAACAAAgg           >  1:1021494/4‑99 (MQ=255)
                                                                                          tatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTCTTTTCCGCCTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAgg           >  1:244218/4‑99 (MQ=255)
                                                                                           atggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTg         >  1:1786868/3‑100 (MQ=255)
                                                                                           atggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTg         >  1:2174135/3‑100 (MQ=255)
                                                                                           atggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTg         >  1:1390655/3‑100 (MQ=255)
                                                                                           atggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTg         >  1:2984417/3‑100 (MQ=255)
                                                                                            cggcggCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGa        >  1:1242092/1‑100 (MQ=255)
                                                                                             ggcggcGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg       <  2:649275/100‑1 (MQ=255)
                                                                                             ggcggcGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAg       <  2:2845105/100‑1 (MQ=255)
                                                                                              gcggcgTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTg         >  1:1891031/1‑97 (MQ=255)
                                                                                                 gcgTAAGCCGCCATAAACATCCAGGAAAAAAGCTTATATTCAGCTATAATTAACTAACCACAACATTATTCTGCTGCTAATAACAAAAAAGGTGAGTgg    >  1:2357421/1‑99 (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 < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: