Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,778,372 C→G 15.4% intergenic (+25/+18) ydiO → / ← ydiP putative acyl‑CoA dehydrogenase/putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,778,3720CG15.4% 400.2 / 58.9 143intergenic (+25/+18)ydiO/ydiPputative acyl‑CoA dehydrogenase/putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (65/55);  new base G (11/11);  total (76/67)
Fisher's exact test for biased strand distribution p-value = 8.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

Only 100 of 145 total aligned reads displayed.

ATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTA  >  NC_000913/1778273‑1778470
                                                                                                   |                                                                                                  
aTCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAAc                                                                                                    >  1:2505383/1‑100 (MQ=255)
 tCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACa                                                                                                   <  1:2819496/100‑1 (MQ=255)
  cggcggcggTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACaa                                                                                                  >  2:2537402/1‑100 (MQ=255)
  cggcggcggTACAGACGAAATTATGATTTACGGAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCGTCTTAACa                                                                                                   >  2:2473139/1‑99 (MQ=255)
   ggcggcggTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACaaa                                                                                                 >  1:643178/1‑100 (MQ=255)
    gcggcggTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAAc                                                                                                <  1:3283905/100‑1 (MQ=255)
       gcggTACAGACGAAATTATGATTTACGTAGCAGGTCGGAACATCCTGAAAGGCTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTgc                                                                                             <  1:47187/100‑1 (MQ=255)
         ggTACAGACGAAATTATGATTTACGTAGCTGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTgcgc                                                                                           >  2:1301849/1‑100 (MQ=255)
         ggTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTgcgc                                                                                           <  1:686868/100‑1 (MQ=255)
          gTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGcgcc                                                                                          <  1:782958/100‑1 (MQ=38)
          gTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGcgcc                                                                                          <  1:1011510/100‑1 (MQ=255)
           tACAGACGAAATTATGATTTACGTAGCAGGTCGGCATATCCTGAAAGACTATCAGAACAAATAGTCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGcgccg                                                                                         <  2:1709806/100‑1 (MQ=37)
           tACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGcgccg                                                                                         <  1:616617/100‑1 (MQ=38)
           tACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGcgccg                                                                                         <  1:3698408/100‑1 (MQ=38)
           tACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGcgccg                                                                                         <  2:1989573/100‑1 (MQ=38)
           tACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGcgccg                                                                                         <  1:287677/100‑1 (MQ=255)
           tACAGACGAAATTATGATTTACGTAGCAGATCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGcgccg                                                                                         <  1:1654663/100‑1 (MQ=38)
             cAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGATGCTGCgccgcc                                                                                       <  2:2489978/100‑1 (MQ=38)
               gACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAg                                                                                     >  1:2480788/1‑100 (MQ=255)
               gACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAACCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGCgccgcctg                                                                                     <  1:1296783/100‑3 (MQ=37)
                 cGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGTTTGTTAAGAAGCTGCgccgcctgca                                                                                   <  1:2457699/100‑5 (MQ=25)
                  gAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGAtt                                                                                  <  2:3212488/100‑1 (MQ=255)
                  gAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGAtt                                                                                  <  2:3211714/100‑1 (MQ=255)
                    aaTTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTa                                                                                >  1:1376191/1‑100 (MQ=255)
                     aTTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTAt                                                                               >  1:2155820/1‑100 (MQ=255)
                      ttATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATc                                                                              <  2:1286174/100‑1 (MQ=255)
                      ttATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATc                                                                              <  2:1382677/100‑1 (MQ=255)
                      ttATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATc                                                                              >  1:974258/1‑100 (MQ=255)
                      ttATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATc                                                                              <  2:3727820/100‑1 (MQ=255)
                       tATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATcc                                                                             >  1:670197/1‑100 (MQ=255)
                       tATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATcc                                                                             <  2:3774121/100‑1 (MQ=255)
                        aTGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCca                                                                            <  1:1004997/100‑1 (MQ=255)
                        aTGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCca                                                                            >  1:1230972/1‑100 (MQ=255)
                         tGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCcaa                                                                           <  1:2599008/100‑1 (MQ=255)
                          gATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCcaac                                                                          <  2:1230963/100‑1 (MQ=255)
                           aTTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCcaaca                                                                         <  2:2145884/100‑1 (MQ=255)
                           aTTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCcaaca                                                                         <  1:1585816/100‑1 (MQ=255)
                              tACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGa                                                                      >  1:1359162/1‑100 (MQ=255)
                              tACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGAGTTATCCAACAAg                                                                       >  1:2215907/1‑99 (MQ=255)
                                   aGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGGAGCCGGCGCAGCTTCGTAACAAACTGAGCCGCCAGAGTTATCCAACAAGAATTAc                                                                 >  2:2257661/1‑100 (MQ=255)
                                   aGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTAc                                                                 >  1:1498173/1‑100 (MQ=255)
                                     cAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCg                                                               >  2:439139/1‑100 (MQ=255)
                                     cAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCg                                                               >  1:3178479/1‑100 (MQ=255)
                                      aGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCgg                                                              >  1:1178685/1‑100 (MQ=255)
                                       ggTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGt                                                             >  1:765749/1‑100 (MQ=255)
                                            gctgATCCTGAAAGACTAGCAGAACAGATAGTCTGCCGGCGTCGCAGCTTCGTAACAAACGGCGCCGCCACATTTATCCAACAAGACTTACCTGTTg                                                           <  2:1693187/94‑1 (MQ=255)
                                             cAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGa                                                       <  1:1699877/100‑1 (MQ=255)
                                               gATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGaaa                                                     <  1:3530121/100‑1 (MQ=255)
                                                 tCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAAtt                                                   >  1:3292619/1‑100 (MQ=255)
                                                    tGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAg                                                <  2:2767899/100‑1 (MQ=255)
                                                      aaaGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGa                                               >  2:1178311/1‑99 (MQ=255)
                                                      aaaGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGAt                                              >  1:3775140/1‑100 (MQ=255)
                                                       aaGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATa                                             <  2:741399/100‑1 (MQ=255)
                                                       aaGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATa                                             <  1:2178321/100‑1 (MQ=255)
                                                        aGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAg                                            >  2:2288883/1‑100 (MQ=255)
                                                        aGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAg                                            <  1:3780600/100‑1 (MQ=255)
                                                           cTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGAt                                         <  1:1835898/100‑1 (MQ=255)
                                                             aTCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCCGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGAttt                                       >  2:2655561/1‑100 (MQ=255)
                                                             aTCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGAttt                                       >  2:1614225/1‑100 (MQ=255)
                                                                aGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCtt                                    >  1:3202615/1‑100 (MQ=255)
                                                                  aaCAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGt                                  >  2:793360/1‑100 (MQ=255)
                                                                  aaCAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGt                                  >  2:793506/1‑100 (MQ=255)
                                                                    cAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTc                                >  1:1496582/1‑100 (MQ=255)
                                                                    cAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTc                                >  2:3070315/1‑100 (MQ=255)
                                                                      aataatCTGCAGGCAGCGCAGCTTCTTAACAAACTGCGCCGCGAGGTTTATCCAACAAGACTTACCGGTTGATGAAATTCAGATAGGATTTCTTGTTCCg                              <  2:1201669/100‑1 (MQ=255)
                                                                       ggataaatctggcggcGCAGTTTGTTAAGAAGCTGCGCCGCCTGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGc                             >  2:3127964/11‑100 (MQ=18)
                                                                          atCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCaaa                          >  1:39996/1‑100 (MQ=255)
                                                                            cTGCAGGCGGCGCAGCTTCTTAACTATCTGCGCCGCCATATTTATCCAACAAGACTTACCGGTTGAGGAATGTCAGATAGGATTTCTTGTTCCGCAAAAc                        <  1:389347/100‑1 (MQ=255)
                                                                            cTGCAGGCGGCGCAGCTTCTTAACAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAAc                        >  1:1925311/1‑100 (MQ=255)
                                                                            cTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAAc                        <  2:643171/100‑1 (MQ=255)
                                                                            cTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAAc                        >  1:2042152/1‑100 (MQ=255)
                                                                            cTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAAc                        <  2:2389471/100‑1 (MQ=255)
                                                                             tGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACa                       >  2:542566/1‑100 (MQ=255)
                                                                               ctggcggcGCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGt                     >  1:818552/3‑100 (MQ=37)
                                                                               cAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGt                     >  2:3662062/1‑100 (MQ=255)
                                                                               cAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGt                     >  1:1047167/1‑100 (MQ=255)
                                                                                tggcggcGCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtt                    >  2:28157/2‑100 (MQ=37)
                                                                                 ggcggcGCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGttt                   >  2:2518916/1‑100 (MQ=38)
                                                                                 ggcggcGCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGtt                    >  2:697985/1‑99 (MQ=38)
                                                                                 ggcggcGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGttt                   <  2:746353/100‑1 (MQ=255)
                                                                                 ggcggcGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGttt                   >  2:1521822/1‑100 (MQ=255)
                                                                                   cggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGttttt                 >  1:3031326/1‑100 (MQ=38)
                                                                                    ggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTa                >  1:2500390/1‑100 (MQ=38)
                                                                                    ggcgCAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTa                >  1:1058297/1‑100 (MQ=38)
                                                                                    ggcgCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTa                >  1:785519/1‑100 (MQ=255)
                                                                                    ggagcAGCTTCTTAAAAACCTGCGCCGCCAGAATTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTa                <  1:2473155/97‑1 (MQ=255)
                                                                                     gcgcAGTTTGTTAAGAAGCTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTaa               >  2:3450635/1‑100 (MQ=38)
                                                                                       gcAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAAt             <  1:3568633/100‑1 (MQ=255)
                                                                                         aGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTg           <  2:2128546/100‑1 (MQ=255)
                                                                                           cttcttAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTc         >  1:544779/1‑100 (MQ=255)
                                                                                           catcttAACAAACAGCGCCGCCAGATTTATCAAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTc         >  1:2576538/3‑100 (MQ=255)
                                                                                            ttcttAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCt        <  2:601597/100‑1 (MQ=255)
                                                                                            ttcttAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCt        >  2:618413/1‑100 (MQ=255)
                                                                                               ttAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACg     >  1:885376/1‑100 (MQ=255)
                                                                                                tAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTAccg    <  1:2579384/100‑3 (MQ=255)
                                                                                                tAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACgg    <  2:3075425/100‑1 (MQ=255)
                                                                                                  aCAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTa  >  2:470618/1‑100 (MQ=255)
                                                                                                  aCAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTa  >  1:2639816/1‑100 (MQ=255)
                                                                                                  aCAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTa  >  2:2910115/1‑100 (MQ=255)
                                                                                                  aCAAACTACGCAGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTCGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTa  >  2:2000134/1‑100 (MQ=255)
                                                                                                   |                                                                                                  
ATCGGCGGCGGTACAGACGAAATTATGATTTACGTAGCAGGTCGGCAGATCCTGAAAGACTATCAGAACAAATAATCTGCAGGCGGCGCAGCTTCTTAACAAACTGCGCCGCCAGATTTATCCAACAAGACTTACCGGTTGAGGAAATTCAGATAGGATTTCTTGTTCCGCAAAACAGTTTTTAAATTGTCTACGGTA  >  NC_000913/1778273‑1778470

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