Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,576,393 G→C 5.1% Y409* (TAC→TAG yddB ← putative TonB‑dependent receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,576,3930GC 5.1% 263.0 / 4.2 78Y409* (TAC→TAGyddBputative TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base G (41/33);  new base C (2/2);  total (43/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.38e-01

ATGCGATCCGCCATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAATATCACCATA  >  NC_000913/1576244‑1576539
                                                                                                                                                     |                                                                                                                                                  
aTGCGATCCGCCATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATg                                                                                                                                                    >  1:1724785/1‑150 (MQ=255)
 tGCGATCCGCCATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGt                                                                                                                                                   >  1:1576591/1‑150 (MQ=255)
          ccATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATACATCTATTCcgatg                                                                                                                                          <  1:289869/150‑4 (MQ=255)
             tataCAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGa                                                                                                                                       >  2:1227254/1‑150 (MQ=255)
                aCAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGt                                                                                                                                    <  1:1181712/150‑1 (MQ=255)
                   gCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAgg                                                                                                                                 <  1:1936407/150‑1 (MQ=255)
                   gCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAgg                                                                                                                                 >  1:213501/1‑150 (MQ=255)
                   gCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAgg                                                                                                                                 >  1:507771/1‑150 (MQ=255)
                          gTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGAGTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATg                                                                                                                          <  2:401013/150‑1 (MQ=255)
                          gTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATg                                                                                                                          <  2:1576591/150‑1 (MQ=255)
                           ttGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGc                                                                                                                         <  1:513119/150‑1 (MQ=255)
                           ttGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGc                                                                                                                         <  2:1207686/150‑1 (MQ=255)
                           ttGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGc                                                                                                                         >  1:1502490/1‑150 (MQ=255)
                              tCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGaa                                                                                                                      <  2:1480600/150‑1 (MQ=255)
                                  ttCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACAt                                                                                                                  <  1:1739727/150‑1 (MQ=255)
                                    cccAGGCTGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTa                                                                                                                >  2:1073069/1‑150 (MQ=255)
                                    cccAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTa                                                                                                                >  2:1951918/1‑150 (MQ=255)
                                    cccAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTa                                                                                                                >  1:101639/1‑150 (MQ=255)
                                     ccAGGCGGCCTTTACCTTTATGGTAAATGGTATTGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTa                                                                                                                <  1:167748/149‑1 (MQ=255)
                                     ccAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTa                                                                                                                >  2:167748/1‑149 (MQ=255)
                                     ccAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTAc                                                                                                               <  2:123876/150‑1 (MQ=255)
                                                ttacctttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAAt                                                                                                    <  1:882342/150‑1 (MQ=255)
                                                 tacctttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAAtt                                                                                                   <  1:679857/150‑1 (MQ=255)
                                                    ctttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGNGCCGAAGTAGGGTTGATGCGAAACATTAcccn                                                                                                            <  2:1073630/138‑2 (MQ=255)
                                                    ctttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCa                                                                                                            >  1:1073630/1‑138 (MQ=255)
                                                     tttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTg                                                                                               <  2:400314/150‑1 (MQ=255)
                                                     tttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTg                                                                                               >  2:27525/1‑150 (MQ=255)
                                                     tttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTg                                                                                               <  2:1043594/150‑1 (MQ=255)
                                                     tttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTg                                                                                               <  1:34855/150‑1 (MQ=255)
                                                     tttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTg                                                                                               >  1:140739/1‑150 (MQ=255)
                                                       taTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGcc                                                                                             <  1:1165684/150‑1 (MQ=255)
                                                         tGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAg                                                                                           >  2:1475526/1‑150 (MQ=255)
                                                          ggTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGt                                                                                          >  2:1946390/1‑150 (MQ=255)
                                                          ggTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGt                                                                                          >  2:739538/1‑150 (MQ=255)
                                                                aggtatggtTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGc                                                                                    >  2:1640845/2‑150 (MQ=255)
                                                                   tatggtTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCgtg                                                                                 >  1:212463/1‑150 (MQ=255)
                                                                    atggtTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCgtgt                                                                                >  1:1491624/1‑150 (MQ=255)
                                                                      ggtTAGTTATCTTTCCGGCAGCATTAATCACATAGGATACAGACTGGTTATCGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGttt                                                                              >  2:29828/1‑150 (MQ=255)
                                                                            ttttCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAgg                                                                        <  1:1073069/150‑1 (MQ=255)
                                                                                  ttcCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAAt                                                                  >  1:1951068/1‑150 (MQ=255)
                                                                                    ccGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTa                                                                <  1:1227254/150‑1 (MQ=255)
                                                                                    ccGGCAGCATTAATCACATAGGATTCAGACTGGTGATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTa                                                                <  2:1787723/150‑1 (MQ=255)
                                                                                     cGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTAt                                                               >  1:1999363/1‑150 (MQ=255)
                                                                                         agcaTTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTAc                                                           <  2:901989/150‑1 (MQ=255)
                                                                                           caTTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCt                                                             <  2:757337/146‑1 (MQ=255)
                                                                                           caTTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCt                                                             <  2:750558/146‑1 (MQ=255)
                                                                                           caTTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCt                                                             >  1:757337/1‑146 (MQ=255)
                                                                                           caTTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCt                                                             >  1:750558/1‑146 (MQ=255)
                                                                                            aTTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGc                                                        >  1:790697/1‑150 (MQ=255)
                                                                                                aTCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTgg                                                    >  2:496777/1‑150 (MQ=255)
                                                                                                aTCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTgg                                                    >  1:1163832/1‑150 (MQ=255)
                                                                                                 tCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTggg                                                   <  2:1491624/150‑1 (MQ=255)
                                                                                                 tCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGGGTAATTATCTACAGCCTggg                                                   <  2:1119648/150‑1 (MQ=255)
                                                                                                   acaTAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGaa                                                 >  1:683418/1‑150 (MQ=255)
                                                                                                   acaTAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGaa                                                 >  1:679195/1‑150 (MQ=255)
                                                                                                    caTAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGaaa                                                >  1:1695445/1‑150 (MQ=255)
                                                                                                      tAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAAtg                                              <  2:219299/150‑1 (MQ=255)
                                                                                                       aGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCT‑‑GAAATgtgt                                           <  1:641539/150‑1 (MQ=255)
                                                                                                          aTTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTc                                          >  2:1311835/1‑150 (MQ=255)
                                                                                                          aTTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTc                                          <  1:1802306/150‑1 (MQ=255)
                                                                                                          aTTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGACTGGGAAATGTGTc                                          >  1:1105635/1‑150 (MQ=255)
                                                                                                           ttCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTcc                                         >  2:128689/1‑150 (MQ=255)
                                                                                                              aGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTaa                                      >  2:1401020/1‑150 (MQ=255)
                                                                                                              aGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTaa                                      >  1:1170171/1‑150 (MQ=255)
                                                                                                                            gTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGc                        >  1:1715825/1‑150 (MQ=255)
                                                                                                                                aGTCAACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCAACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACg                    >  2:1178075/1‑150 (MQ=255)
                                                                                                                                 gTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGc                   >  1:266344/1‑150 (MQ=255)
                                                                                                                                   ccACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGccc                 <  2:140739/150‑1 (MQ=255)
                                                                                                                                       ggcgcGGAATACATCTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTa             >  2:1063035/5‑150 (MQ=255)
                                                                                                                                        cATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTaa            >  2:676561/1‑150 (MQ=255)
                                                                                                                                           cGGAATACATCTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAatat         >  1:1718335/1‑150 (MQ=255)
                                                                                                                                            ggAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAATATc        >  2:981803/1‑150 (MQ=255)
                                                                                                                                             gAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAATATCa       <  1:1946390/150‑1 (MQ=255)
                                                                                                                                             gAATACATCTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAATATCa       <  1:1063035/150‑1 (MQ=255)
                                                                                                                                                 tgATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCTACGCCCTGTAATATCACCat   <  2:101639/149‑1 (MQ=255)
                                                                                                                                                  gATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCGTAATCCGCCACGGGTGCAACGGCCTGTAATATCACCatg  >  1:899855/1‑149 (MQ=255)
                                                                                                                                                  gATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAATATCACCata  <  2:1999363/150‑1 (MQ=255)
                                                                                                                                                  gATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAATATCACCata  >  2:1914036/1‑150 (MQ=255)
                                                                                                                                                  gATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAATATCACCata  <  1:27525/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                  
ATGCGATCCGCCATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAATATCACCATA  >  NC_000913/1576244‑1576539

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