Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,885,808 T→A 5.2% intergenic (+19/‑37) 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,8080TA 5.2% 287.8 / 4.0 95intergenic (+19/‑37)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 = 7.69e-01

ACAGAAAGGGCTGAAATCGCGTGGCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGCAAACCATTACCTGCG  >  NC_000913/1885660‑1885952
                                                                                                                                                    |                                                                                                                                                
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                                                                                 >  1:1134042/1‑149 (MQ=255)
                                                                 gTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACa                                                                                 <  2:1134042/149‑1 (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:567848/150‑1 (MQ=255)
                                                                                  ggTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAAcc                                                               <  2:1698697/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: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                 <  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)
                                                                                                                                                    |                                                                                                                                                
ACAGAAAGGGCTGAAATCGCGTGGCTATCGTGGTCAGGGGCGCATCATGGCCGACAGTAGCGGTAGTGGCATTTCCGAACATGGTATCGCCCATTTCCATAATCTGCTGGCGCAGGTGTTTAAGGACTAATGACATCGGCGGCGGTATTTTCCGCCGCTGGGCTGATTTTTGATGGAGTACAGCAATGTCAGACTATCAAATGTTTGAAGTACAGGTGAGCCAGGTTGAACCCCTTACCGAACAGGTGAAACGCTTCACGCTGGTGGCAACCGATGGCAAACCATTACCTGCG  >  NC_000913/1885660‑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: