Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,098,518 C→T 7.8% intergenic (+115/+40) yggW → / ← yggM HemN family putative oxidoreductase/DUF1202 family putative secreted protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,098,5180CT 7.8% 200.8 / 5.9 64intergenic (+115/+40)yggW/yggMHemN family putative oxidoreductase/DUF1202 family putative secreted protein
Reads supporting (aligned to +/- strand):  ref base C (29/30);  new base T (2/3);  total (31/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01

AAATTCGCTGCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTT  >  NC_000913/3098370‑3098665
                                                                                                                                                    |                                                                                                                                                   
aaaTTCGCTGCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTc                                                                                                                                                     >  1:417133/1‑149 (MQ=255)
      gctgctGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAAcgcc                                                                                                                                               >  4:99668/1‑149 (MQ=255)
            ggAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATcc                                                                                                                                         <  1:542889/149‑1 (MQ=255)
               gCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCAGGATGCGGCGTTGACCCCGCATCCGGc                                                                                                                                      <  2:401137/149‑1 (MQ=25)
                      cTGGCTGAGTAAACTTGTATTGCCGGATGCGGCATGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAAt                                                                                                                               <  4:92191/149‑1 (MQ=255)
                         gCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGt                                                                                                                            <  4:584729/149‑1 (MQ=255)
                         gCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGt                                                                                                                            <  3:816099/149‑1 (MQ=255)
                         gCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGt                                                                                                                            <  2:417133/149‑1 (MQ=255)
                          cTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTg                                                                                                                           <  3:596184/149‑1 (MQ=255)
                           tGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGc                                                                                                                          >  4:584509/1‑149 (MQ=255)
                              gTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACt                                                                                                                       >  4:555160/1‑149 (MQ=255)
                               tAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACta                                                                                                                      <  1:309535/149‑1 (MQ=255)
                                aaaCTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACtat                                                                                                                     <  1:856494/149‑1 (MQ=255)
                                  aCTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTatat                                                                                                                   >  3:87070/1‑149 (MQ=255)
                                    ttGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCa                                                                                                                 >  3:545214/1‑149 (MQ=255)
                                    ttGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCa                                                                                                                 >  4:860147/1‑149 (MQ=255)
                                    ttGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCa                                                                                                                 >  4:620443/1‑149 (MQ=255)
                                    ttGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCa                                                                                                                 >  3:120098/1‑149 (MQ=255)
                                    ttGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCa                                                                                                                 >  2:210708/1‑149 (MQ=255)
                                      gTATTGCCGGATGCGGCGTAAACGCCTTATCCATCCGACCTGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAACCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCaga                                                                                                               <  2:179197/149‑3 (MQ=255)
                                       tATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAac                                                                                                              <  2:823312/149‑1 (MQ=255)
                                           gCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAacttac                                                                                                          >  2:780458/1‑149 (MQ=255)
                                            ccGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAacttact                                                                                                         >  1:441165/1‑149 (MQ=255)
                                                aTGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAAt                                                                                                     <  3:620443/149‑1 (MQ=255)
                                                  gcggcgTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTACACTATATCAAACTTACTTAATCa                                                                                                   <  1:745345/149‑1 (MQ=255)
                                                    ggcgTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCAGGATGCGGCGTTGACCCCGCATCCGGca                                                                                                                                     >  4:104726/1‑112 (MQ=18)
                                                    ggcgTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCAGGATGCGGCGTTGACCCCGCATCCGGca                                                                                                                                     <  3:104726/113‑2 (MQ=18)
                                                     gcgTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCAccc                                                                                                <  1:210708/149‑1 (MQ=255)
                                                     gcgTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCAccc                                                                                                <  4:574887/149‑1 (MQ=255)
                                                          aaCGCCTTATCCAGCCGACACGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTg                                                                                           >  2:100401/1‑149 (MQ=255)
                                                             gCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTgcgc                                                                                        >  4:311505/1‑149 (MQ=255)
                                                             gCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTgcgc                                                                                        >  2:15230/1‑149 (MQ=255)
                                                                 tATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGATGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCa                                                                                    <  4:325630/149‑1 (MQ=255)
                                                                   tCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAgt                                                                                  >  2:838293/1‑149 (MQ=255)
                                                                      aGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAgttgt                                                                               >  3:114698/1‑149 (MQ=255)
                                                                        ccGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTc                                                                             <  2:789107/149‑1 (MQ=255)
                                                                         cGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCg                                                                            <  3:95126/149‑1 (MQ=255)
                                                                          gACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCgg                                                                           >  2:123977/1‑149 (MQ=255)
                                                                            cATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTa                                                                         <  4:120098/149‑1 (MQ=255)
                                                                               gtgGCAGCGGTAGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCAGCCTGCTACAAATCGTGCACGATATCAATCTTACTTAATCACCCCATTGCGCTGCAGTTGTTCGGTATCt                                                                      <  4:535231/149‑1 (MQ=38)
                                                                                  gCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCtttt                                                                   >  3:69101/1‑149 (MQ=255)
                                                                                     gCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTc                                                                >  2:318882/1‑149 (MQ=255)
                                                                                       ggTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCa                                                              <  3:680793/149‑1 (MQ=255)
                                                                                       ggTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCa                                                              <  3:860147/149‑1 (MQ=255)
                                                                                         ttGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTAAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGc                                                            <  2:859479/149‑1 (MQ=255)
                                                                                            tAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGtt                                                         >  2:349948/1‑149 (MQ=255)
                                                                                                tCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCAtt                                                     >  2:37823/1‑149 (MQ=255)
                                                                                                   gATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGc                                                  >  1:651617/1‑149 (MQ=255)
                                                                                                     tAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCtt                                                <  4:87070/149‑1 (MQ=255)
                                                                                                      aaGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTg                                               <  1:783544/149‑1 (MQ=255)
                                                                                                       aGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTgg                                              >  4:696468/1‑149 (MQ=255)
                                                                                                            cgcAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGat                                         <  2:254023/149‑1 (MQ=255)
                                                                                                                    tCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCaa                                 >  1:577701/1‑149 (MQ=255)
                                                                                                                         tCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTtc                            >  1:823867/1‑149 (MQ=255)
                                                                                                                          cAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTtca                           >  2:33860/1‑149 (MQ=255)
                                                                                                                           aGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTtcat                          >  1:350582/1‑149 (MQ=255)
                                                                                                                               gTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATAGTCAATTTTCATCaaa                      >  2:873754/1‑149 (MQ=255)
                                                                                                                                    ttGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATcgc                 <  2:609644/149‑1 (MQ=255)
                                                                                                                                      gCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGa               <  3:52605/149‑1 (MQ=255)
                                                                                                                                         ggATGCGGCGTTGACCCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGtt                                            >  3:646869/1‑117 (MQ=255)
                                                                                                                                         ggATGCGGCGTTGACCCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGtt                                            <  4:646869/117‑1 (MQ=255)
                                                                                                                                            tGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACttttt         <  3:99668/149‑1 (MQ=255)
                                                                                                                                            tGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACttttt         <  1:349948/149‑1 (MQ=255)
                                                                                                                                                gggTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACt     >  3:418672/1‑149 (MQ=255)
                                                                                                                                                   tCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGtt  <  1:318882/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AAATTCGCTGCTGGAGCTTTTTCTGGCTGAGTAAACTTGTATTGCCGGATGCGGCGTGAACGCCTTATCCAGCCGACATGTGGCAGCGGTTGTAGGTCTGATAAGACGCGCAAGCGTCGCATCAGACGTTGATTGCCGGATGCGGGGTCAACGCCGCATCCTGCTACAAATCGTGCACTATATCAAACTTACTTAATCACCCCATTGCGCCGCAGTTGTTCGGTATCTTTTTTCCAGCGTTCATTAGCTTGGTTGATATCGTCAATTTTCATCAAAATCGCGACTTTTTTACTGTT  >  NC_000913/3098370‑3098665

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑