Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 478,750 C→T 16.5% intergenic (‑133/+31) ylaB ← / ← ylaC putative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913478,7500CT16.5% 281.8 / 42.7 124intergenic (‑133/+31)ylaB/ylaCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (57/45);  new base T (11/9);  total (69/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.37e-01

Only 100 of 129 total aligned reads displayed.

TATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGA  >  NC_000913/478651‑478848
                                                                                                    |                                                                                                  
gatatCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGGATTTCATGAAGTTATGGG‑TGGAAACGATCATTGTCATTTGACCAGgcgc                                                                                                    >  1:2768761/2‑100 (MQ=255)
 atatCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCa                                                                                                   <  2:2655533/100‑1 (MQ=255)
   atCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCagn                                                                                                 <  1:2407014/100‑2 (MQ=255)
   atCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCaga                                                                                                 >  2:2244242/1‑100 (MQ=255)
     cTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGc                                                                                               >  1:1268458/1‑100 (MQ=255)
     cTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGc                                                                                               <  2:260076/100‑1 (MQ=255)
       gATTGTGGAAAAAGGCTGGCTGGAGCTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCtt                                                                                             >  1:2368595/1‑100 (MQ=255)
       gATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCtt                                                                                             <  2:803703/100‑1 (MQ=255)
       gATTGTGAAAAAAGGCTGGCCGAAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAAAGACAAATGTCATTCGAACACACGCagag                                                                                                >  1:1141952/1‑97 (MQ=255)
         ttctgGAAAAATGTTGGCTGGGGTGGAAGAGGTAGAGCTGGTATTACCTTAAGTAATGGGCTGGAAA‑AATAATTGTCATTTGACCAGGCTGAGAGCTTTa                                                                                           <  1:1818007/97‑1 (MQ=255)
           gtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                         >  1:2714565/1‑100 (MQ=255)
           gtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                         >  1:1820555/1‑100 (MQ=255)
           gtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                         >  2:553945/1‑100 (MQ=255)
            tgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGc                                                                                        <  1:553942/100‑1 (MQ=255)
            tgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGc                                                                                        >  2:517406/1‑100 (MQ=255)
             ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                       >  2:1741436/1‑100 (MQ=255)
             ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                       >  1:2391189/1‑100 (MQ=255)
              gAAAAAGGCTGGCTGGAGATAAGGAGGGGGAGCTGGTGTTTCATGAAGTAATGGG‑AAGAAACGATCATAGTTATTTGACCAGGCGAAGAGCTTTACGCCt                                                                                      >  1:877725/1‑100 (MQ=255)
               aaagaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACTATCATTGTCATTTGACCAGGCGCATAGCTTTACGCCTg                                                                                     <  1:1801747/100‑1 (MQ=255)
                aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTgg                                                                                    <  1:1292628/100‑1 (MQ=255)
                 aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGt                                                                                   <  2:2222516/100‑1 (MQ=38)
                 aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGt                                                                                   >  1:2522324/1‑100 (MQ=255)
                 aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGt                                                                                   <  2:2522339/100‑1 (MQ=255)
                  aaGGTTTGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGt                                                                                   <  1:873004/99‑1 (MQ=255)
                  aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtc                                                                                  <  2:164933/100‑2 (MQ=38)
                  aaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTa                                                                                  <  2:2220290/100‑1 (MQ=255)
                   aGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtca                                                                                 <  1:1234385/100‑3 (MQ=38)
                   aGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTAt                                                                                 >  2:1342554/1‑100 (MQ=255)
                    ggctggctGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtcaa                                                                                <  2:1712979/100‑4 (MQ=38)
                     gctggctgGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTTACCAGGCGTAAAGCTCTGCGCCTGGtcaaa                                                                               <  1:1913926/100‑5 (MQ=37)
                       tggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtcaaatg                                                                             <  2:1718990/100‑7 (MQ=37)
                        ggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTaa                                                                             >  1:2253560/1‑99 (MQ=255)
                        ggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                            <  1:1610533/100‑1 (MQ=255)
                         gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAcc                                                                           <  1:106454/100‑1 (MQ=255)
                         gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAcc                                                                           >  2:2204370/1‑100 (MQ=255)
                          ctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                          >  1:1477357/1‑100 (MQ=255)
                            ggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTa                                                                        >  2:1087716/1‑100 (MQ=255)
                            ggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTa                                                                        >  1:1064091/1‑100 (MQ=255)
                            ggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTa                                                                        >  1:33167/1‑100 (MQ=255)
                              agaTAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCGTTACGCCTGGTATGTAACCCTAAc                                                                      <  1:2456383/100‑1 (MQ=255)
                                 tAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTa                                                                   <  1:2509328/100‑1 (MQ=255)
                                  aaGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTaa                                                                  <  2:2695264/100‑1 (MQ=255)
                                    ggaggTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGa                                                                <  2:1935661/100‑1 (MQ=255)
                                        gTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGtt                                                            <  2:372215/100‑1 (MQ=255)
                                        gTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGtt                                                            <  2:1512195/100‑1 (MQ=255)
                                         tAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTg                                                           <  2:2911991/100‑1 (MQ=255)
                                         tAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTg                                                           >  1:1756862/1‑100 (MQ=255)
                                          agagCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGa                                                          <  1:2840523/100‑1 (MQ=255)
                                              cTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCg                                                      <  2:2689090/100‑1 (MQ=255)
                                                ggTATTTCATGAAGTGATGGG‑TGGCAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgc                                                    >  1:105171/1‑100 (MQ=255)
                                                ggTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgc                                                    >  2:1388287/1‑100 (MQ=255)
                                                ggTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgc                                                    >  2:1386315/1‑100 (MQ=255)
                                                  tATTTCATGAAGTGATGGG‑TGTAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGgcgc                                                  <  2:1094417/100‑1 (MQ=255)
                                                     ttCATGAAGTGATGGG‑TGGAAACGATCATTGTCCTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgag                                               >  1:1210496/1‑100 (MQ=255)
                                                     ttCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgag                                               >  2:487558/1‑100 (MQ=255)
                                                      tCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagc                                              >  2:2103699/1‑100 (MQ=255)
                                                          gAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                          <  2:82484/100‑1 (MQ=255)
                                                          gAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                          >  2:2298456/1‑100 (MQ=255)
                                                            aGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGt                                        >  2:285502/1‑100 (MQ=255)
                                                            aGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGt                                        >  1:1658921/1‑100 (MQ=255)
                                                             gTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTa                                       >  1:1111128/1‑100 (MQ=255)
                                                               gATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaa                                     >  1:2060274/1‑100 (MQ=255)
                                                                aTGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTaaaa                                    <  2:1821687/100‑1 (MQ=255)
                                                                    g‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAatat                                <  2:1544005/100‑1 (MQ=255)
                                                                       ggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                               >  2:1175295/1‑99 (MQ=255)
                                                                       ggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                              >  2:2597603/1‑100 (MQ=255)
                                                                           aCGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaa                          >  2:2786876/1‑100 (MQ=255)
                                                                            cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaa                         >  1:405760/1‑100 (MQ=255)
                                                                            cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaa                         >  1:1715004/1‑100 (MQ=255)
                                                                              aTCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAg                       >  1:663335/1‑100 (MQ=255)
                                                                              aTCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAg                       <  2:547494/100‑1 (MQ=255)
                                                                              aTCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAg                       <  1:1790257/100‑1 (MQ=255)
                                                                               tCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGa                      >  2:2916933/1‑100 (MQ=255)
                                                                                cATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGAc                     <  1:1386310/100‑1 (MQ=255)
                                                                                 aTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACa                    >  1:1512019/1‑100 (MQ=255)
                                                                                   tGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACGTAAGAGGTGGATGCGGCGCGAGCGAGGGTGAAAATATGGTAGAAATACAGt                  >  2:2323073/1‑100 (MQ=255)
                                                                                   ggggtacataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGt                  >  2:701596/10‑100 (MQ=25)
                                                                                     tCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTc                <  1:789061/100‑1 (MQ=255)
                                                                                     gttacataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTc                >  2:338508/8‑100 (MQ=25)
                                                                                       aTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAc              >  2:2890209/1‑100 (MQ=255)
                                                                                       aTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAc              >  1:1655214/1‑100 (MQ=255)
                                                                                       aTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAc              <  1:535672/100‑1 (MQ=255)
                                                                                          ataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt           >  2:1310070/3‑100 (MQ=38)
                                                                                          ataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt           >  2:779315/3‑100 (MQ=38)
                                                                                           taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt          >  1:790851/2‑100 (MQ=38)
                                                                                           taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt          >  2:2526989/2‑100 (MQ=38)
                                                                                            aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa         >  2:81918/1‑100 (MQ=38)
                                                                                            aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa         >  2:3852/1‑100 (MQ=38)
                                                                                            aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAAACCTAACTTAAGATGTTGATTCGGCGCGAGTGAGGGTAAAAATATCGTAAAAGGACAGTTCACCttt          >  1:877412/1‑99 (MQ=38)
                                                                                            aCCAGGCGCAGTGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGAAGCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt          <  1:1615329/99‑1 (MQ=255)
                                                                                            aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa         <  2:212528/100‑1 (MQ=255)
                                                                                             ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc        <  1:1034698/100‑1 (MQ=255)
                                                                                              cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAACAATATCGTAAAAAGTCAGTTCACGTTTACg       >  1:2053055/1‑100 (MQ=255)
                                                                                               aGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACgg      >  2:2504580/1‑100 (MQ=255)
                                                                                                ggCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa     <  2:1064095/100‑1 (MQ=255)
                                                                                                ggCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa     >  2:2398806/1‑100 (MQ=255)
                                                                                                  cgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg   <  1:2607997/100‑1 (MQ=255)
                                                                                                  cgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg   <  1:1283158/100‑1 (MQ=255)
                                                                                                  cgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg   <  2:610801/100‑1 (MQ=255)
                                                                                                   gcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGa  >  2:2772212/1‑100 (MQ=255)
                                                                                                    |                                                                                                  
TATATCTGATTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGA  >  NC_000913/478651‑478848

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