Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,252,812 A→T 6.9% intergenic (+24/+83) yeiI → / ← nupX putative kinase/nucleoside permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,252,8120AT 6.9% 243.5 / 7.3 72intergenic (+24/+83)yeiI/nupXputative kinase/nucleoside permease
Reads supporting (aligned to +/- strand):  ref base A (35/32);  new base T (3/2);  total (38/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.57e-01

ATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTCGCCGCC  >  NC_000913/2252665‑2252958
                                                                                                                                                   |                                                                                                                                                  
aTACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTaa                                                                                                                                                   <  4:744221/149‑1 (MQ=255)
 tACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAAc                                                                                                                                                  >  1:332124/1‑149 (MQ=255)
    aGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAAcgcg                                                                                                                                               >  3:644329/1‑149 (MQ=255)
         ccGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCGGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTa                                                                                                                                          >  2:275330/1‑149 (MQ=255)
          cgcgACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAAACTTGCCGGATGCGACGCTTAACGCGTCTTAt                                                                                                                                         >  4:52566/1‑149 (MQ=255)
           gcgACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATAAGACGCGTTAAGCGTcgcatc                                                                                                                                        <  1:584485/149‑6 (MQ=255)
             gACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCg                                                                                                                                      >  1:757608/1‑149 (MQ=255)
               cAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGc                                                                                                                                    <  1:481438/149‑1 (MQ=255)
                   gcCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTAc                                                                                                                                >  2:195589/1‑149 (MQ=255)
                       gTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAAt                                                                                                                            >  4:182251/1‑149 (MQ=255)
                         tttGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCa                                                                                                                          <  3:425633/149‑1 (MQ=255)
                           tGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAgg                                                                                                                        <  3:221242/149‑1 (MQ=255)
                             cGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATAAGACGCGTTAAGCGTCGCATCCGGcaac                                                                                                                                 <  1:431948/138‑4 (MQ=38)
                             cGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATAAGACGCGTTAAGCGTCGCATCCGGcaac                                                                                                                                 >  2:431948/1‑135 (MQ=38)
                              ggTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGTCCTACAAATCAGGCaa                                                                                                                     >  3:89500/1‑149 (MQ=255)
                                tAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAAtt                                                                                                                   <  4:743641/149‑1 (MQ=255)
                                           gCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTaa                                                                                                        >  2:207769/1‑149 (MQ=255)
                                             aaTTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAAt                                                                                                      <  3:602736/149‑1 (MQ=255)
                                             aaTTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAAt                                                                                                      <  4:248334/149‑1 (MQ=255)
                                               tttCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGt                                                                                                    <  1:167461/149‑1 (MQ=255)
                                                 tCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAg                                                                                                  <  1:265632/149‑1 (MQ=255)
                                                   gcgcgcCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGc                                                                                                <  2:45708/149‑1 (MQ=255)
                                                      cgcCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTg                                                                                             <  3:320565/149‑1 (MQ=255)
                                                       gcCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGAGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAggtatga                                                                                            >  2:517161/1‑144 (MQ=255)
                                                          aGCGGCAGCTCAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCACGCAATTTCAGCAAGTAAATGTAGGCCTGATaa                                                                                         <  4:483660/149‑1 (MQ=255)
                                                            cGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGc                                                                                       <  1:298766/149‑1 (MQ=255)
                                                             ggCAGCTTAAACAACCCTACCCAGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCt                                                                                      <  3:464968/149‑1 (MQ=255)
                                                               cAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTg                                                                                    <  1:302866/149‑1 (MQ=255)
                                                                aGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTgc                                                                                   <  4:207153/149‑1 (MQ=255)
                                                                aGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTgc                                                                                   <  1:689896/149‑1 (MQ=255)
                                                                 gCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTgcg                                                                                  >  1:261227/1‑149 (MQ=255)
                                                                   ttAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGc                                                                                                <  2:268143/133‑1 (MQ=255)
                                                                   ttAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGc                                                                                                >  1:268143/1‑133 (MQ=255)
                                                                   ttAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCa                                                                                >  2:682326/1‑149 (MQ=255)
                                                                     aaaCAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATc                                                                              <  1:195589/149‑1 (MQ=255)
                                                                            cctaccctGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAACGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAAt                                                                       >  2:661439/1‑149 (MQ=255)
                                                                                 cctGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAAc                                                                                                                                                  >  1:512834/1‑69 (MQ=255)
                                                                                 cctGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAAc                                                                                                                                                  <  2:512834/69‑1 (MQ=255)
                                                                                         ccGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTa                                                          >  1:121953/1‑149 (MQ=255)
                                                                                            aTAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAAc                                                       <  1:661439/149‑1 (MQ=255)
                                                                                               aCGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCaa                                                    >  4:229030/1‑149 (MQ=255)
                                                                                                 gcgcTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATa                                                  >  2:295563/1‑149 (MQ=255)
                                                                                                   gcTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATaaa                                                <  3:182251/149‑1 (MQ=255)
                                                                                                    cTTTCATTAGTGCCAATGGTGTAACAATGATGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAg                                               <  2:692608/149‑1 (MQ=255)
                                                                                                     tttCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGa                                              >  2:565951/1‑149 (MQ=255)
                                                                                                     tttCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGa                                              >  4:265905/1‑149 (MQ=255)
                                                                                                      ttCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGaa                                             <  2:164231/149‑1 (MQ=255)
                                                                                                         aTTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAAccc                                          >  4:696601/1‑149 (MQ=255)
                                                                                                          ttAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAAc                                                       <  4:755936/135‑1 (MQ=255)
                                                                                                          ttAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAAc                                                       >  3:755936/1‑135 (MQ=255)
                                                                                                             gTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCa                                      >  4:389547/1‑149 (MQ=255)
                                                                                                                  aaTGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCaa                                                              <  4:64176/120‑1 (MQ=255)
                                                                                                                  aaTGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCaa                                                              >  3:64176/1‑120 (MQ=255)
                                                                                                                  aaTGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGt                                 <  1:682326/149‑1 (MQ=255)
                                                                                                                      gtgtAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCa                             >  1:645278/1‑149 (MQ=255)
                                                                                                                       tgtAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCAc                            >  1:726466/1‑149 (MQ=255)
                                                                                                                         tAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACTAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACtc                          >  3:627435/1‑149 (MQ=255)
                                                                                                                         tAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACtc                          >  4:128456/1‑149 (MQ=255)
                                                                                                                          aaCAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACtca                         >  3:413723/1‑149 (MQ=255)
                                                                                                                             aaTGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACtcatca                      >  3:137688/1‑149 (MQ=255)
                                                                                                                                gTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGt                   >  1:91935/1‑149 (MQ=255)
                                                                                                                                 ttGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGtt                  <  1:395455/149‑1 (MQ=255)
                                                                                                                                 ttGCCCGATGCGACGCTTAACGCGTCTTATACGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGtt                  <  4:392576/149‑1 (MQ=255)
                                                                                                                                 aggCCGGATAAGACGCGTTAAGCGTCGCATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGtt                  >  3:634529/3‑149 (MQ=38)
                                                                                                                                   gCCGGATAAGACGCGTTAAGCGTCGCATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTgg                >  2:77812/1‑149 (MQ=255)
                                                                                                                                    ccGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGa               >  2:415328/1‑149 (MQ=255)
                                                                                                                                      ggATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGaaa             <  1:207769/149‑1 (MQ=255)
                                                                                                                                           cGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTc        <  3:128456/149‑1 (MQ=255)
                                                                                                                                            gACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTcg       >  2:117949/1‑149 (MQ=255)
                                                                                                                                             aCGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTcgc      >  4:154351/1‑149 (MQ=255)
                                                                                                                                               gCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAAccc                                          >  1:23219/1‑111 (MQ=255)
                                                                                                                                               gCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAAccc                                          <  2:23219/111‑1 (MQ=255)
                                                                                                                                                 ttAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTcgccgcc  <  2:261227/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATACAGTTTCCGCGACAGCGCCCGTTTTGCGGTAGCCTGCGCGGCAATTTCGCGCGCCAGCGGCAGCTTAAACAACCCTACCCTGTCTGCCGATAACGCGCTTTCATTAGTGCCAATGGTGTAACAATGTTGCCGGATGCGACGCTTAACGCGTCTTATCCGGCCTACAAATCAGGCAATTTCAGCAAGTAAATGTAGGCCTGATAAGCTTGCGCATCAGGCAATGACACTCAAGCTAAACCAATAAAGAACCCGGCAATGGTCGCACTCATCAAGTTGGAAAGCGTCGCCGCC  >  NC_000913/2252665‑2252958

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: