Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,885,807 T→A 5.2% intergenic (+18/‑38) yeaW → / → yeaX putative YeaWX dioxygenase alpha subunit; 2Fe‑2S cluster/putative YeaWX dioxygenase beta subunit, reductase component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,885,8070TA 5.2% 275.4 / 2.7 95intergenic (+18/‑38)yeaW/yeaXputative YeaWX dioxygenase alpha subunit; 2Fe‑2S cluster/putative YeaWX dioxygenase beta subunit, reductase component
Reads supporting (aligned to +/- strand):  ref base T (43/47);  new base A (3/2);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

GTACAGAAAGGGCTGAAATCGCGTGGCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGCAAACCATTACCTGCG  >  NC_000913/1885658‑1885952
                                                                                                                                                     |                                                                                                                                                 
gTACAGAATGGCCTGAAATCGCGTGGCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTAt                                                                                                                                                   <  2:1095082/150‑1 (MQ=255)
  aaaGAAAGCGCTGAAATCGCGTGGCTATCGTGGTCAGGGGCGCATCATGGACGACAGTAGCGGTAGTGGCATTTCCGTACATGGTATCGCCCAATTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTAttt                                                                                                                                                 >  2:861837/3‑150 (MQ=255)
     caaaGGGCTGAAATCGCGTGGCTATCGTGGGCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTCTCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTcc                                                                                                                                              <  1:834887/149‑1 (MQ=255)
         gggCTGAAATCGCGTGGCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTccgccg                                                                                                                                          >  2:1442190/1‑150 (MQ=255)
           gCTGAAATCGCGTGGCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGAAAATACcgccgcc                                                                                                                                        <  2:1152667/150‑2 (MQ=255)
             tGAAATCGCGTGGCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTgg                                                                                                                                      >  1:283161/1‑150 (MQ=255)
             tGAAATCGCGTGGCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGAAAATACcgccgccga                                                                                                                                      <  1:1495062/150‑4 (MQ=255)
                 aTCGCGTGGCTATCGTGGTCAGGGGCGCATCATGGCCAACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTg                                                                                                                                  >  1:319339/1‑150 (MQ=255)
                   cgcgTGGCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGAt                                                                                                                                <  1:958301/150‑1 (MQ=255)
                   cgcgTGGCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGAtt                                                                                                                               >  1:216534/1‑150 (MQ=255)
                     cgTGGCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTTATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGAttt                                                                                                                              >  1:808098/1‑150 (MQ=255)
                      gTGGCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGAtttt                                                                                                                             >  2:509455/1‑150 (MQ=255)
                        ggCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTTGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTg                                                                                                                           <  1:1584660/150‑1 (MQ=255)
                        ggCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTCGTGGCATTTCCGAACATGGTATCGCCCATTTCCATCATCTGCGGGCGCAGGTGGTTATGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTg                                                                                                                           <  1:193649/150‑1 (MQ=255)
                        ggCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTg                                                                                                                           >  2:1208665/1‑150 (MQ=255)
                          cTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGAt                                                                                                                         <  1:1442190/150‑1 (MQ=255)
                          cTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGGAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGAt                                                                                                                         <  1:167024/150‑1 (MQ=255)
                            aTCGTGGTCAGGGGCGCATCATGGCCTACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATgg                                                                                                                       <  2:319339/150‑1 (MQ=255)
                                 ggTCAGGGGCGCATCATGGCGGACAGTAGCGGTTTGGGCATTTCCGAACATGGTATCGCCCATTTCCATACTCTGCTGGCGCAGGTGTTTAAGGACTAATGTCATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTAc                                                                                                                  <  2:669267/150‑1 (MQ=255)
                                  gTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGAGGGAGTACa                                                                                                                 >  2:138434/1‑150 (MQ=255)
                                    cAGGGGCGCATCATGGCTGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGc                                                                                                               >  2:1585177/1‑150 (MQ=255)
                                     aGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCa                                                                                                              >  1:1402610/1‑150 (MQ=255)
                                      ggggCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCaa                                                                                                             >  1:200542/1‑150 (MQ=255)
                                           gcATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCa                                                                                                        <  1:776834/150‑1 (MQ=255)
                                                atGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTa                                                                                                   >  2:1346/1‑150 (MQ=255)
                                                  ggCTGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATc                                                                                                 <  1:1585177/150‑1 (MQ=255)
                                                    ccGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCaa                                                                                               >  2:888478/1‑150 (MQ=255)
                                                      gACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAAt                                                                                             >  2:212226/1‑150 (MQ=255)
                                                       aCAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATg                                                                                            >  1:76519/1‑150 (MQ=255)
                                                         aGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGTCTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGtt                                                                                          <  2:761156/150‑1 (MQ=255)
                                                         aGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGGTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGtt                                                                                          >  1:1056009/1‑150 (MQ=255)
                                                         aGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGtt                                                                                          <  1:474907/150‑1 (MQ=255)
                                                            tgCGGTAGTGGCATTTGCGAACATCGTATCGCCCATTTCCATATTCTGCGCGCGCAGGTGTTTTAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGAAGGAGTACAGCAATGTCAGACAATCAAATGTTTgt                                                                                       <  2:918510/149‑2 (MQ=255)
                                                            aGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGa                                                                                       <  1:801850/150‑1 (MQ=255)
                                                               ggTAGTGGCATTTCCGAACATGGTATAGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGt                                                                                    >  1:810291/1‑150 (MQ=255)
                                                                   gTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACa                                                                                 <  2:1134042/149‑1 (MQ=255)
                                                                   gTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACa                                                                                 >  1:1134042/1‑149 (MQ=255)
                                                                   gTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAg                                                                                >  2:741485/1‑150 (MQ=255)
                                                                    tGGTATTTCCGAACATGGTATCGCCCATTTCCATAATGTGCTGGGGCAGGTGTTTAAGGACTATTGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAgg                                                                               <  1:1149132/150‑1 (MQ=255)
                                                                     ggCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGt                                                                              <  1:398216/150‑1 (MQ=255)
                                                                      gCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTg                                                                             >  1:423932/1‑150 (MQ=255)
                                                                        aTTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAg                                                                           <  1:1208665/150‑1 (MQ=255)
                                                                         tttCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGc                                                                          >  1:1632567/1‑150 (MQ=255)
                                                                           tCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTTATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCa                                                                        <  2:808098/150‑1 (MQ=255)
                                                                             cGAACATGGGATCGCCAATTTCCATAATCTGCTTGCGCAGGTGTTAAAGGACTAATGAAAGCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGACTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAgg                                                                      <  2:869574/150‑1 (MQ=255)
                                                                                aCATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTg                                                                   <  1:509455/150‑1 (MQ=255)
                                                                                    ggTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAAcc                                                               <  2:1698697/150‑1 (MQ=255)
                                                                                    ggTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAAcc                                                               <  2:567848/150‑1 (MQ=255)
                                                                                      tATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAAcccc                                                             >  1:577554/1‑150 (MQ=255)
                                                                                       aTCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCt                                                            <  1:481468/150‑1 (MQ=255)
                                                                                       aTCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCt                                                            >  1:1667850/1‑150 (MQ=255)
                                                                                       aTCGCCCAATTCCATTATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCt                                                            <  2:447606/150‑1 (MQ=255)
                                                                                          gCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTAc                                                         <  2:1402610/150‑1 (MQ=255)
                                                                                            ccATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCg                                                       <  2:283161/150‑1 (MQ=255)
                                                                                            ccATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCg                                                       <  2:76519/150‑1 (MQ=255)
                                                                                                ttCCATAATCTGCTGGCGCATGTGATTAAGGACTAATGACATCGGCGGCGGTATTTTCGGCCGCGGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAATTACATGTGAGCCAGGTTGTACCCCTTACCGAACa                                                   <  2:1187433/150‑1 (MQ=255)
                                                                                                ttCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACa                                                   >  1:1396640/1‑150 (MQ=255)
                                                                                                ttCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACa                                                   >  2:1247832/1‑150 (MQ=255)
                                                                                                  ccATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTAGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAgg                                                 <  1:981859/150‑1 (MQ=255)
                                                                                                   cATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGt                                                >  2:1439220/1‑150 (MQ=255)
                                                                                                   cATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGt                                                >  1:224910/1‑150 (MQ=255)
                                                                                                    aTAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTg                                               <  2:238045/150‑1 (MQ=255)
                                                                                                    aTAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTg                                               <  1:795424/150‑1 (MQ=255)
                                                                                                      aaTCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGaa                                             <  1:277335/150‑1 (MQ=255)
                                                                                                       aTCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAAc                                           >  2:1343240/1‑150 (MQ=255)
                                                                                                           gctgGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCt                                        <  2:1632567/150‑1 (MQ=255)
                                                                                                              ggCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCa                                     >  2:1356698/1‑150 (MQ=255)
                                                                                                                cgcAGGTGTTTATGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACg                                   >  2:991738/1‑150 (MQ=255)
                                                                                                                cgcAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACg                                   >  2:952247/1‑150 (MQ=255)
                                                                                                                  cAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCt                                 >  2:1504696/1‑150 (MQ=255)
                                                                                                                   aGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCtg                                <  1:245461/150‑1 (MQ=255)
                                                                                                                    ggTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCAAACAGGTGAAACGCTTCACGCtgg                               >  2:29041/1‑150 (MQ=255)
                                                                                                                    ggTGTGTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCGGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCtgg                               <  2:1700222/150‑1 (MQ=255)
                                                                                                                     gtgtTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCtggt                              <  1:1163501/150‑1 (MQ=255)
                                                                                                                     gtgtTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCAAACAGGTGAAACGCTTCACGCtggt                              >  2:574799/1‑150 (MQ=255)
                                                                                                                        tttAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGc                           >  1:947675/1‑150 (MQ=255)
                                                                                                                        tttAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGc                           >  1:941523/1‑150 (MQ=255)
                                                                                                                           aaGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGAATTTTGATTGAGTACAGCAATGTCAGACTATCAACAGTTTCAAGTACAGGTGACCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAAc                        <  1:861837/150‑1 (MQ=255)
                                                                                                                              gccTAATGACATCGGCGGCGGTATTTTCCGCCGCTGTGCTGATTTTTTATGGAGTACAGCAATTTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGa                     <  1:465139/148‑1 (MQ=255)
                                                                                                                              gACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGa                     <  1:458755/150‑1 (MQ=255)
                                                                                                                               aCTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGAt                    >  1:765129/1‑150 (MQ=255)
                                                                                                                                 tAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGGGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATgg                  <  1:1358435/150‑1 (MQ=255)
                                                                                                                                 tAATGACATCGGCGGCGGTATTTTACGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATgg                  >  1:1682863/1‑150 (MQ=255)
                                                                                                                                  aaTGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGc                 <  1:587748/150‑1 (MQ=255)
                                                                                                                                  aaTGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGc                 >  2:1249082/1‑150 (MQ=255)
                                                                                                                                   aTGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGCa                <  2:1396640/150‑1 (MQ=255)
                                                                                                                                   aTGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGCa                <  1:212226/150‑1 (MQ=255)
                                                                                                                                     gACATCGGCGGCGCTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGCaaa              <  1:589173/150‑1 (MQ=255)
                                                                                                                                      ccATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATTTTTGAAGTACAGCTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGCAAAc             <  1:593954/149‑1 (MQ=255)
                                                                                                                                        aTCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGCAAACCa           >  1:197598/1‑150 (MQ=255)
                                                                                                                                         ccagcggcggAAAATACCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGCCTATCAAATGTTTGAAGTACCGGGGAGCCAGGTTGCACCCCTTACCGAACAGGGGAAACGCTTCACGCTGGTGGCACCCGAAGGCAAACCAt          >  1:1061811/4‑150 (MQ=255)
                                                                                                                                         ccagcggcggAAAATACCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGATCCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGCAAACCAt          >  1:661904/4‑150 (MQ=255)
                                                                                                                                         ccagcggcggAAAATACCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGCAAACCaa          >  1:604590/4‑149 (MQ=255)
                                                                                                                                            gcggcggTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGCAAACCATTAc       >  2:378356/1‑150 (MQ=255)
                                                                                                                                                cggTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGACCAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGCAAACCATTACCTGc   <  1:307569/150‑1 (MQ=255)
                                                                                                                                                cggTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACATGTGCAACGCTTCACGCTGGTGGCAACCGATGGCAAACCATTACCTGc   >  1:890581/1‑150 (MQ=255)
                                                                                                                                                 ggTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGCAAACCATTACCTGCg  >  2:3825/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
GTACAGAAAGGGCTGAAATCGCGTGGCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGCAAACCATTACCTGCG  >  NC_000913/1885658‑1885952

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