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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,520,0830AC15.0% 414.2 / 66.5 165intergenic (+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 (73/67);  new base C (7/18);  total (80/85)
Fisher's exact test for biased strand distribution p-value = 3.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.17e-01

Only 100 of 169 total aligned reads displayed.

ATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATG  >  NC_000913/1519985‑1520182
                                                                                                  |                                                                                                   
aTGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAg                                                                                                    >  1:409976/1‑100 (MQ=255)
aTGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAg                                                                                                    >  2:171205/1‑100 (MQ=255)
 tGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgc                                                                                                   <  1:1581632/100‑1 (MQ=255)
   cgcCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccg                                                                                                 >  2:3583626/1‑100 (MQ=255)
    gcCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgc                                                                                                <  2:3100343/100‑1 (MQ=255)
     ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgcc                                                                                               <  1:3245912/100‑1 (MQ=255)
     ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgcc                                                                                               <  2:3044131/100‑1 (MQ=255)
     ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgcc                                                                                               <  1:704268/100‑1 (MQ=255)
     ccACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAgccgcc                                                                                               >  2:2883954/1‑100 (MQ=255)
       acaGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGCTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgt                                                                                             <  2:2103571/100‑3 (MQ=255)
        caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                            <  1:1078812/100‑4 (MQ=255)
        caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                            <  1:2832711/100‑4 (MQ=255)
        caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                            <  2:3358763/100‑4 (MQ=255)
        caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgtt                                                                                            <  2:197518/100‑4 (MQ=255)
        caGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGAGCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATa                                                                                            <  2:2266787/100‑1 (MQ=255)
         aGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgttg                                                                                           <  1:269553/100‑5 (MQ=255)
          gACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATaaa                                                                                          >  2:2955961/1‑100 (MQ=255)
          gACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATaaa                                                                                          <  1:2686425/100‑1 (MQ=255)
          gACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATCAAAGTTTCAACGGTGGCGTAAGCCGCCATaaa                                                                                          <  1:3740636/100‑1 (MQ=255)
           aCGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgttgaa                                                                                         <  1:196235/100‑7 (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)
               ttATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCTTACgccgccgttgaaactt                                                                                     <  2:2794099/100‑11 (MQ=38)
               ttATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATcc                                                                                     <  1:1228627/100‑1 (MQ=255)
                tATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCa                                                                                    >  1:1604733/1‑100 (MQ=255)
                   cGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAgg                                                                                  <  2:2695904/99‑1 (MQ=255)
                    ggCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCa                                                                                >  1:1863424/1‑100 (MQ=255)
                       ctgctgAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCaaaa                                                                             >  1:3421033/1‑100 (MQ=255)
                        tgctgAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCaaaaa                                                                            >  2:3759497/1‑100 (MQ=255)
                          ctgAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAgn                                                                          <  1:805303/100‑2 (MQ=255)
                          ctgAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGc                                                                          >  1:666993/1‑100 (MQ=255)
                             aaGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCtt                                                                        <  1:592695/99‑1 (MQ=255)
                             aaGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTa                                                                       <  2:3710079/100‑1 (MQ=255)
                              aGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCACAAAGCTTa                                                                       >  2:3355053/1‑99 (MQ=255)
                               gggTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAtt                                                                     >  1:2947763/1‑100 (MQ=255)
                                ggTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAttt                                                                    >  2:2174539/1‑100 (MQ=255)
                                 gTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAtttt                                                                   >  1:2971592/1‑100 (MQ=255)
                                 gTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAtttt                                                                   >  1:3215059/1‑100 (MQ=255)
                                 gTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAtttt                                                                   >  2:1622415/1‑100 (MQ=255)
                                 gTAAAAACTTCGGCAAAGTGGTGATCAGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTAtttt                                                                   >  1:2406375/1‑100 (MQ=255)
                                 gTAAAAACTTCGGCAAAGTGGTAATCCGCGTGGCGAGTAAAGATTAAAATTTCAACAGCGGCGTAAGCCGCCATAAAAAAAAAAAAAAAAAGCTTAttt                                                                    <  2:837871/99‑1 (MQ=255)
                                  tAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTc                                                                  >  2:1555693/1‑100 (MQ=255)
                                  tAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTc                                                                  <  1:565148/100‑1 (MQ=255)
                                   aaaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCa                                                                 >  1:2173384/1‑100 (MQ=255)
                                   aaaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCa                                                                 <  2:1731441/100‑1 (MQ=255)
                                   aaaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCa                                                                 <  1:619161/100‑1 (MQ=255)
                                   aaaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTATCa                                                                 >  2:1090359/1‑100 (MQ=255)
                                    aaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAg                                                                >  1:2471536/1‑100 (MQ=255)
                                    aaaaCTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAg                                                                >  1:1521497/1‑100 (MQ=255)
                                        cTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCttt                                                            >  1:1842257/1‑100 (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                                                       >  2:2485105/1‑100 (MQ=255)
                                             gCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaatta                                                       <  1:2028726/100‑1 (MQ=255)
                                             gCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaatta                                                       >  2:3908842/1‑100 (MQ=255)
                                             gCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaatta                                                       >  1:2834748/1‑100 (MQ=255)
                                              cAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTttaattaa                                                      >  1:1472921/1‑100 (MQ=255)
                                               aaaGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaac                                                     <  2:1007673/100‑1 (MQ=255)
                                               aaaGTGGTGATCCGCGTGGCGGGTGATGACTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAAAAAAAAAAAAAAGCTTATTTTCAGCTTTAATtaac                                                     <  2:3539602/100‑1 (MQ=255)
                                                   tggtgATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATtaactaac                                                 <  1:823730/100‑1 (MQ=255)
                                                       gATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcaca                                             >  2:591077/1‑100 (MQ=255)
                                                       gATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcaca                                             >  2:2365733/1‑100 (MQ=255)
                                                       gATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACcaca                                             >  2:1224896/1‑100 (MQ=255)
                                                            gcgTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCAtt                                        >  1:3061336/1‑100 (MQ=255)
                                                              gTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGt                                      >  2:2824024/1‑100 (MQ=255)
                                                              gTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGt                                      >  1:1002842/1‑100 (MQ=255)
                                                               tGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGtt                                     >  2:1432991/1‑100 (MQ=255)
                                                                 gCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTct                                   >  1:2039928/1‑100 (MQ=255)
                                                                    ggTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCATTAATTAACTAACCACATCATTGTTctgc                                 >  1:1169976/1‑99 (MQ=255)
                                                                     gTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTctgctg                               <  2:1993370/100‑1 (MQ=255)
                                                                      tgatgaTTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCGGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCtgctgc                              <  1:1690097/100‑1 (MQ=255)
                                                                      tgatgaTTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCtgctgc                              <  2:3127176/100‑1 (MQ=255)
                                                                       gatgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTgctgct                             <  1:1224905/100‑1 (MQ=255)
                                                                       gatgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTgctgct                             <  1:2706429/100‑1 (MQ=255)
                                                                        atgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAACGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCta                            <  2:3433773/100‑1 (MQ=255)
                                                                        aagatTAAATTTTAAACGGCGGCGTAATCCGCCATAAAAAAAAAAAAAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTgctgct                             <  2:2552026/97‑1 (MQ=255)
                                                                         tgatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaa                           <  1:1843616/100‑1 (MQ=255)
                                                                          gatTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaat                          >  1:2758803/1‑100 (MQ=255)
                                                                            ttAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaataa                        >  2:2259133/1‑100 (MQ=255)
                                                                            ttAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaataa                        <  1:3471152/100‑1 (MQ=255)
                                                                            ttAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCtaataa                        >  1:2268787/1‑100 (MQ=255)
                                                                             tAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATaac                       >  2:1029084/1‑100 (MQ=255)
                                                                              aaaGTTTCAGCGGGGTCGTAAGCTGCGATAAAGAGCCAGGCAAAAAGTTTATTTGCAGCTTTAATTGACTAACCACATCTTTGTTCTGCTGCtaa                           <  2:2155072/95‑1 (MQ=255)
                                                                                 gTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAacaaca                   >  1:483451/1‑100 (MQ=255)
                                                                                   ttttatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACaaa                 >  1:1349758/7‑100 (MQ=255)
                                                                                   ttCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACaaa                 >  1:825550/1‑100 (MQ=255)
                                                                                    tttatggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAg                >  1:2471148/6‑100 (MQ=255)
                                                                                       atggcggcTTACGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAggag             >  1:2321656/3‑98 (MQ=255)
                                                                                            ggcgTAATCCGCCATAAAAAAAAAAAAAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTg         <  2:468358/99‑1 (MQ=255)
                                                                                             gcgTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGt       >  1:1273324/1‑100 (MQ=255)
                                                                                             gcgTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAAATAACAAAATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTgg        >  2:3147051/1‑99 (MQ=255)
                                                                                              cgTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGtt      <  1:755306/100‑1 (MQ=255)
                                                                                               gTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGttt     >  1:3661069/1‑100 (MQ=255)
                                                                                                 aaGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTAt   <  2:3244013/100‑1 (MQ=255)
                                                                                                 aaGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTAt   >  1:3788771/1‑100 (MQ=255)
                                                                                                  |                                                                                                   
ATGCGCCACAGACGTTTATCGGCCTGCTGAAGGGTAAAAACTTCGGCAAAGTGGTGATCCGCGTGGCGGGTGATGATTAAAGTTTCAACGGCGGCGTAAGCCGCCATAAAAATCCAGGCAAAAAGCTTATTTTCAGCTTTAATTAACTAACCACATCATTGTTCTGCTGCTAATAACAACAAAGGTGAGTGGTTTATG  >  NC_000913/1519985‑1520182

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: