Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,076,391 C→A 5.1% T28T (ACC→ACA yeeS → CP4‑44 prophage; putative DNA repair protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,076,3910CA 5.1% 273.7 / 2.7 79T28T (ACC→ACAyeeSCP4‑44 prophage; putative DNA repair protein
Reads supporting (aligned to +/- strand):  ref base C (37/38);  new base A (2/2);  total (39/40)
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.53e-01

TCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCA  >  NC_000913/2076245‑2076537
                                                                                                                                                  |                                                                                                                                                  
tCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTc                                                                                                                                                  <  7:390406/149‑1 (MQ=255)
   gTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCAc                                                                                                                                               <  3:192352/149‑1 (MQ=255)
      tGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCAccc                                                                                                                                             <  8:68343/148‑1 (MQ=255)
          ggCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTgca                                                                                                                                        <  8:226199/149‑1 (MQ=255)
               tGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACg                                                                                                                                   >  7:193199/1‑149 (MQ=255)
               tGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACg                                                                                                                                   >  2:80226/1‑149 (MQ=255)
                 ttGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACgg                                                                                                                                  <  1:80226/148‑1 (MQ=255)
                 ttGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACg                                                                                                                                   <  6:408308/147‑1 (MQ=255)
                   gAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGaa                                                                                                                               >  3:365902/1‑149 (MQ=255)
                   gAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGaa                                                                                                                               >  1:305900/1‑149 (MQ=255)
                    aacccctGAGGAATAAGATCCGCGAACGTAAAGTAGAGCAGAGAAGACGCCCGGCGAGCGCCGCCCCATTCCACCGGCCCCGCAAACCCAGGACCGCCATCATCCTGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGaa                                                                                                                               <  6:74562/142‑1 (MQ=255)
                     aaCCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACATCCACCCGTgca                                                                                                                                        <  2:418552/138‑1 (MQ=255)
                     aaCCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACATCCACCCGTgca                                                                                                                                        >  1:412676/1‑138 (MQ=255)
                     aaCCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACATCCACCCGTgca                                                                                                                                        >  1:418552/1‑138 (MQ=255)
                     aaCCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACATCCACCCGTgca                                                                                                                                        <  2:412676/138‑1 (MQ=255)
                       ccATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGc                                                                                                                           <  4:365902/149‑1 (MQ=255)
                          tGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGCGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGa                                                                                                                        >  4:10748/1‑149 (MQ=255)
                          tGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGa                                                                                                                        <  4:301885/149‑1 (MQ=255)
                            aGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGAtt                                                                                                                      <  4:290788/149‑1 (MQ=255)
                               aaTAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTg                                                                                                                   <  7:282428/149‑1 (MQ=255)
                                aTAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGCCCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGa                                                                                                                  <  4:207951/149‑1 (MQ=255)
                                 tAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGa                                                                                                                  <  5:91892/148‑1 (MQ=255)
                                 tAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGa                                                                                                                  >  6:91892/1‑148 (MQ=255)
                                 tAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCCCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGCACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATGCTGaa                                                                                                                 >  5:184038/1‑149 (MQ=255)
                                   aGATCAGCGAAGGTAAAGTAGAGTAGAGTTGACGCCCTGCGAGCGCAGCCTCTTTCAACTGGGCCTGAACACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCGACCCGTGCAGCACGGGAATGGCTGATTCTGAACa                                                                                                               <  4:333869/149‑1 (MQ=255)
                                    gATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACAt                                                                                                              >  1:52597/1‑149 (MQ=255)
                                    gATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATGCTGAACAt                                                                                                              >  1:200393/1‑149 (MQ=255)
                                     aTCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATg                                                                                                             >  7:188110/1‑149 (MQ=255)
                                         gcgAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAACCCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCgg                                                                                                         >  1:178008/1‑149 (MQ=255)
                                         gcgAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCgg                                                                                                         <  8:193199/149‑1 (MQ=255)
                                          cgAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCgg                                                                                                         >  3:127505/1‑148 (MQ=255)
                                               gTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGATCATGGCGGGACTg                                                                                                    >  1:57242/1‑148 (MQ=255)
                                               gTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTgg                                                                                                   >  5:271877/1‑149 (MQ=255)
                                                tAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGa                                                                                                  <  8:283555/149‑1 (MQ=255)
                                                tAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGa                                                                                                  <  3:386075/149‑1 (MQ=255)
                                                 aaaGTAGAGGAGAGATGACTCCCGGCGATCGCAGCCTCATTCAACTGGCACTGAAAACCCTGGACAGACATCTTAATGAACCCGGCGTGGCCTTCACCTCCTCTCGTGCAGCACGGGAATGGCTGATTCTTAACATGGCGGTAATGGa                                                                                                  <  3:223103/148‑1 (MQ=255)
                                                  aaGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGc                                                                                                >  1:231363/1‑149 (MQ=255)
                                                       gaggagAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTgaa                                                                                           <  2:178008/149‑1 (MQ=255)
                                                         ggagAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTg                                                                                             >  1:216168/1‑145 (MQ=255)
                                                         ggagAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTg                                                                                             <  2:216168/145‑1 (MQ=255)
                                                          gagagaTGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTgaagaa                                                                                        >  6:114064/1‑149 (MQ=255)
                                                            gagaTGCCGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAAtt                                                                                      <  3:10748/149‑1 (MQ=255)
                                                                tGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCgg                                                                                  <  5:114064/149‑1 (MQ=255)
                                                                 gACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGa                                                                                                       >  7:34127/1‑127 (MQ=255)
                                                                 gACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGa                                                                                                       <  8:34127/127‑1 (MQ=255)
                                                                    gCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGc                                                                              >  8:187541/1‑149 (MQ=255)
                                                                        ggCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTCCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTa                                                                          >  5:116117/1‑149 (MQ=255)
                                                                             gcgcAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGa                                                                     <  6:271877/149‑1 (MQ=255)
                                                                                 aGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGaacaa                                                                 >  5:304526/1‑149 (MQ=255)
                                                                                   ccTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATCAACCCGGCGTGGCCTTCACCTCCACCCCTGCAGCACGGGAATGGCTGAATCTGAACAAGGCGGGACTGGAGCGTGAAGAATACCGGGTGCTGTAACTGaacaac                                                                >  8:35681/1‑148 (MQ=255)
                                                                                     tCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGa                                                                                                       <  1:415470/107‑1 (MQ=255)
                                                                                     tCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGa                                                                                                       >  2:415470/1‑107 (MQ=255)
                                                                                       aTTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGaa                                                           >  2:399518/1‑149 (MQ=255)
                                                                                            aCGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGGGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGc                                                      <  4:127505/149‑1 (MQ=255)
                                                                                               ggCCCTGAAAACCCTGGACCGCCTTCTTCATGAACCCGGCGTGGCCTTCAGCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGa                                                   <  1:248193/149‑1 (MQ=255)
                                                                                                   cTGAATACCCTGGACCGCGATCTTCATGAACCCGGCGTGGCCTTCACCTCGACCCGTGCATCACGGGTATGGCTGATTCTGAACATTGAGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGc                                               <  4:350074/149‑1 (MQ=255)
                                                                                                    tGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACCTGGCGGGACTGGAGCGTGAAGAAGTCCGGGTGCTGTATCGGACCACCCAGAAGCAGCTGATTGcc                                              >  7:193123/1‑149 (MQ=255)
                                                                                                    tGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGcc                                              >  2:394229/1‑149 (MQ=255)
                                                                                                       aaaCCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGt                                           <  8:188110/149‑1 (MQ=255)
                                                                                                           ccTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCg                                             <  7:318203/143‑1 (MQ=255)
                                                                                                           ccTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCg                                             >  8:318203/1‑143 (MQ=255)
                                                                                                           ccTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGaaa                                       >  1:310087/1‑149 (MQ=255)
                                                                                                             tGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAAcc                                     >  5:139373/1‑149 (MQ=255)
                                                                                                               gACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAg                                                       <  5:80888/129‑1 (MQ=255)
                                                                                                               gACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAg                                                       >  6:80888/1‑129 (MQ=255)
                                                                                                                 ccGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCctct                                 <  4:173340/149‑1 (MQ=255)
                                                                                                                   gCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTc                               <  3:168107/149‑1 (MQ=255)
                                                                                                                    ccATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGCGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCa                              >  7:377787/1‑149 (MQ=255)
                                                                                                                         ttCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCgg                          >  3:193692/1‑148 (MQ=255)
                                                                                                                          tCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGaacaa                                                                 >  4:239099/1‑108 (MQ=255)
                                                                                                                          tCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGaacaa                                                                 <  3:239099/108‑1 (MQ=255)
                                                                                                                           cATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCGGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCg                           <  2:231363/145‑1 (MQ=255)
                                                                                                                                aCCCGGCGTGGCCTTCACCTCCTCGCGTGCAGCTCGGGAATGGCTGATTCTGAACATGGCGGGACTGGATATTGAAGAATTACTGGTGCTGTATCTGAACAACCAGATTCAGCTGATTTACGGTGAAACACTCTTCACCGGCACCATCa                  <  6:184038/149‑1 (MQ=255)
                                                                                                                                 cccGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCaa                 >  1:251652/1‑149 (MQ=255)
                                                                                                                                     gCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCCACCg              >  1:300920/1‑148 (MQ=255)
                                                                                                                                           ccTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACgg         >  3:187044/1‑147 (MQ=255)
                                                                                                                                            cTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGACTTCCGGGTGCTGTATCTGAACAACCAGATTCAGCTGATTGCCGGTGAAAGCCTCTTCACCGGCACCATCAACCGCAGGGAAg      <  2:57242/149‑1 (MQ=255)
                                                                                                                                                aCCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGTGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTcc   <  2:52597/148‑1 (MQ=255)
                                                                                                                                                aCCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCa  >  3:233538/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCATTCAACGGGCCCTGAAAACCCTGGACCGCCATCTTCATGAACCCGGCGTGGCCTTCACCTCCACCCGTGCAGCACGGGAATGGCTGATTCTGAACATGGCGGGACTGGAGCGTGAAGAATTCCGGGTGCTGTATCTGAACAACCAGAATCAGCTGATTGCCGGTGAAACCCTCTTCACCGGCACCATCAACCGCACGGAAGTCCA  >  NC_000913/2076245‑2076537

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