Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 138,664 T→G 8.4% intergenic (+31/+171) cueO → / ← gcd multicopper oxidase (laccase)/glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913138,6640TG 8.4% 305.9 / 15.9 107intergenic (+31/+171)cueO/gcdmulticopper oxidase (laccase)/glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (47/51);  new base G (6/3);  total (53/54)
Fisher's exact test for biased strand distribution p-value = 3.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

Only 100 of 117 total aligned reads displayed.

ATGTCAGCGAAGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACG  >  NC_000913/138515‑138812
                                                                                                                                                     |                                                                                                                                                    
aTGTCAGCGAAGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGAtct                                                                                                                                                      <  2:464786/150‑1 (MQ=255)
  gTCAGCGAAGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGAtctct                                                                                                                                                    >  1:94795/1‑150 (MQ=255)
  gTCAGCGAAGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGAtctct                                                                                                                                                    >  2:855658/1‑150 (MQ=255)
  gTCAGCGAAGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGAtctct                                                                                                                                                    >  2:1593677/1‑150 (MQ=255)
  gTCAGCGAAGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGAtctct                                                                                                                                                    >  1:1084293/1‑150 (MQ=255)
  gTCAGCGAAGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGAtctct                                                                                                                                                    <  2:538456/150‑1 (MQ=255)
  gTCAGCGAAGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTAAATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATACCCGGAtctct                                                                                                                                                    >  1:1544835/1‑150 (MQ=255)
   tCAGCGAAGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGCGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGAtctctc                                                                                                                                                   >  2:201999/1‑150 (MQ=255)
         aaGTGCTGGTGAATTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCggg                                                                                                                                               <  2:1293267/148‑1 (MQ=255)
         aaGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCggg                                                                                                                                               >  1:1293267/1‑148 (MQ=255)
          aGTGCTGGTTCAGTTTAATCACGCTGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGAGAGATCCGGGCAt                                                                                                                                            <  2:292325/150‑1 (MQ=255)
          aGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGAGAGATCCGGGCAt                                                                                                                                            <  2:1258228/150‑1 (MQ=255)
           gTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGAGAGATCCGGGCAta                                                                                                                                           <  1:61466/150‑2 (MQ=255)
               tGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAAc                                                                                                                                       <  1:380731/150‑1 (MQ=255)
               tGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAAc                                                                                                                                       >  1:1536591/1‑150 (MQ=255)
               tGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAAc                                                                                                                                       >  2:736497/1‑150 (MQ=255)
                   gAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACtgat                                                                                                                                   >  1:507273/1‑150 (MQ=255)
                   gAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACtgat                                                                                                                                   >  1:1592361/1‑150 (MQ=255)
                    taGTTTAATCACGATGCACCGAAAGTACATGCTTATATGGCTCACTGCCATCTGCTGGAGCATGAAGATACGGGTATTTTGTTAGGTTTTACTGTATAAAAAGACCGTATTCGGAAATATGCCCGGTTCTCTCCGGACATTAAACtgatg                                                                                                                                  <  2:1316442/149‑1 (MQ=255)
                    aaGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACtgatg                                                                                                                                  <  2:407685/150‑1 (MQ=255)
                      gTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGAc                                                                                                                                >  2:393864/1‑150 (MQ=255)
                      gTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGAc                                                                                                                                <  1:971446/150‑1 (MQ=255)
                      gTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGAc                                                                                                                                >  2:176453/1‑150 (MQ=255)
                        ttAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACaa                                                                                                                              >  2:14737/1‑150 (MQ=255)
                         tAATCACGATGGTCCGAAAGAACATGCTTATTTGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGTTGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACaaa                                                                                                                             <  1:163142/150‑1 (MQ=255)
                          aaTCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAAc                                                                                                                            <  1:1083181/150‑1 (MQ=255)
                           aTCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACg                                                                                                                           <  2:1084293/150‑1 (MQ=255)
                            tCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGc                                                                                                                          >  1:1537082/1‑150 (MQ=255)
                                   gCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTg                                                                                                                   <  1:44535/150‑1 (MQ=255)
                                      ccGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCt                                                                                                                <  1:761860/150‑1 (MQ=255)
                                        gAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGTATTAAACTGATGACAAACGCAAATCTGCCTGa                                                                                                              >  2:505996/1‑150 (MQ=255)
                                        gAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGa                                                                                                              >  1:922500/1‑150 (MQ=255)
                                        gAAACAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATAAGGGGATGATGTTAGGGTTTACGGTATAAAAAGAACGTAATCGGAAATATGCCCTGATCTCTCCTGGCATTAAACTGATGACAAAAGCAAAACTGCATGa                                                                                                              >  2:205687/1‑150 (MQ=255)
                                         aaaGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGAt                                                                                                             <  2:573953/150‑1 (MQ=255)
                                           aGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATAAGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATgc                                                                                                           <  1:125584/150‑1 (MQ=255)
                                             aaCATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATgcgc                                                                                                         >  2:890744/1‑150 (MQ=255)
                                               cATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTa                                                                                                       <  1:1421650/150‑1 (MQ=255)
                                                 tGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACg                                                                                                     >  1:1405492/1‑150 (MQ=255)
                                                     tataTGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTATGCTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTa                                                                                                 <  1:1111581/150‑1 (MQ=255)
                                                          ggCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAgg                                                                                            <  2:1321012/150‑1 (MQ=255)
                                                           gcgcACTGCCATCTGCTGGAGCATGAAGAGACGGGGATGAGGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTATCCGGGCATTAAACTGATGACAAACTCAAATCTGCCTGATGCGCTCCGCTTATCAGGc                                                                                           >  2:257153/1‑150 (MQ=255)
                                                            cgcTCTGCCATGTGCTGGAGCATGAAGATACGTTGATGTTGTTAGGGTTTACGGTATAATAAGAACGTATTCGGAAATATGCCCGGAGCTCTCCGGGCATTAAACTGATTACAAACGCAAATCTGCCTGATGCGCTACTCTTATCAGGcc                                                                                          <  1:105487/150‑1 (MQ=255)
                                                            cgcACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGcc                                                                                          <  1:638370/150‑1 (MQ=255)
                                                                  gCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCt                                                                                    <  2:1004484/150‑1 (MQ=255)
                                                                  gCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCt                                                                                    >  2:179467/1‑150 (MQ=255)
                                                                     aTCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGc                                                                                 <  1:634318/150‑1 (MQ=255)
                                                                          ctgGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGc                                                                            <  1:138605/150‑1 (MQ=255)
                                                                              aGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATc                                                                        <  2:1288246/150‑1 (MQ=255)
                                                                              aGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATc                                                                        <  2:1281537/150‑1 (MQ=255)
                                                                               gCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATcc                                                                       >  1:1080384/1‑150 (MQ=255)
                                                                                cATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCa                                                                      >  1:124166/1‑150 (MQ=255)
                                                                                cATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCa                                                                      <  2:503948/150‑1 (MQ=255)
                                                                                cATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCa                                                                      >  2:870568/1‑150 (MQ=255)
                                                                                cATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCa                                                                      >  2:21711/1‑150 (MQ=255)
                                                                                  tGAAGATACGGGGATGATGTTAGGGTTTACGGTATAACAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCTGGCCTACGCTAGCTGTGCAATCCAtt                                                                    <  1:1062394/150‑1 (MQ=255)
                                                                                  tGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCAtt                                                                    <  2:618534/150‑1 (MQ=255)
                                                                                  tGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCAtt                                                                    <  1:1149513/150‑1 (MQ=255)
                                                                                        tACCGGCATGATGCTAGGGTTTACGGAATAACAAGCCCGTATGCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCACATCTGCCTGATGCGCTACGCTTATCAGGCCTTCGCTAGCTGTGCAATCCATTGAtttt                                                              <  1:299603/150‑1 (MQ=38)
                                                                                         aCGGGGTTGATGTTAGGGTGTAGGTTATAAATAGACCGTATTCGGATTTATGACCGTATCTATCCGGGCATTAAACTGATGAGAATCGCAAATCTGCCTGTTGCGCTACGCTTATCAGGCCTACGCTTGCTGTGCATTCCATTGATTTTg                                                             <  1:433746/150‑1 (MQ=25)
                                                                                         aCGGGGATGATGTTAGGGTTTACTTTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCGTTGATTTTg                                                             <  1:675504/150‑1 (MQ=255)
                                                                                           ggggATGATGTTATGGTTTACGTTATAAAAATACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAATCTGATGACTAACGCTAATCTGCCTGTTGCGCTTCGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTggt                                                           <  1:205687/150‑3 (MQ=255)
                                                                                             ggATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCGTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGcaca                                                         >  2:564131/1‑150 (MQ=255)
                                                                                             ggATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACTCTAGCTGTGCAATCCATTGATTTTGcaca                                                         <  1:1261368/150‑1 (MQ=255)
                                                                                             ggATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGcaca                                                         >  2:213248/1‑150 (MQ=255)
                                                                                                        gggTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCATTCCATTCATTTTGCACAATTTTGTAGGc                                              >  1:840470/1‑150 (MQ=255)
                                                                                                           tttACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTgg                                           <  2:1536591/150‑1 (MQ=255)
                                                                                                            ttACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCAGTGATTTTGCACAATTTTGTAGGCTTGa                                          >  2:362836/1‑150 (MQ=255)
                                                                                                             tACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGAATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGAt                                         <  2:1282492/150‑1 (MQ=255)
                                                                                                                 gTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAgg                                     >  1:1381808/1‑150 (MQ=255)
                                                                                                                 gTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAgg                                     >  1:504019/1‑150 (MQ=255)
                                                                                                                      aaaaGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTc                                >  2:1358366/1‑150 (MQ=255)
                                                                                                                        aaGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCAc                              >  2:1498213/1‑150 (MQ=255)
                                                                                                                         aGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCAcg                             >  2:711955/1‑150 (MQ=255)
                                                                                                                         aGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCAcg                             >  2:938475/1‑150 (MQ=255)
                                                                                                                         aGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCAcg                             >  2:1069911/1‑150 (MQ=255)
                                                                                                                         aGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCAcg                             <  2:1080384/150‑1 (MQ=255)
                                                                                                                           aCCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCAcgcc                           >  2:1580795/1‑150 (MQ=255)
                                                                                                                             cGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCAcgccgc                         >  2:527869/1‑150 (MQ=255)
                                                                                                                             cGTATTCGGAAATATGCCCGGATCTCCCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCAcgccgc                         >  2:40048/1‑150 (MQ=255)
                                                                                                                               tATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCAt                       <  1:436600/150‑1 (MQ=255)
                                                                                                                                   cGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCgg                   >  2:1601985/1‑150 (MQ=255)
                                                                                                                                    ggAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGc                  >  2:355332/1‑150 (MQ=255)
                                                                                                                                    ggAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGc                  <  2:1537082/150‑1 (MQ=255)
                                                                                                                                     gAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCa                 <  2:124166/150‑1 (MQ=255)
                                                                                                                                        taatGCCCGGAGAGATCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGa              >  1:203314/3‑150 (MQ=255)
                                                                                                                                         aatGCCCGGAGAGATCCGGGCATTAAACTGATGACAAACGCATATCTGCCTGATGCGCTACGCTTATCTGGCCTACGGTAGCTGTGCAATCCTTTGTTTTTGCAGAATTTTGTAGGCTGGATTAGGCGTTCACGCCGCGGACGGCATGaa             >  1:1524850/2‑150 (MQ=21)
                                                                                                                                          atGCCCGGAGAGATCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaac            >  2:511990/1‑150 (MQ=255)
                                                                                                                                          atGCCCGGAGAGATCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaac            >  2:893173/1‑150 (MQ=255)
                                                                                                                                          atGCCCGGAGAGATCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaac            >  2:1618361/1‑150 (MQ=255)
                                                                                                                                          atGCCCGGAGAGATCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaac            >  1:1393440/1‑150 (MQ=255)
                                                                                                                                           tGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaaca           >  1:559909/1‑150 (MQ=255)
                                                                                                                                           tGCCCGGAGCTCTCCGGGCATTAAACTGATGACACACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTCGGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaaca           <  1:201999/150‑1 (MQ=255)
                                                                                                                                            gCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaacaa          >  2:48162/1‑150 (MQ=255)
                                                                                                                                            gCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaacaa          <  1:27543/150‑1 (MQ=255)
                                                                                                                                             cccGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaacaac         >  1:1400442/1‑150 (MQ=255)
                                                                                                                                              ccGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcg        <  2:1405492/150‑1 (MQ=255)
                                                                                                                                              ccGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcg        <  1:1015590/150‑1 (MQ=255)
                                                                                                                                              ccGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAcg        >  1:729876/1‑150 (MQ=255)
                                                                                                                                                   tctctcCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCAc   >  2:268902/1‑150 (MQ=255)
                                                                                                                                                    ctctcCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACg  <  1:176453/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
ATGTCAGCGAAGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACG  >  NC_000913/138515‑138812

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