Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,039,480 T→C 8.5% V62A (GTC→GCC) ‡ cbdB → cytochrome bd‑II oxidase, subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,039,4800TC 8.5% 292.8 / 12.1 94V62A (GTC→GCC) ‡cbdBcytochrome bd‑II oxidase, subunit II
Reads supporting (aligned to +/- strand):  ref base T (42/44);  new base C (4/4);  total (46/48)
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.02e-01

TGCTGATTGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGAAAAATCGC  >  NC_000913/1039333‑1039627
                                                                                                                                                   |                                                                                                                                                   
tgctgATTGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGt                                                                                                                                                     >  1:389034/1‑148 (MQ=255)
  ctgATTGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCt                                                                                                                                                   >  3:74644/1‑148 (MQ=255)
    gATTGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTgg                                                                                                                                                 <  4:74644/148‑1 (MQ=255)
        ggCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGa                                                                                                                                             <  1:865316/148‑1 (MQ=255)
            tGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCt                                                                                                                                         >  4:901102/1‑148 (MQ=255)
                ccTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCt                                                                                                                                     >  3:397420/1‑148 (MQ=255)
                ccTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCt                                                                                                                                     <  3:901102/148‑1 (MQ=255)
                     tggtCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTg                                                                                                                                                  >  1:798369/1‑130 (MQ=255)
                     tggtCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTg                                                                                                                                                  <  2:121834/130‑1 (MQ=255)
                     tggtCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTg                                                                                                                                                  >  1:121834/1‑130 (MQ=255)
                     tggtCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTg                                                                                                                                                  <  2:798369/130‑1 (MQ=255)
                      ggtCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTggg                                                                                                                               <  4:397420/148‑1 (MQ=255)
                           ttATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACCGCGTTGGTGCACACTGGGAAGGCAACCAGGTc                                                                                                                                                    <  1:320667/122‑1 (MQ=255)
                           ttATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTc                                                                                                                                                    >  2:320667/1‑122 (MQ=255)
                           ttATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGAt                                                                                                                                            >  4:499995/1‑130 (MQ=255)
                           ttATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGAt                                                                                                                                            <  3:499995/130‑1 (MQ=255)
                                aTCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCtgg                                                                                                                                   <  2:732028/134‑1 (MQ=255)
                                aTCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCAGGTTGATCCTCGCtgg                                                                                                                                   >  1:732028/1‑134 (MQ=255)
                                 tctcCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTg                                                                                                                                                  <  4:798578/118‑1 (MQ=255)
                                 tctcCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTg                                                                                                                                                  >  3:798578/1‑118 (MQ=255)
                                    ccGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGACCTGGTTgcctt                                                                                                                                          >  4:294424/1‑119 (MQ=255)
                                    ccGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGACCTGGTTgcctt                                                                                                                                          <  3:294424/123‑5 (MQ=255)
                                            tttGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGACCTGGTTg                                                                                                                                              <  4:565969/111‑1 (MQ=255)
                                            tttGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGACCTGGTTg                                                                                                                                              >  3:565969/1‑111 (MQ=255)
                                             ttGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCag                                                                                                        <  2:389034/148‑1 (MQ=255)
                                               gACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCagag                                                                                                      >  2:662811/1‑148 (MQ=255)
                                                 cATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAgtg                                                                                                    >  1:695286/1‑148 (MQ=255)
                                                   tGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAgtgt                                                                                                   <  3:240212/147‑1 (MQ=255)
                                                   tGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAgtgt                                                                                                   >  4:240212/1‑147 (MQ=255)
                                                   tGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTa                                                                                                  >  3:455510/1‑148 (MQ=255)
                                                     gggATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCt                                                                                                                                         >  1:448706/1‑107 (MQ=255)
                                                     gggATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCt                                                                                                                                         <  2:448706/107‑1 (MQ=255)
                                                        aTCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCtggtg                                                                                                                                 <  3:98896/112‑1 (MQ=255)
                                                        aTCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCtggtg                                                                                                                                 >  4:98896/1‑112 (MQ=255)
                                                           ggCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAgcgg                                                                                          <  1:662811/148‑1 (MQ=255)
                                                              tGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGcc                                                                                                           >  2:124389/1‑128 (MQ=255)
                                                              tGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGcc                                                                                                           <  1:124389/128‑1 (MQ=255)
                                                              tGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGACCTGGTTgc                                                                                                                                             <  1:1030726/94‑2 (MQ=255)
                                                              tGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGACCTGGTTgc                                                                                                                                             >  2:1030726/1‑93 (MQ=255)
                                                               gTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGa                                                                                                                                             <  4:306676/93‑1 (MQ=255)
                                                               gTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGa                                                                                                                                             >  3:306676/1‑93 (MQ=255)
                                                               gTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGtt                                                                                      <  2:695286/148‑1 (MQ=255)
                                                                   gctgcCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGACCTGGTTg                                                                                                                                              <  4:852806/88‑1 (MQ=255)
                                                                   gctgcCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGACCTGGTTg                                                                                                                                              >  3:852806/1‑88 (MQ=255)
                                                                     tgcCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCag                                                                                                        <  4:851065/124‑1 (MQ=255)
                                                                     tgcCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCag                                                                                                        >  3:851065/1‑124 (MQ=255)
                                                                        cGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCtt                                                                             <  4:455510/148‑1 (MQ=255)
                                                                         gCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCttt                                                                            >  1:21940/1‑148 (MQ=255)
                                                                                         gatgatgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTgccgc                                                                                                                 <  1:1010857/95‑1 (MQ=255)
                                                                                         gatgatgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTgccgc                                                                                                                 >  2:1010857/1‑95 (MQ=255)
                                                                                         gatgatgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCTCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGGTGATcc                                                            >  2:346058/1‑148 (MQ=255)
                                                                                             atgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCag                                                                                                        >  3:1131610/1‑100 (MQ=255)
                                                                                             atgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCag                                                                                                        <  4:1131610/100‑1 (MQ=255)
                                                                                              tgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTggg                                                                                                                               <  2:335837/76‑1 (MQ=255)
                                                                                              tgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTggg                                                                                                                               >  1:335837/1‑76 (MQ=255)
                                                                                              tgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTAt                                                                                                 >  4:36698/1‑106 (MQ=255)
                                                                                              tgaACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTAt                                                                                                 <  3:36698/106‑1 (MQ=255)
                                                                                                       gATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTgg                                                                                                                                                 >  4:457944/1‑49 (MQ=255)
                                                                                                       gATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTgg                                                                                                                                                 <  3:457944/49‑1 (MQ=255)
                                                                                                       gATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCa                                              <  1:346058/148‑1 (MQ=255)
                                                                                                            tGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGtt                                         <  2:21940/148‑1 (MQ=255)
                                                                                                             gATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAgtgt                                                                                                   >  1:610053/1‑89 (MQ=255)
                                                                                                             gATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAgtgt                                                                                                   <  2:610053/89‑1 (MQ=255)
                                                                                                              aTAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCt                                                                                                               >  4:864222/1‑76 (MQ=255)
                                                                                                              aTAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCt                                                                                                               <  3:864222/76‑1 (MQ=255)
                                                                                                              aTAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCag                                                                                                        >  1:185012/1‑83 (MQ=255)
                                                                                                              aTAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCag                                                                                                        <  2:185012/83‑1 (MQ=255)
                                                                                                              aTAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGttct                                       >  1:548115/1‑148 (MQ=255)
                                                                                                                 aaCAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTtctc                                                                                   <  2:861529/101‑1 (MQ=255)
                                                                                                                 aaCAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTc                                                                                   >  1:861529/1‑101 (MQ=255)
                                                                                                                  aCAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCg                                                                                                                                       <  3:13971/48‑1 (MQ=255)
                                                                                                                  aCAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCg                                                                                                                                       >  4:13971/1‑48 (MQ=255)
                                                                                                                         tGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTg                            >  1:810548/1‑148 (MQ=255)
                                                                                                                           gTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGc                          <  2:548115/148‑1 (MQ=255)
                                                                                                                                cacTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGgtgct                                                     >  1:690924/1‑116 (MQ=255)
                                                                                                                                cacTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGgtgct                                                     <  2:690924/116‑1 (MQ=255)
                                                                                                                                  cTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGAt                   >  1:839389/1‑148 (MQ=255)
                                                                                                                                      gAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCgg                                                                                <  4:386540/83‑1 (MQ=255)
                                                                                                                                      gAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCgg                                                                                >  3:386540/1‑83 (MQ=255)
                                                                                                                                      gAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTg                     <  1:840840/142‑1 (MQ=255)
                                                                                                                                      gAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTg                     >  2:840840/1‑142 (MQ=255)
                                                                                                                                       aaGTCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCtt                                                                             >  2:630881/1‑85 (MQ=255)
                                                                                                                                       aaGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCtt                                                                             <  1:630881/85‑1 (MQ=255)
                                                                                                                                        aGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATcgc             <  2:810548/148‑1 (MQ=255)
                                                                                                                                         ggCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAgtgt                                                                                                   >  2:199583/1‑61 (MQ=255)
                                                                                                                                         ggCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAgtgt                                                                                                   <  1:199583/61‑1 (MQ=255)
                                                                                                                                          gCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGTGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATcg              <  2:713882/145‑1 (MQ=255)
                                                                                                                                          gCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCa                                              <  1:181703/113‑1 (MQ=255)
                                                                                                                                          gCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCa                                              >  2:181703/1‑113 (MQ=255)
                                                                                                                                          gCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGATGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATcg              >  1:713882/1‑145 (MQ=255)
                                                                                                                                           cAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGgtgctg                                                    <  4:386938/106‑1 (MQ=255)
                                                                                                                                           cAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGgtgctg                                                    >  3:386938/1‑106 (MQ=255)
                                                                                                                                             aCCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAgc                                                                                            <  2:33847/64‑1 (MQ=255)
                                                                                                                                             aCCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAgc                                                                                            >  1:33847/1‑64 (MQ=255)
                                                                                                                                               cAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGaaaaa      <  2:839389/148‑1 (MQ=255)
                                                                                                                                                   tCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGAAAAATCGc  >  3:741385/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGCTGATTGGCGTGATCCTGGTGGTCTTTATGATCTCCGACGGATTTGACATGGGGATCGGCTGTCTGCTGCCGCTGGTGGCGCGTAATGATGATGAACGCCGGATAGTGATAAACAGCGTTGGTGCACACTGGGAAGGCAACCAGGTCTGGTTGATCCTCGCTGGTGGGGCATTATTTGCCGCCTGGCCCAGAGTGTATGCAGCGGCGTTTTCCGGCTTTTATGTGGCGATGATCCTGGTGCTGTGCTCACTGTTCTTCCGCCCGCTGGCCTTTGATTATCGCGGAAAAATCGC  >  NC_000913/1039333‑1039627

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: