Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,885,805 T→A 5.0% intergenic (+16/‑40) 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,8050TA 5.0% 278.8 / 3.9 100intergenic (+16/‑40)yeaW/yeaXputative YeaWX dioxygenase alpha subunit; 2Fe‑2S cluster/putative YeaWX dioxygenase beta subunit, reductase component
Reads supporting (aligned to +/- strand):  ref base T (46/49);  new base A (3/2);  total (49/51)
Fisher's exact test for biased strand distribution p-value = 6.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.40e-01

CGTACAGAAAGGGCTGAAATCGCGTGGCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGCAAACCATTACCTGCG  >  NC_000913/1885657‑1885952
                                                                                                                                                    |                                                                                                                                                   
cGTACAGAAAGGGCTGAAATCGCGTGGCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTa                                                                                                                                                    >  2:474907/1‑150 (MQ=255)
cGTACAGAAAGGGCTGAAATCGCGTGGCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTa                                                                                                                                                    >  2:245461/1‑150 (MQ=255)
cGTACAGAAAGGGCTGAAATCGCGTGGCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTa                                                                                                                                                    <  1:518942/150‑1 (MQ=255)
 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:1667850/1‑150 (MQ=255)
                                                                                        aTCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCt                                                            <  1:481468/150‑1 (MQ=255)
                                                                                        aTCGCCCAATTCCATTATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCt                                                            <  2:447606/150‑1 (MQ=255)
                                                                                           gCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTAc                                                         <  2:1402610/150‑1 (MQ=255)
                                                                                             ccATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCg                                                       <  2:76519/150‑1 (MQ=255)
                                                                                             ccATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCg                                                       <  2:283161/150‑1 (MQ=255)
                                                                                                 ttCCATAATCTGCTGGCGCATGTGATTAAGGACTAATGACATCGGCGGCGGTATTTTCGGCCGCGGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAATTACATGTGAGCCAGGTTGTACCCCTTACCGAACa                                                   <  2:1187433/150‑1 (MQ=255)
                                                                                                 ttCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACa                                                   >  2:1247832/1‑150 (MQ=255)
                                                                                                 ttCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACa                                                   >  1:1396640/1‑150 (MQ=255)
                                                                                                   ccATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTAGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAgg                                                 <  1:981859/150‑1 (MQ=255)
                                                                                                    cATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGt                                                >  1:224910/1‑150 (MQ=255)
                                                                                                    cATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGt                                                >  2:1439220/1‑150 (MQ=255)
                                                                                                     aTAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTg                                               <  1:795424/150‑1 (MQ=255)
                                                                                                     aTAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTg                                               <  2:238045/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                 >  2:1249082/1‑150 (MQ=255)
                                                                                                                                   aaTGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGc                 <  1:587748/150‑1 (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)
                                                                                                                                                    |                                                                                                                                                   
CGTACAGAAAGGGCTGAAATCGCGTGGCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGCAAACCATTACCTGCG  >  NC_000913/1885657‑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: