Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 478,759 A→G 15.9% intergenic (‑142/+22) 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,7590AG15.9% 291.8 / 43.3 125intergenic (‑142/+22)ylaB/ylaCputative membrane‑anchored cyclic‑di‑GMP phosphodiesterase/DUF1449 family inner membrane protein
Reads supporting (aligned to +/- strand):  ref base A (60/45);  new base G (11/9);  total (71/54)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.23e-01

Only 100 of 126 total aligned reads displayed.

TTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTT  >  NC_000913/478660‑478857
                                                                                                    |                                                                                                  
ttctgGAAAAATGTTGGCTGGGGTGGAAGAGGTAGAGCTGGTATTACCTTAAGTAATGGGCTGGAAA‑AATAATTGTCATTTGACCAGGCTGAGAGCTTTa                                                                                                    <  1:1818007/97‑1 (MQ=255)
  gtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                                  >  1:1820555/1‑100 (MQ=255)
  gtgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACg                                                                                                  >  1:2714565/1‑100 (MQ=255)
   tgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGc                                                                                                 <  1:553942/100‑1 (MQ=255)
   tgGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGc                                                                                                 >  2:517406/1‑100 (MQ=255)
    ggAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGcc                                                                                                >  1:2391189/1‑100 (MQ=255)
      aaagaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACTATCATTGTCATTTGACCAGGCGCATAGCTTTACGCCTg                                                                                              <  1:1801747/100‑1 (MQ=255)
       aaaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTgg                                                                                             <  1:1292628/100‑1 (MQ=255)
        aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGt                                                                                            <  2:1080633/100‑1 (MQ=38)
        aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGt                                                                                            <  2:2222516/100‑1 (MQ=38)
        aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGt                                                                                            <  2:2522339/100‑1 (MQ=255)
        aaaGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGt                                                                                            >  1:2522324/1‑100 (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)
          aGGCTGGGTGGAGATAAGGAGGTATAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTa                                                                                           <  1:2480151/99‑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)
            gctggctgGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtcaa                                                                                         <  1:791572/99‑4 (MQ=38)
              tggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGTAAAGCTCTGCGCCTGGtcaaatg                                                                                      <  2:1718990/100‑7 (MQ=37)
               ggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTaa                                                                                      >  1:2253560/1‑99 (MQ=255)
               ggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                                     >  1:1543999/1‑100 (MQ=255)
               ggctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAc                                                                                     <  1:1610533/100‑1 (MQ=255)
                gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAcc                                                                                    >  2:2204370/1‑100 (MQ=255)
                gctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAcc                                                                                    <  1:106454/100‑1 (MQ=255)
                 ctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                   >  1:1477357/1‑100 (MQ=255)
                 ctggAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAAccc                                                                                   >  1:1234347/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)
                         aaGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTaa                                                                           <  2:2695264/100‑1 (MQ=255)
                           ggaggTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGa                                                                         <  2:1935661/100‑1 (MQ=255)
                              ggTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGt                                                                      >  2:391785/1‑100 (MQ=255)
                               gTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGtt                                                                     <  2:1512195/100‑1 (MQ=255)
                               gTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGtt                                                                     <  2:372215/100‑1 (MQ=255)
                                tAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTg                                                                    <  2:2911991/100‑1 (MQ=255)
                                tAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTg                                                                    >  1:1756862/1‑100 (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)
                                            ttCATGAAGTGATGGG‑TGGAAACGATCATTGTCCTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgag                                                        >  1:1210496/1‑100 (MQ=255)
                                            ttCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgag                                                        >  2:487558/1‑100 (MQ=255)
                                             tCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCgcgagc                                                       >  2:2103699/1‑100 (MQ=255)
                                                tGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGcgagcgag                                                    >  2:783039/1‑100 (MQ=255)
                                                tGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGcgagcgag                                                    <  2:2473225/100‑1 (MQ=255)
                                                 tccgatctGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                   <  2:59675/93‑1 (MQ=255)
                                                 gAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                   >  2:2298456/1‑100 (MQ=255)
                                                 gAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAgg                                                   <  2:82484/100‑1 (MQ=255)
                                                   aGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGt                                                 >  2:285502/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)
                                                         ggg‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAggagatcgga                                           >  1:59674/1‑92 (MQ=255)
                                                             tGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                        <  2:2824314/100‑1 (MQ=255)
                                                              ggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATc                                        >  2:1175295/1‑99 (MQ=255)
                                                              ggAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCg                                       <  2:1429652/100‑1 (MQ=255)
                                                                   cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaa                                  >  1:1715004/1‑100 (MQ=255)
                                                                   cGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaa                                  >  1:405760/1‑100 (MQ=255)
                                                                    gATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTaaaaa                                 <  2:1668375/100‑1 (MQ=255)
                                                                     aTCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAg                                >  1:663335/1‑100 (MQ=255)
                                                                     aTCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAg                                <  2:547494/100‑1 (MQ=255)
                                                                     aTCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAg                                <  1:1790257/100‑1 (MQ=255)
                                                                       cATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGAc                              <  1:1386310/100‑1 (MQ=255)
                                                                        aTTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACa                             <  2:2548068/100‑1 (MQ=255)
                                                                          tGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAg                            <  1:2627614/99‑1 (MQ=255)
                                                                          tGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACGTAAGAGGTGGATGCGGCGCGAGCGAGGGTGAAAATATGGTAGAAATACAGt                           >  2:2323073/1‑100 (MQ=255)
                                                                          ggggtacataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGt                           >  2:701596/10‑100 (MQ=25)
                                                                            gttacataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTc                         >  2:338508/8‑100 (MQ=25)
                                                                             cATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCa                        <  2:1216456/100‑1 (MQ=255)
                                                                              aTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAc                       >  1:1655214/1‑100 (MQ=255)
                                                                              aTTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCAc                       >  2:2890209/1‑100 (MQ=255)
                                                                                 ataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt                    >  2:1310070/3‑100 (MQ=38)
                                                                                 ataCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCtt                    >  2:1834782/3‑100 (MQ=38)
                                                                                  taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                   >  1:790851/2‑100 (MQ=38)
                                                                                  taCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                   >  2:2526989/2‑100 (MQ=38)
                                                                                  taCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCttt                   >  2:296245/2‑100 (MQ=255)
                                                                                   aCCAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                  >  2:81918/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)
                                                                                   aCCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTa                  <  1:1377260/100‑1 (MQ=255)
                                                                                    ccAGGCGTAAAGCTCTGCGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                 >  2:269199/1‑100 (MQ=255)
                                                                                    ccAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTAc                 <  1:1034698/100‑1 (MQ=255)
                                                                                     cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAACAATATCGTAAAAAGTCAGTTCACGTTTACg                >  1:2053055/1‑100 (MQ=255)
                                                                                     cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg                >  1:1660480/1‑100 (MQ=255)
                                                                                     cAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACg                >  1:2898621/1‑100 (MQ=255)
                                                                                       ggCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa              <  2:1064095/100‑1 (MQ=255)
                                                                                       ggCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGa              >  2:2398806/1‑100 (MQ=255)
                                                                                         cgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg            <  2:610801/100‑1 (MQ=255)
                                                                                         cgcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACg            <  1:2607997/100‑1 (MQ=255)
                                                                                          gcAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGa           >  2:2772212/1‑100 (MQ=255)
                                                                                             gagCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCTAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGCACGATc         >  1:656015/1‑99 (MQ=255)
                                                                                               acTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTGACCTTTACGGACGATCAtt      >  1:481447/2‑100 (MQ=255)
                                                                                                cTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttt     >  1:2502268/1‑100 (MQ=255)
                                                                                                 tttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAtttt    >  2:55931/1‑100 (MQ=255)
                                                                                                  ttACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCAttttt   <  2:33168/100‑1 (MQ=255)
                                                                                                    |                                                                                                  
TTGTGGAAAAAGGCTGGCTGGAGATAAGGAGGTAGAGCTGGTATTTCATGAAGTGATGGG‑TGGAAACGATCATTGTCATTTGACCAGGCGCAGAGCTTTACGCCTGGTATGTAACCCTAACTTAAGATGTTGATTCGGCGCGAGCGAGGGTAAAAATATCGTAAAAAGACAGTTCACCTTTACGGACGATCATTTTTT  >  NC_000913/478660‑478857

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: ‑