Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,576,393 G→C 6.2% 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 6.2% 239.6 / 2.8 73Y409* (TAC→TAGyddBputative TonB‑dependent receptor
Reads supporting (aligned to +/- strand):  ref base G (36/32);  new base C (3/2);  total (39/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGATCCGCCATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAATATCACCATATGT  >  NC_000913/1576247‑1576542
                                                                                                                                                  |                                                                                                                                                      
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)
                                              tacctttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAAtt                                                                                                      >  1:1334415/1‑150 (MQ=255)
                                               acctttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAAttt                                                                                                     >  2:283808/1‑150 (MQ=255)
                                                cctttaTGTTAAATTGTATGGTTATTTTTCTTTCCGGCAGCATTAATCACTTAGGATTCAGACTGGTTATGGAGTTCAGTCCACTCATCTGAATAGATGTATTCCGAGCCGAAGTAGGTTTTATTCTAAACATTACCGACGGCGAAtttt                                                                                                    <  2:477412/150‑1 (MQ=255)
                                                 ctttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAAttttt                                                                                                   >  2:400289/1‑150 (MQ=255)
                                                 ctttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAAttttt                                                                                                   <  2:1110126/150‑1 (MQ=255)
                                                  tttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTg                                                                                                  >  1:425381/1‑150 (MQ=255)
                                                  tttaTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTg                                                                                                  <  1:629843/150‑1 (MQ=255)
                                                    taTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGcc                                                                                                <  1:529622/150‑1 (MQ=255)
                                                     aTGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCa                                                                                               <  1:846578/150‑1 (MQ=255)
                                                        gTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTACGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTc                                                                                            >  1:1151605/1‑150 (MQ=255)
                                                         tAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTcc                                                                                           >  1:555768/1‑150 (MQ=255)
                                                          aaaTGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCa                                                                                          >  2:961768/1‑150 (MQ=255)
                                                              ggtatggtTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCg                                                                                      >  2:1226972/1‑150 (MQ=255)
                                                                 atggtTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCgtgt                                                                                   <  1:1187429/150‑1 (MQ=255)
                                                                   ggtTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGACGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCAAGTCCAGGCGTCttt                                                                                 >  2:1065021/1‑150 (MQ=255)
                                                                      tAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGa                                                                              >  2:146120/1‑150 (MQ=255)
                                                                        gTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAg                                                                            >  2:341182/1‑150 (MQ=255)
                                                                         ttttCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAgg                                                                           <  2:1195212/150‑1 (MQ=255)
                                                                           ttctttcCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGgtg                                                                         <  1:373207/150‑1 (MQ=255)
                                                                             ctttcCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTa                                                                       >  2:591460/1‑150 (MQ=255)
                                                                              tttcCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTATGTTTTATGCGAATGATTACACAGTGATAATTTTTGTAATTCAAGGTTTGTTTTGAAGTTGTaa                                                                      <  2:629140/150‑1 (MQ=255)
                                                                                              tCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTggg                                                      <  2:14960/150‑1 (MQ=255)
                                                                                               cacaTAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCTAATTATTAGCCACGGCGAATTTTTGCCAGTACATTCTTGTTTTGAATGTGTAATTATCTACTGCGTGGGa                                                     <  2:645697/150‑1 (MQ=255)
                                                                                                   tAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAAtg                                                 >  2:430126/1‑150 (MQ=255)
                                                                                                     ggATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAAtgtg                                               <  2:1334415/150‑1 (MQ=255)
                                                                                                     ggATTCAGACTGGTTATGCCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGAAGAAGTAGCGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAAtgtg                                               >  2:124075/1‑150 (MQ=255)
                                                                                                      gATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATgtgt                                              >  1:1014748/1‑150 (MQ=255)
                                                                                                         tCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGGTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCt                                           <  1:209460/150‑1 (MQ=255)
                                                                                                           aGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTaa                                         <  1:1086646/150‑1 (MQ=255)
                                                                                                              cTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATc                                       >  2:641808/1‑150 (MQ=255)
                                                                                                              cTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATcc                                      >  1:830930/1‑150 (MQ=255)
                                                                                                                ggTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGc                                    <  1:400289/150‑1 (MQ=255)
                                                                                                                     tGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTggg                                                      <  1:379803/127‑1 (MQ=255)
                                                                                                                     tGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTggg                                                      >  2:379803/1‑127 (MQ=255)
                                                                                                                        cGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTg                            <  2:75889/150‑1 (MQ=255)
                                                                                                                          ttCAGTCCACGCATCGGAATAGATGTATTCCGCTCCGAAGTAGGTTTGATGCTAAACTTTACCCACGGCTAATTTTTTCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTTCAGCCTGGGAAATGTGTGCTAATCCGCCACGGGTGCa                          <  2:842894/150‑1 (MQ=255)
                                                                                                                           tCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCaa                         <  2:169270/150‑1 (MQ=255)
                                                                                                                            cAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAAc                        >  2:1050676/1‑150 (MQ=255)
                                                                                                                            aaGTCCACGCATCGGAATAGATGTTTTCCGCGGCGAAGTAGGGTTGATGCGAAACATTACCCATTGCGAATTTTTTCCAGTACAGTCGTGTTTTGAATTTTTAATTATCTACAGCCTGGGAAATGTGTCGTTATCAGCCACGGGTGCAAc                        <  2:425381/149‑1 (MQ=255)
                                                                                                                               tCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGcc                     <  1:146120/150‑1 (MQ=255)
                                                                                                                                ccACTCATCGTAATAGATGTATTCCGCGCCGAAGTATGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCATTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGccc                    <  2:1151605/150‑1 (MQ=255)
                                                                                                                                ccACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGccc                    <  2:546001/150‑1 (MQ=255)
                                                                                                                                  aCGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTg                  <  2:663394/150‑1 (MQ=255)
                                                                                                                                   cGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGt                 <  1:510151/150‑1 (MQ=255)
                                                                                                                                   cGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGt                 <  1:515025/150‑1 (MQ=255)
                                                                                                                                    ggcgcGGAATACCTCTATTCGGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTGGAAGGTGTAATTATCTACAGCCGGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCGTGTa                >  2:1079917/5‑150 (MQ=255)
                                                                                                                                    ggcgcGGAATACATCTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTa                >  1:512314/5‑150 (MQ=255)
                                                                                                                                     gcgcGGAATACATCTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTaa               <  2:72937/147‑1 (MQ=255)
                                                                                                                                     gcgcGGAATACATCTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTaa               >  1:72937/4‑150 (MQ=255)
                                                                                                                                         ggAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAATATc           >  2:1245825/1‑150 (MQ=255)
                                                                                                                                            attgATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCATAATCCGCCACGGGTGCAACGCCCTGTAATATCAcc        >  2:822022/4‑150 (MQ=255)
                                                                                                                                               gATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAATATCACCata     >  2:935871/1‑150 (MQ=255)
                                                                                                                                                  gTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAATATCACCATATGt  <  2:555768/150‑1 (MQ=255)
                                                                                                                                                  gTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAATATCACCATATGt  <  1:341182/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
CGATCCGCCATATACAGCGTGTAGTTGTCAATTCCCAGGCGGCCTTTACCTTTATGGTAAATGGTATGGTTAGTTTTCTTTCCGGCAGCATTAATCACATAGGATTCAGACTGGTTATGGCGTTCAGTCCACGCATCGGAATAGATGTATTCCGCGCCGAAGTAGGGTTGATGCGAAACATTACCCACGGCGAATTTTTGCCAGTCCAGGCGTGTTTTGAAGGTGTAATTATCTACAGCCTGGGAAATGTGTCCTAATCCGCCACGGGTGCAACGCCCTGTAATATCACCATATGT  >  NC_000913/1576247‑1576542

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: