Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,678,937 T→A 56.3% Q84L (CAG→CTG)  yhjH ← cyclic‑di‑GMP phosphodiesterase, FlhDC‑regulated

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,678,9370TA56.3% 60.9 / 220.9 168Q84L (CAG→CTG) yhjHcyclic‑di‑GMP phosphodiesterase, FlhDC‑regulated
Reads supporting (aligned to +/- strand):  ref base T (33/40);  new base A (44/50);  total (77/91)
Fisher's exact test for biased strand distribution p-value = 8.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.46e-01

Only 100 of 169 total aligned reads displayed.

GCTGACGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACC  >  NC_000913/3678839‑3679035
                                                                                                  |                                                                                                  
gCTGACGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTg                                                                                                   >  2:1801886/1‑100 (MQ=255)
 cTGACGCAGGGCGATGAGCGTAGGGCCATCAATATTGCCCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGc                                                                                                  <  1:3172118/100‑1 (MQ=255)
 cTGACGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGc                                                                                                  <  1:1867015/100‑1 (MQ=255)
 cTGACGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGc                                                                                                  >  2:389894/1‑100 (MQ=255)
    aCGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGctct                                                                                               >  1:1712374/1‑100 (MQ=255)
     cGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCtt                                                                                              >  2:3556879/1‑100 (MQ=255)
     cGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCtt                                                                                              >  2:748499/1‑100 (MQ=255)
      gCAGGGCGATGGGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCt                                                                                                    <  1:3272139/93‑1 (MQ=255)
       cAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTc                                                                                            <  2:2467956/100‑1 (MQ=255)
        aGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTca                                                                                           <  1:422036/100‑1 (MQ=255)
        aGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTca                                                                                           <  1:1160262/100‑1 (MQ=255)
         gggCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTcac                                                                                          >  2:1385096/1‑100 (MQ=255)
         gggCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTcac                                                                                          >  2:1383934/1‑100 (MQ=255)
               tGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTc                                                                                    <  2:2791852/100‑1 (MQ=255)
                 aGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCa                                                                                  >  2:1601865/1‑100 (MQ=255)
                 aGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCa                                                                                  <  2:3254612/100‑1 (MQ=255)
                 aGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCa                                                                                  >  1:4157024/1‑100 (MQ=255)
                 aGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCa                                                                                  >  2:326158/1‑100 (MQ=255)
                  gCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCAt                                                                                 >  1:2893956/1‑100 (MQ=255)
                    ntAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATcc                                                                               >  2:3486000/2‑100 (MQ=255)
                     tAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCg                                                                              >  2:1128322/1‑100 (MQ=255)
                     tAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCg                                                                              <  2:1189516/100‑1 (MQ=255)
                      agcgcccttaaTATTGACCGATGCCAGCAGGCCGTGCTCTCTTAAGTAGTCGCCTTTTTGCGCCAGCCAATCAATCAGCTCTTTCACATCTTCCTTCCGa                                                                             <  1:1903178/91‑1 (MQ=255)
                           cATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTTTTTTTGTATTTCAATCTGCTCTTTCACAACCTCCATCCGATGGCt                                                                        >  1:1558358/1‑100 (MQ=255)
                            aTCAATATTGACTGATGTCAGCAGTCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGc                                                                         <  1:2581194/98‑1 (MQ=255)
                             tCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTACGCCAGCAAATCAATCAGCACTTTCACAACCTCCATCCGATGGCTg                                                                       >  1:2000406/1‑99 (MQ=255)
                               aaTATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACg                                                                    >  2:2964019/1‑100 (MQ=255)
                                atatTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACgg                                                                   <  1:4086665/100‑1 (MQ=255)
                                atatTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACgg                                                                   <  1:255384/100‑1 (MQ=255)
                                 tatTGACCGATGCCAGCAGGCCGTGCTATATAAAGAAGTCGGCTGTTTGAGCCAGCAAATCAATAAGCTCTTTCACAACCTCCCTCCGATGGCTGAggg                                                                   >  1:1570818/1‑99 (MQ=255)
                                     gACCGATGCCAGCAGGCCGTGCTCTATAAAGCAGTCGGCTTTGTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGAtt                                                              >  2:2881509/1‑100 (MQ=255)
                                     gACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGTCAGTAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGAtt                                                              >  1:709797/1‑100 (MQ=255)
                                     gACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGAtt                                                              >  1:3081081/1‑100 (MQ=255)
                                     gACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGAtt                                                              <  1:2876859/100‑1 (MQ=255)
                                      aCCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGAttt                                                             <  2:555516/100‑1 (MQ=255)
                                       ccGATGCCAGCAGGCCGTGCTCTATAAAGACGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTc                                                            <  1:3812409/100‑1 (MQ=255)
                                       ccGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTc                                                            >  1:3117807/1‑100 (MQ=255)
                                       ccGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCCCTTTCACAACCTCCAAACGATGGCTGACGGTGATTTc                                                            >  1:2077038/1‑100 (MQ=255)
                                       ccGACGCCAGCAGGCCGTGCCCTATAAAGAAGTCGGCTATTTGCGCCAGCAACTCAATCAGCTCATTCCCCCCCTCCATCCGATGGCTAACCGTGATTTc                                                            >  1:2698728/1‑100 (MQ=255)
                                         gATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAg                                                          >  2:4250450/1‑100 (MQ=255)
                                            gCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTaa                                                       <  2:409318/100‑1 (MQ=255)
                                             ccAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTaaa                                                      >  1:1037461/1‑100 (MQ=255)
                                             ccAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTaaa                                                      >  2:253965/1‑100 (MQ=255)
                                             ccAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTaaa                                                      >  1:1467228/1‑100 (MQ=255)
                                             ccAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTaaa                                                      >  1:2786438/1‑100 (MQ=255)
                                               agcagGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAAt                                                    <  2:751351/100‑1 (MQ=255)
                                               agcagGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAAt                                                    <  1:4202743/100‑1 (MQ=255)
                                                gcagGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATa                                                   <  2:2667881/100‑1 (MQ=255)
                                                gcagGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATa                                                   >  2:1748898/1‑100 (MQ=255)
                                                 cagGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAg                                                  >  2:525956/1‑100 (MQ=255)
                                                 cagGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAg                                                  <  2:2565057/100‑1 (MQ=255)
                                                  agGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGc                                                 <  1:1607778/100‑1 (MQ=255)
                                                  agGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCCTCCGATGGCTGACGGTGATTTCAGTAAAATAGc                                                 <  1:2799904/100‑1 (MQ=255)
                                                  agGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGc                                                 >  2:1179582/1‑100 (MQ=255)
                                                  agGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGc                                                 <  2:412961/100‑1 (MQ=255)
                                                  agGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTACATCCGATGGCTGACGGTGATTTCAGTAAAATAGc                                                 >  1:449871/1‑100 (MQ=255)
                                                   ggCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCg                                                >  1:637128/1‑100 (MQ=255)
                                                       gTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATcc                                            <  1:3453941/100‑1 (MQ=255)
                                                          ctctATAAAGAAGTGGGCTCTTTGCGCTAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATc                                             >  1:4084081/1‑96 (MQ=255)
                                                            ctATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCcggcgg                                       <  2:1010076/100‑1 (MQ=255)
                                                             tataAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCCATCGCTTACGGTCATTTCAGTAAAATAGCGATCCggcggc                                      >  1:1633210/1‑100 (MQ=255)
                                                             tataAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAGTAGCGATCCggcggc                                      >  2:4213949/1‑100 (MQ=255)
                                                              ataAAGAAGTCGGCTTTTTTTTTTAGTTAATCAATCAGCTCTTTCACAACCTCCATCCTATGGCTGACGGTGATTTCAGTAAAATAGCGATCTGGCGGCa                                     >  1:3868658/1‑100 (MQ=255)
                                                              ataAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCAGCCGATGGCTGAAGGTGATTTCACTAAAATAGCGATACGGCGGCa                                     >  1:529486/1‑100 (MQ=255)
                                                               taAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAg                                    <  2:189485/100‑1 (MQ=255)
                                                               taAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAg                                    >  1:547406/1‑100 (MQ=255)
                                                                 aagaagTCTGCTTTTTGCGCCAGTAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGc                                  <  1:943037/100‑1 (MQ=255)
                                                                   gaagTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGt                                <  1:1228022/100‑1 (MQ=255)
                                                                   gaagTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGt                                <  2:414435/100‑1 (MQ=255)
                                                                     agTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTg                              >  2:3299166/1‑100 (MQ=255)
                                                                      gTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGc                             >  1:828048/1‑100 (MQ=255)
                                                                      gTCGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGc                             <  1:1007492/100‑1 (MQ=255)
                                                                        cGGCTTTTTGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGa                           >  2:3734886/1‑100 (MQ=255)
                                                                          gCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGCTCCGGCGGCAGGCGTTGCGAAg                         <  1:3556866/100‑1 (MQ=255)
                                                                           cTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAgg                        >  2:4064243/1‑100 (MQ=255)
                                                                             ttttGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGt                      <  1:4207121/100‑1 (MQ=255)
                                                                             ttttGCGCCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGt                      <  2:1467232/100‑1 (MQ=255)
                                                                                 gcgcCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCaa                  >  2:3074819/1‑100 (MQ=255)
                                                                                  cgcCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAg                 >  1:2602015/1‑100 (MQ=255)
                                                                                  cgcCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAg                 >  2:760280/1‑100 (MQ=255)
                                                                                  cgcCAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAg                 <  2:3990370/100‑1 (MQ=255)
                                                                                  cgcCAGCAAATCAATCAGCTCTTTCACAACCGCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAg                 <  1:1044547/100‑1 (MQ=255)
                                                                                    ccAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAggg               <  1:4236686/100‑1 (MQ=255)
                                                                                    ccAGCAAATCAATCTGATCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAggg               >  1:3503253/1‑100 (MQ=255)
                                                                                    ccAGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGCTCCGGCGGCAGGCGTTGCGAAGGGTTCAAggg               >  1:39191/1‑100 (MQ=255)
                                                                                     cAGCAAATAAATCCGCTGTTTCACAACCTCCCTCCGATGGCTGAAGGTGATTTCAGTAAAATAGCGATCCGTCGGCAGGCGTTGCGAAGGGTTCAAGGGa              <  1:786886/100‑1 (MQ=255)
                                                                                      aGCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAGGGGTTCAAGTGAt             <  1:2472379/100‑1 (MQ=255)
                                                                                       gCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATg            >  2:3749590/1‑100 (MQ=255)
                                                                                       gCAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATg            <  2:1642650/100‑1 (MQ=255)
                                                                                        cAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGc           >  2:1314930/1‑100 (MQ=255)
                                                                                        cAAATCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGc           >  1:3161502/1‑100 (MQ=255)
                                                                                         aaaTCAATCAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGCTTCCAGTAAAATAGCGATCCGGCGGCAGGCGGTGCGAAGGGTTCACGGGATGCg          >  2:4133484/1‑100 (MQ=255)
                                                                                           atcaatcTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTg        <  2:3362191/100‑1 (MQ=255)
                                                                                           atcaatcTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTg        <  2:2097267/100‑1 (MQ=255)
                                                                                            tcaatcTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGa       <  2:839855/100‑1 (MQ=255)
                                                                                             caatcTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGac      >  2:2404092/1‑100 (MQ=255)
                                                                                             caatcAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGac      <  1:2405918/100‑1 (MQ=255)
                                                                                             caatcAGCTCTTTCACAACCCCCATCCGGCGTCTGTCGCTGATTTAAGTAAAATAGCGATCCGGCCGCAGGCGTTGCGTAGGGTTCAAGGGATGCGTGat      <  1:188824/100‑2 (MQ=255)
                                                                                              aatcAGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGacc     >  2:3630249/1‑100 (MQ=255)
                                                                                                 cTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGaccacc  <  1:2953560/100‑1 (MQ=255)
                                                                                                  |                                                                                                  
GCTGACGCAGGGCGATGAGCGTAGGGCCATCAATATTGACCGATGCCAGCAGGCCGTGCTCTATAAAGAAGTCGGCTTTTTGCGCCAGCAAATCAATCTGCTCTTTCACAACCTCCATCCGATGGCTGACGGTGATTTCAGTAAAATAGCGATCCGGCGGCAGGCGTTGCGAAGGGTTCAAGGGATGCGTGACCACC  >  NC_000913/3678839‑3679035

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