Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 138,665 C→A 9.4% intergenic (+32/+170) 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,6650CA 9.4% 321.3 / 12.8 107intergenic (+32/+170)cueO/gcdmulticopper oxidase (laccase)/glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (46/50);  new base A (7/4);  total (53/54)
Fisher's exact test for biased strand distribution p-value = 3.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 117 total aligned reads displayed.

GTCAGCGAAGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTT  >  NC_000913/138517‑138814
                                                                                                                                                    |                                                                                                                                                     
gTCAGCGAAGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGAtctct                                                                                                                                                      >  2:1593677/1‑150 (MQ=255)
gTCAGCGAAGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGAtctct                                                                                                                                                      >  1:1084293/1‑150 (MQ=255)
gTCAGCGAAGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGAtctct                                                                                                                                                      >  1:94795/1‑150 (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                                                                                                                                         >  2:736497/1‑150 (MQ=255)
             tGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAAc                                                                                                                                         >  1:1536591/1‑150 (MQ=255)
             tGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAAc                                                                                                                                         <  1:380731/150‑1 (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                                                                                                                                  >  2:176453/1‑150 (MQ=255)
                    gTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGAc                                                                                                                                  <  1:971446/150‑1 (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)
                                                               tGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGc                                                                                       >  2:793284/1‑150 (MQ=255)
                                                                gCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCt                                                                                      >  2:179467/1‑150 (MQ=255)
                                                                gCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCt                                                                                      <  2:1004484/150‑1 (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                                                                        >  2:21711/1‑150 (MQ=255)
                                                                              cATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCa                                                                        >  2:870568/1‑150 (MQ=255)
                                                                              cATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCa                                                                        >  1:124166/1‑150 (MQ=255)
                                                                              cATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCa                                                                        <  2:503948/150‑1 (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:53141/150‑1 (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              >  1:1393440/1‑150 (MQ=255)
                                                                                                                                        atGCCCGGAGAGATCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaac              >  2:511990/1‑150 (MQ=255)
                                                                                                                                        atGCCCGGAGAGATCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaac              >  2:1618361/1‑150 (MQ=255)
                                                                                                                                        atGCCCGGAGAGATCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGaac              >  2:893173/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)
                                                                                                                                                    ctcCGGGCATTAAACTGATGACAAACGCAAATCTGCCGGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGtt  >  2:903720/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GTCAGCGAAGTGCTGGTGAAGTTTAATCACGATGCACCGAAAGAACATGCTTATATGGCGCACTGCCATCTGCTGGAGCATGAAGATACGGGGATGATGTTAGGGTTTACGGTATAAAAAGACCGTATTCGGAAATATGCCCGGATCTCTCCGGGCATTAAACTGATGACAAACGCAAATCTGCCTGATGCGCTACGCTTATCAGGCCTACGCTAGCTGTGCAATCCATTGATTTTGCACAATTTTGTAGGCTGGATAAGGCGTTCACGCCGCATCCGGCATGAACAACGCGCACGTT  >  NC_000913/138517‑138814

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: