Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,987,858 G→C 5.7% intergenic (‑100/+15) yecJ ← / ← azuC DUF2766 family protein/acid‑inducible small membrane‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,987,8580GC 5.7% 739.7 / 28.1 210intergenic (‑100/+15)yecJ/azuCDUF2766 family protein/acid‑inducible small membrane‑associated protein
Reads supporting (aligned to +/- strand):  ref base G (104/94);  new base C (3/9);  total (107/103)
Fisher's exact test for biased strand distribution p-value = 7.84e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.78e-01

Only 100 of 214 total aligned reads displayed.

CAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTG  >  NC_000913/1987710‑1988002
                                                                                                                                                     |                                                                                                                                                 
caaCCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑Cg                                                                                                                                                   <  4:1337891/149‑1 (MQ=255)
 aaCCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑Cgg                                                                                                                                                  >  4:1167575/1‑149 (MQ=255)
       aTCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAgg                                                                                                                                            <  1:1671724/149‑1 (MQ=255)
          aGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑TCCGGCAGG‑ttt                                                                                                                                        <  1:1055454/149‑1 (MQ=255)
          aGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑ttt                                                                                                                                        >  4:1053161/1‑149 (MQ=255)
           gAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑TCCGGCAGG‑tttt                                                                                                                                       <  3:419995/149‑1 (MQ=255)
            aaaCTAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑ttttt                                                                                                                                      >  2:360463/1‑149 (MQ=255)
            aaaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑TCCGGCAGG‑ttttt                                                                                                                                      <  2:1542004/149‑1 (MQ=255)
            aaaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑ttttt                                                                                                                                      <  4:1292482/149‑1 (MQ=255)
             aaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑tttttt                                                                                                                                     <  4:1559122/149‑1 (MQ=255)
             aaCCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACAC‑AAAAAACCTGCTCCG‑GCAGG‑ttttttt                                                                                                                                    <  4:1610349/149‑2 (MQ=255)
             aaCCAATTCCTGATCCGCATTCCGCAGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCCTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑tttttt                                                                                                                                     <  2:1141893/149‑1 (MQ=255)
              aCCAATTCGTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTg                                                                                                                                    >  3:1212885/1‑149 (MQ=255)
               ccAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGt                                                                                                                                   >  1:1580098/1‑149 (MQ=255)
               ccAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGt                                                                                                                                   >  2:883404/1‑149 (MQ=255)
               ccAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGt                                                                                                                                   >  1:1199285/1‑149 (MQ=255)
                 aaTTCCTGATCAGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGACGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTa                                                                                                                                 >  1:1432108/1‑149 (MQ=255)
                   ttCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATc                                                                                                                               <  1:691185/149‑1 (MQ=255)
                   ttCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATc                                                                                                                               >  1:820232/1‑149 (MQ=255)
                   ttCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATc                                                                                                                               <  4:1219990/149‑1 (MQ=255)
                      cTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑ttt                                                                                                                                        <  3:1165127/137‑1 (MQ=255)
                       tGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCGGG‑TTTTTTGTTATCGGaa                                                                                                                           <  2:1027767/149‑1 (MQ=255)
                       tGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGaa                                                                                                                           <  2:985654/149‑1 (MQ=255)
                           cGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACata                                                                                                                       <  4:167048/149‑1 (MQ=255)
                                ttCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGcc                                                                                                                  >  1:648024/1‑149 (MQ=255)
                                 tCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCt                                                                                                                 <  3:1479873/149‑1 (MQ=255)
                                         tGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTAc                                                                                                         <  3:1499475/149‑1 (MQ=255)
                                         tGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTAc                                                                                                         <  2:1648758/149‑1 (MQ=255)
                                         tGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTAc                                                                                                         <  4:836156/149‑1 (MQ=255)
                                         tGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTAc                                                                                                         <  2:1391693/149‑1 (MQ=255)
                                         tGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTAc                                                                                                         <  2:147872/149‑1 (MQ=255)
                                          gggACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACg                                                                                                        >  2:408716/1‑149 (MQ=255)
                                          gggACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACg                                                                                                        >  1:1743235/1‑149 (MQ=255)
                                           ggACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGt                                                                                                       >  4:907242/1‑149 (MQ=255)
                                           ggACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTAACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAATACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTActt                                                                                                       <  2:1199285/149‑3 (MQ=255)
                                               aTCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCttt                                                                                                   <  4:1409277/149‑1 (MQ=255)
                                               aTCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCttg                                                                                                   >  1:1551551/1‑148 (MQ=255)
                                                tCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTACAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTg                                                                                                  <  2:425333/149‑1 (MQ=255)
                                                 cGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGa                                                                                                 <  4:1683003/149‑1 (MQ=255)
                                                  gTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑TCCGGCAGG‑tttttt                                                                                                                                     >  4:346904/1‑112 (MQ=255)
                                                  gTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑tttttt                                                                                                                                     <  1:1345719/112‑1 (MQ=255)
                                                  gTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATc                                                                                                                               >  2:1062105/1‑118 (MQ=255)
                                                   tAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAAc                                                                                               <  3:1123838/149‑1 (MQ=255)
                                                        tCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCtt                                                                                          >  1:1169515/1‑149 (MQ=255)
                                                        tCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCtt                                                                                          >  3:923106/1‑149 (MQ=255)
                                                             aTAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAAt                                                                                     >  1:1112768/1‑149 (MQ=255)
                                                                gCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGt                                                                                  <  1:940335/149‑1 (MQ=255)
                                                                   ttttGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTAt                                                                               >  2:123922/1‑149 (MQ=255)
                                                                         aaGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACa                                                                         >  2:1455795/1‑149 (MQ=255)
                                                                          aGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACAt                                                                        >  2:330956/1‑149 (MQ=255)
                                                                            ccGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATAc                                                                      >  4:681114/1‑149 (MQ=255)
                                                                             cGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACt                                                                     >  1:26082/1‑149 (MQ=255)
                                                                              gTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACtt                                                                    >  1:1012554/1‑149 (MQ=255)
                                                                                  tagtagCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCa                                                                <  3:1053161/149‑1 (MQ=255)
                                                                                    gtagCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAgg                                                              >  3:960682/1‑149 (MQ=255)
                                                                                      agCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTAtt                                                                              <  1:596186/131‑1 (MQ=255)
                                                                                       gCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGAtt                                                           >  3:94455/1‑149 (MQ=255)
                                                                                       gCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGAtt                                                           >  2:1314080/1‑149 (MQ=255)
                                                                                         gACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGAtttt                                                         >  2:595865/1‑149 (MQ=255)
                                                                                         gACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGATATTGAACATACTTTTCAGGAtttt                                                         <  4:144374/149‑1 (MQ=255)
                                                                                               acaTCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTg                                                                                         <  2:67904/111‑1 (MQ=255)
                                                                                                caTCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGt                                                  >  4:142981/1‑149 (MQ=255)
                                                                                                 aTCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGtt                                                 <  4:825220/149‑1 (MQ=255)
                                                                                                 aTCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGtt                                                 >  4:558401/1‑149 (MQ=255)
                                                                                                 aTCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGtt                                                 >  3:6200/1‑149 (MQ=255)
                                                                                                   cTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTc                                                                                            <  3:214490/104‑1 (MQ=255)
                                                                                                   cTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTc                                                                                            >  4:214490/1‑104 (MQ=255)
                                                                                                   cTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTc                                               <  1:1517529/149‑1 (MQ=255)
                                                                                                      cAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATc                                            >  3:1443249/1‑149 (MQ=255)
                                                                                                           aTAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGc                                       >  1:67734/1‑149 (MQ=255)
                                                                                                                 acaTCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCg                                 >  3:1262360/1‑149 (MQ=255)
                                                                                                                   aTCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGAc                               >  4:405732/1‑149 (MQ=255)
                                                                                                                    tCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGAcc                              >  2:600283/1‑149 (MQ=255)
                                                                                                                       aCCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCAtt                           <  3:1379293/149‑1 (MQ=255)
                                                                                                                        ccGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCAttt                          <  3:681114/149‑1 (MQ=255)
                                                                                                                         cGTCATCAAGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTg                         <  2:1432108/149‑1 (MQ=255)
                                                                                                                          gTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGt                        <  3:1339592/149‑1 (MQ=255)
                                                                                                                           tcatcaGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGtt                       <  3:907242/149‑1 (MQ=255)
                                                                                                                           tcatcaGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATACTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGtt                       <  4:1057134/149‑1 (MQ=255)
                                                                                                                              tcaGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTtata                    >  3:815358/1‑149 (MQ=255)
                                                                                                                                  gACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAgg                >  1:1325226/1‑149 (MQ=255)
                                                                                                                                   acacAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTTAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGgt               <  2:820232/149‑1 (MQ=255)
                                                                                                                                   acacAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGgt               <  3:1359914/149‑1 (MQ=255)
                                                                                                                                     aacaaaaaaCCTGC‑TCCGGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGAAGCACTCCGACCATTTGTTATACAGgtgt             >  3:487436/4‑149 (MQ=255)
                                                                                                                                      caAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGtt            <  2:1112768/149‑1 (MQ=255)
                                                                                                                                       aaaaaaaCCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTa           >  4:946815/1‑149 (MQ=255)
                                                                                                                                       aaaaaaaCCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTa           <  1:960254/149‑1 (MQ=255)
                                                                                                                                       aaaaaaaCCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTa           >  3:269219/1‑149 (MQ=255)
                                                                                                                                        aaaaaaCCTGCTCCG‑GCAGGTTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTa           >  3:1553599/1‑149 (MQ=255)
                                                                                                                                        aaaaaaCCTGC‑TCCGGCAGGTTTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGtt                       <  3:76910/137‑1 (MQ=255)
                                                                                                                                        aaaaaaCCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAt          <  3:558401/149‑1 (MQ=255)
                                                                                                                                         aaaaaCCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTg                                                                             >  2:490957/1‑81 (MQ=255)
                                                                                                                                         aaaaaCCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTg                                                                             <  1:490957/81‑1 (MQ=255)
                                                                                                                                         aaaaaCCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTAtt         <  2:648024/149‑1 (MQ=255)
                                                                                                                                         aaaaaCCTGC‑CGGAGCAGC‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTTATCGCGGCACTCCGACCATTTGt                        >  4:134853/1‑134 (MQ=255)
                                                                                                                                          aaaaCCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTg        >  1:1134470/1‑149 (MQ=255)
                                                                                                                                               cTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCttt   <  3:144984/149‑1 (MQ=255)
                                                                                                                                               cTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTAATTGTCttt   <  4:1485877/149‑1 (MQ=255)
                                                                                                                                                tGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATATTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTg  >  2:289606/1‑149 (MQ=255)
                                                                                                                                                tGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTg  <  1:422369/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CAACCACATCAGAAACCAATTCCTGATCGGCATTCAGCGGCTGGGACATCGTAAACTCCTGATAGCATTTTGAAAGCCGTTATAGTAGCGACTTCACATCTTCAGCGATAGTCACATCCACCGTCATCAGGACACAAAAAAACCTGC‑CGGAGCAGG‑TTTTTTGTTATCGGAACATATTGCCTGGCGGTACGTCTTTGAACGTCTTGCAATAGTTATTGAACATACTTTTCAGGATTTTGCGCAGTTTCATCGCGGCACTCCGACCATTTGTTATACAGGTGTTATTGTCTTTG  >  NC_000913/1987710‑1988002

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