Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,805,670 A→T 5.1% V741D (GTC→GAC)  PP_4234 ← bifunctional aldehyde oxidase/xanthine dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,805,6700AT 5.1% 246.7 / 2.4 79V741D (GTC→GAC) PP_4234bifunctional aldehyde oxidase/xanthine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (40/35);  new base T (2/2);  total (42/37)
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.14e-01

AGCTGGGTGTCGGGGTAGGCCGCATTCTGCCGGTCGAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAATCGGCGCCAGGCCGTTGCCGGTCAGCTCACCGTTGACCCACACC  >  NC_002947/4805523‑4805816
                                                                                                                                                   |                                                                                                                                                   
aGCTGGGTGTCGGGGTAGGCCGCATTCTGCCGGTCGAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGa                                                                                                                                                     <  6:181675/148‑1 (MQ=255)
 gCTGGGTGTCGGGGTAGGCCGCATTCTGCCGGTCGAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACa                                                                                                                                                   <  7:65808/149‑1 (MQ=255)
  cTGGGTGTCGGGGTAGGCCGCATTCTGCCGGTCGAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACa                                                                                                                                                   >  4:96129/1‑148 (MQ=255)
       tgtCGGGGTAGGCCGCATTCTGCCGGTCGAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGtcg                                                                                                                                             >  1:239891/1‑149 (MQ=255)
            gggTAGGCCGCATTCTGCCGGTCGAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGtcggcctc                                                                                                                                        >  2:296723/1‑149 (MQ=255)
              gTAGGCCGCATTCTGCCGGTCGAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGtcggcctcgg                                                                                                                                      <  2:58473/149‑1 (MQ=255)
              gTAGGCCGCATTCTGCCGGTCGAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGACATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGtcggcctcgg                                                                                                                                      >  7:248140/1‑149 (MQ=255)
                   ccGCATTCTGCCGGTCGAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAg                                                                                                                                 <  5:332883/149‑1 (MQ=255)
                    cGCATTCTGCCGGTCGAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAg                                                                                                                                 >  3:127314/1‑148 (MQ=255)
                    cGCATTCTGCCGGTCGAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGc                                                                                                                                >  2:121266/1‑149 (MQ=255)
                      cATTCTGCCGGTCGAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTc                                                                                                                              >  5:135617/1‑149 (MQ=255)
                         tCTGCCGGTCGAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAg                                                                                                                          >  8:315198/1‑149 (MQ=255)
                             ccGGTCGAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGt                                                                                                                      <  5:84889/149‑1 (MQ=255)
                             ccGGTCGAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGt                                                                                                                      >  8:213523/1‑149 (MQ=255)
                                 tCGAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGaa                                                                                                                  <  7:287167/149‑1 (MQ=255)
                                 tCGAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGaa                                                                                                                  <  7:315198/149‑1 (MQ=255)
                                   gAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAg                                                                                                                 <  7:213523/148‑1 (MQ=255)
                                   gAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGc                                                                                                                >  5:340057/1‑149 (MQ=255)
                                    aGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGcc                                                                                                               <  2:214686/149‑1 (MQ=255)
                                    aGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGcc                                                                                                               >  8:167601/1‑149 (MQ=255)
                                        ggTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGc                                                                                                           >  6:96023/1‑149 (MQ=255)
                                           gATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACaaa                                                                                                                                                 >  6:254648/1‑109 (MQ=255)
                                           gATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACaaa                                                                                                                                                 <  5:254648/109‑1 (MQ=255)
                                               gCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCa                                                                                                    <  6:135617/149‑1 (MQ=255)
                                                ccATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGTCAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCAc                                                                                                   >  7:326524/1‑149 (MQ=255)
                                                  aTGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGc                                                                                                 <  3:328704/149‑1 (MQ=255)
                                                     gCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGTCAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCg                                                                                               <  8:326524/148‑1 (MQ=255)
                                                        gCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTaa                                                                                           >  8:30926/1‑149 (MQ=255)
                                                         cccATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTa                                                                                            <  3:96129/147‑1 (MQ=255)
                                                         cccATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAAc                                                                                          <  7:237824/149‑1 (MQ=255)
                                                          ccATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTAACGCCGGTAAc                                                                                          >  1:203586/1‑148 (MQ=255)
                                                           cATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGt                                                                                                                      <  3:336702/119‑1 (MQ=255)
                                                           cATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGt                                                                                                                      >  4:336702/1‑119 (MQ=255)
                                                           cATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCa                                                                                        >  6:168867/1‑149 (MQ=255)
                                                           cATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCa                                                                                        <  5:219054/149‑1 (MQ=255)
                                                               tGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGccc                                                                                    >  6:89227/1‑149 (MQ=255)
                                                                gcgcCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGccc                                                                                    >  7:142473/1‑148 (MQ=255)
                                                                 cgcCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGccc                                                                                    <  7:167601/147‑1 (MQ=255)
                                                                 cgcCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGcc                                                                                     <  1:296723/146‑1 (MQ=255)
                                                                    cttcttGGCATTCACCGCACAATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGATCACGCCGGGAACCCGCCCCATGt                                                                               >  2:293584/1‑149 (MQ=255)
                                                                        ttGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTg                                                                           <  6:73203/149‑1 (MQ=255)
                                                                        ttGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTg                                                                           >  5:73203/1‑149 (MQ=255)
                                                                          ggCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGccc                                                                                    >  3:51635/1‑138 (MQ=255)
                                                                          ggCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGccc                                                                                    <  4:51635/138‑1 (MQ=255)
                                                                          ggCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGt                                                                            <  2:239891/146‑1 (MQ=255)
                                                                          ggCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGt                                                                            <  7:288362/146‑1 (MQ=255)
                                                                          ggCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTgg                                                                          <  5:168867/148‑1 (MQ=255)
                                                                           gCATTCACCGCACCATCGCCGTATTTCCCCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGc                                                                        >  7:180558/1‑149 (MQ=255)
                                                                              ttCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCtt                                                                     >  3:200366/1‑149 (MQ=255)
                                                                                 aCCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCt                                                                   <  7:30926/148‑1 (MQ=255)
                                                                                  ccGCACGATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTgc                                                                 >  4:64551/1‑149 (MQ=255)
                                                                                         aTCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCt                                                          >  6:211474/1‑149 (MQ=255)
                                                                                          tCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTg                                                         <  5:89227/149‑1 (MQ=255)
                                                                                           cGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGc                                                                                                 <  1:291885/108‑1 (MQ=255)
                                                                                           cGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGc                                                                                                 >  2:291885/1‑108 (MQ=255)
                                                                                             ccGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTggg                                                       >  3:98358/1‑148 (MQ=255)
                                                                                              cGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGc                                                      >  7:201553/1‑148 (MQ=255)
                                                                                                  tttCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACg                                                 <  6:340057/149‑1 (MQ=255)
                                                                                                  attCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAAc                                                  <  3:64551/147‑1 (MQ=255)
                                                                                                      aCCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAAt                                             >  2:242282/1‑149 (MQ=255)
                                                                                                       ccGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAATc                                            >  1:161969/1‑149 (MQ=255)
                                                                                                           cATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCCTGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAATCGgcg                                        >  2:312753/1‑149 (MQ=255)
                                                                                                            aTGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAATCGgcgc                                       <  4:200366/149‑1 (MQ=255)
                                                                                                                cATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAATCGGCGCCAgg                                   >  7:235613/1‑149 (MQ=255)
                                                                                                                    gAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAATCGGCGCCAGGCCGt                               >  3:243657/1‑149 (MQ=255)
                                                                                                                      gCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAATCGGCGCCAGGCCGTTg                             >  3:101363/1‑149 (MQ=255)
                                                                                                                      gCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAATCGGCGCCAGGCCGTTg                             <  8:248140/149‑1 (MQ=255)
                                                                                                                       cGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCa                                                     <  6:241452/124‑1 (MQ=255)
                                                                                                                       cGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCa                                                     >  5:241452/1‑124 (MQ=255)
                                                                                                                         ggAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAATCGGCGCCAGGCCGTTGCCg                          >  1:229697/1‑149 (MQ=255)
                                                                                                                              cgcTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAATCGGCGCCAGGCCGTTGCCGGTCAg                     <  8:201553/149‑1 (MQ=255)
                                                                                                                               gcTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCGGAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGt                                                                            >  1:264749/1‑94 (MQ=255)
                                                                                                                               gcTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGt                                                                            <  2:264749/93‑1 (MQ=255)
                                                                                                                                cTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGGAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAATCGGCGCCAGGCCGTTGCCGGTCAGCt                   <  8:142473/149‑1 (MQ=255)
                                                                                                                                    cgcACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAATCGGCGCCAGGCCGTTGCCGGTCAGCTCAcc               <  4:98358/149‑1 (MQ=255)
                                                                                                                                           cATCGATGTCAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAATCGGCGCCAGGCCGTTg                             <  8:274350/128‑1 (MQ=255)
                                                                                                                                           cATCGATGTCAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAATCGGCGCCAGGCCGTTg                             >  7:274350/1‑128 (MQ=255)
                                                                                                                                           cATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAATCGGCGCCAGGCCGTTGCCGGTCAGCTCACCGTTGAcc        <  4:101363/149‑1 (MQ=255)
                                                                                                                                                aTGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAATCGGCGCCAGGCCGTTGCCGGTCAGCTCACCGTTGACCCacac   >  3:73445/1‑149 (MQ=255)
                                                                                                                                                 tGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAATCGGCGCCAGGCCGTTGCCGGTCAGCTCACCGTTGACCCACAcc  >  1:222281/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGCTGGGTGTCGGGGTAGGCCGCATTCTGCCGGTCGAGCAGGTGATAGCCATGGCTGCCCATCTGCGCCTTCTTGGCATTCACCGCACCATCGCCGTATTTCACCGCCATGGCATCGAGCGGGAAGCGCTCGCGCACGCCATCGATGACAAAGTCGGCCTCGGCCCAGCTCCAGCGGTTGAAGCCGTGCACGCTCACGCCGGTAACCAGCCCCATGTCGTGGGCCTTCTGCGCCAGCTGGGCAAACGGAATCGGCGCCAGGCCGTTGCCGGTCAGCTCACCGTTGACCCACACC  >  NC_002947/4805523‑4805816

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