Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,576,390 G→C 8.0% I410M (ATC→ATG)  yddB ← putative TonB‑dependent receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,576,3900GC 8.0% 215.0 / 8.1 71I410M (ATC→ATG) yddBputative TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base G (33/32);  new base C (4/2);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.06e-01

CGATCCGCCATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAATATCACCATA  >  NC_000913/1576247‑1576539
                                                                                                                                               |                                                                                                                                                      
cGATCCGCCATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTAt                                                                                                                                                  >  1:169270/1‑150 (MQ=255)
 gATCCGCCATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTAtt                                                                                                                                                 >  1:629140/1‑150 (MQ=255)
      gCCATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCcgcg                                                                                                                                            >  1:546001/1‑150 (MQ=255)
      gCCATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCcgcg                                                                                                                                            >  1:75889/1‑150 (MQ=255)
        cATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGcc                                                                                                                                          <  2:1363842/150‑1 (MQ=255)
        cATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATACATCTATTCcgatgc                                                                                                                                          <  1:338312/150‑5 (MQ=255)
           ataCAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGaa                                                                                                                                       >  2:1098126/1‑150 (MQ=255)
                 cGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAggg                                                                                                                                 <  2:1177077/150‑1 (MQ=255)
                  gtgtAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGt                                                                                                                                <  1:688800/150‑1 (MQ=255)
                   tgtAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGtt                                                                                                                               >  1:1195212/1‑150 (MQ=255)
                   tgtAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGtt                                                                                                                               >  1:645697/1‑150 (MQ=255)
                     tAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGa                                                                                                                             <  1:1103085/150‑1 (MQ=255)
                            cAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGaaa                                                                                                                      >  1:753583/1‑150 (MQ=255)
                             aaTTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAAc                                                                                                                     <  2:1035171/150‑1 (MQ=255)
                             aaTTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTACCGCAGAATTAATCACATAGGATTCAGACTGGTTATGGCGTTAAGTCCACGCATCGGAATAGATGTATTCAGCGCCGAAGTAGGGTTGATTCGaaaa                                                                                                                     >  1:842894/1‑149 (MQ=255)
                                  ccAGGCGGCCTTTAGCTTTTTGGTAAATGGTAAGGTTAGTTTTCGTTCCGCCAGCATTAAACACATAGGATTCAGACAGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCCCGCCGAAGAAGGGATGATGCGAAACATTAc                                                                                                                >  2:1181040/1‑150 (MQ=255)
                                    aGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTAccc                                                                                                              <  1:1012784/150‑1 (MQ=255)
                                    aGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTAccc                                                                                                              >  1:531318/1‑150 (MQ=255)
                                     ggcggcCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCa                                                                                                             >  2:373207/1‑150 (MQ=255)
                                      ataccatttaccataaAGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCTAc                                                                                                            >  2:1373321/7‑150 (MQ=255)
                                       cggcCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCCGCAGAATTAATCAAATTGGATTCAGACTGGTTATGGATTTCAGTCCACCCATCGGAATACATGTATTCCGCGCCGATGTAGGTTTCATGCGAAAAATTACCCACg                                                                                                           >  2:368028/1‑150 (MQ=255)
                                              tacctttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑tt                                                                                                   >  1:1334415/1‑150 (MQ=255)
                                               acctttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑ttt                                                                                                  >  2:283808/1‑150 (MQ=255)
                                                cctttaTGTTAAATTGTATGGTTATTTTTCTTTCCGGCAGCATTAATCACTTAGGATTCAGACTGGTTATGGAGTTCAGTCCACTCATCTGAATAGATGTATTCCGAGCCGAAGTAGGTTTTATTCTAAACATTACCGACGGCGAA‑tttt                                                                                                 <  2:477412/150‑1 (MQ=255)
                                                 ctttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑ttttt                                                                                                >  2:400289/1‑150 (MQ=255)
                                                 ctttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑ttttt                                                                                                <  2:1110126/150‑1 (MQ=255)
                                                  tttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTg                                                                                               >  1:425381/1‑150 (MQ=255)
                                                  tttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTg                                                                                               <  1:629843/150‑1 (MQ=255)
                                                    taTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGcc                                                                                             <  1:529622/150‑1 (MQ=255)
                                                     aTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCa                                                                                            <  1:846578/150‑1 (MQ=255)
                                                        gTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTACGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTc                                                                                         >  1:1151605/1‑150 (MQ=255)
                                                         tAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTcc                                                                                        >  1:555768/1‑150 (MQ=255)
                                                          aaaTGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCa                                                                                       >  2:961768/1‑150 (MQ=255)
                                                              ggtatggtTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCg                                                                                   >  2:1226972/1‑150 (MQ=255)
                                                                 atggtTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCgtgt                                                                                <  1:1187429/150‑1 (MQ=255)
                                                                   ggtTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGACGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCAAGTCCAGGCGTCttt                                                                              >  2:1065021/1‑150 (MQ=255)
                                                                      tAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGa                                                                           >  2:146120/1‑150 (MQ=255)
                                                                        gTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAg                                                                         >  2:341182/1‑150 (MQ=255)
                                                                         ttttCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAgg                                                                        <  2:1195212/150‑1 (MQ=255)
                                                                           ttctttcCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGgtg                                                                      <  1:373207/150‑1 (MQ=255)
                                                                             ctttcCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTa                                                                    >  2:591460/1‑150 (MQ=255)
                                                                              tttcCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTATGTTTTATGCGAATGATTACACAGTGATAA‑TTTTTGTAATTCAAGGTTTGTTTTGAAGTTGTaa                                                                   <  2:629140/150‑1 (MQ=255)
                                                                                              tCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTggg                                                   <  2:14960/150‑1 (MQ=255)
                                                                                               cacaTAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCTAATTATTAGCCACGGCGAA‑TTTTTGCCAGTACATTCTTGTTTTGAATGTGTAATTATCTACTGCGTGGGa                                                  <  2:645697/150‑1 (MQ=255)
                                                                                                   tAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAAtg                                              >  2:430126/1‑150 (MQ=255)
                                                                                                     ggATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAAtgtg                                            <  2:1334415/150‑1 (MQ=255)
                                                                                                     ggATTCAGACTGGTTATGCCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGAAGAAGTAGCGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAAtgtg                                            >  2:124075/1‑150 (MQ=255)
                                                                                                      gATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATgtgt                                           >  1:1014748/1‑150 (MQ=255)
                                                                                                         tCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGGTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCt                                        <  1:209460/150‑1 (MQ=255)
                                                                                                           aGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTaa                                      <  1:1086646/150‑1 (MQ=255)
                                                                                                              cTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATc                                    >  2:641808/1‑150 (MQ=255)
                                                                                                              cTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATcc                                   >  1:830930/1‑150 (MQ=255)
                                                                                                                ggTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGc                                 <  1:400289/150‑1 (MQ=255)
                                                                                                                     tGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTggg                                                   >  2:379803/1‑127 (MQ=255)
                                                                                                                     tGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTggg                                                   <  1:379803/127‑1 (MQ=255)
                                                                                                                        cGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTg                         <  2:75889/150‑1 (MQ=255)
                                                                                                                          ttCAGTCCACGCATCGGAATAGATGTATTCCGCTCCGAAGTAGGTTTGATGCTAAACTTTACCCACGGCTAA‑TTTTTTCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTTCAGCCTGGGAAATGTGTGCTAATCCGCCACGGGTGCa                       <  2:842894/150‑1 (MQ=255)
                                                                                                                           tCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCaa                      <  2:169270/150‑1 (MQ=255)
                                                                                                                            cAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAAc                     >  2:1050676/1‑150 (MQ=255)
                                                                                                                            aaGTCCACGCATCGGAATAGATGTTTTCCGCGGCGAAGTAGGGTTGATGCGAAACATTACCCATTGCGAA‑TTTTTTCCAGTACAGTCGTGTTTTGAATTTTTAATTATCTACAGCCTGGGAAATGTGTCGTTATCAGCCACGGGTGCAAc                     <  2:425381/149‑1 (MQ=255)
                                                                                                                               tCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGcc                  <  1:146120/150‑1 (MQ=255)
                                                                                                                                ccACTCATCGTAATAGATGTATTCCGCGCCGAAGTATGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCATTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGccc                 <  2:1151605/150‑1 (MQ=255)
                                                                                                                                ccACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGccc                 <  2:546001/150‑1 (MQ=255)
                                                                                                                                  aCGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTg               <  2:663394/150‑1 (MQ=255)
                                                                                                                                   cGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGt              <  1:510151/150‑1 (MQ=255)
                                                                                                                                   cGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGt              <  1:515025/150‑1 (MQ=255)
                                                                                                                                    ggcgcGGAATACCTCTATTCGGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTGGAAGGTGTAATTATCTACAGCCGGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCGTGTa             >  2:1079917/5‑150 (MQ=255)
                                                                                                                                    ggcgcGGAATACATCTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTa             >  1:512314/5‑150 (MQ=255)
                                                                                                                                     gcgcGGAATACATCTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTaa            <  2:72937/147‑1 (MQ=255)
                                                                                                                                     gcgcGGAATACATCTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTaa            >  1:72937/4‑150 (MQ=255)
                                                                                                                                         ggAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAATATc        >  2:1245825/1‑150 (MQ=255)
                                                                                                                                            attgATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCATAATCCGCCACGGGTGCAACGCCCTGTAATATCAcc     >  2:822022/4‑150 (MQ=255)
                                                                                                                                               gATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAATATCACCata  >  2:935871/1‑150 (MQ=255)
                                                                                                                                               |                                                                                                                                                      
CGATCCGCCATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAA‑TTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAATATCACCATA  >  NC_000913/1576247‑1576539

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