Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,076,188 C→T 5.6% N470N (AAC→AAT yeeR → CP4‑44 prophage; putative membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,076,1880CT 5.6% 222.8 / 2.8 73N470N (AAC→AATyeeRCP4‑44 prophage; putative membrane protein
Reads supporting (aligned to +/- strand):  ref base C (35/33);  new base T (2/2);  total (37/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.20e-01

GACAGAGTCCGTCTGAAAAACGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCAT  >  NC_000913/2076042‑2076333
                                                                                                                                                  |                                                                                                                                                 
gACAGAGTCCGTCTGAAAAACGCACGGATCCTCGGTGATGATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACg                                                                                                                                                  >  7:364195/1‑148 (MQ=255)
gACAGAGTCCGTCTGAAAAACGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACgg                                                                                                                                                 >  8:301176/1‑149 (MQ=255)
gACAGAGTCCGTCTGAAAAACGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACg                                                                                                                                                  <  4:206414/148‑1 (MQ=255)
gACAGAGTCCGTCTGAAAAACGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACg                                                                                                                                                  >  3:206414/1‑148 (MQ=255)
 aCAGAGTCCGTCTGAAAAACGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGa                                                                                                                                                <  5:336726/149‑1 (MQ=255)
 aCAGAGTCCGTCTGAAAAACGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGa                                                                                                                                                <  8:261210/149‑1 (MQ=255)
   agagTCCGTCTGAAAAACGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGa                                                                                                                                                >  6:14514/1‑147 (MQ=255)
     agTCCGTCTGAAAAACGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGCCTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAg                                                                                                                                             >  8:270846/1‑148 (MQ=255)
         cGTCTGAAAAACGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCTACGGACAAAACGGACAGTATc                                                                                                                                         <  7:270846/148‑1 (MQ=255)
         cGTCTGAAAAACGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATc                                                                                                                                         <  5:361326/148‑1 (MQ=255)
             tGAAAAACGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAATGGACAGTATCGTTa                                                                                                                                     >  1:144607/1‑148 (MQ=255)
             tGAAAAACGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAATGGACAGTATCGTTa                                                                                                                                     >  1:149720/1‑148 (MQ=255)
              gAAAAACGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACa                                                                                                                                   >  7:369580/1‑149 (MQ=255)
              gAAAAACGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACa                                                                                                                                   >  1:89850/1‑149 (MQ=255)
               aaaaaCGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTa                                                                                                                                     <  5:195536/146‑1 (MQ=255)
               aaaaaCGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGGGCGGCATTCGCCGGACAAAACGGACAGTATCGTTACAt                                                                                                                                  >  5:263593/1‑149 (MQ=255)
               aaaaaCGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGCAACAGAAATCCACACGAAATATTGATCAACTGCAGCCCGTAGCGTCGGTGCGCCATTCGACGCACAAAACGGACAGTATCGTTACAt                                                                                                                                  >  5:337521/1‑149 (MQ=255)
                 aaaCGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCTTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATgg                                                                                                                                >  8:328679/1‑149 (MQ=255)
                      cACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGaaa                                                                                                                            >  4:267352/1‑148 (MQ=255)
                      cACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGACCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGaaa                                                                                                                            >  4:272020/1‑148 (MQ=255)
                       aCGGATCCTCGGTGATGATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAaata                                                                                                                          <  8:364195/149‑1 (MQ=255)
                        cGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAaataa                                                                                                                         <  6:63090/149‑1 (MQ=255)
                         ggATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAaataat                                                                                                                        <  7:336482/149‑1 (MQ=255)
                         ggATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAaataat                                                                                                                        <  3:272020/149‑1 (MQ=255)
                         ggATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAaataat                                                                                                                        <  3:267352/149‑1 (MQ=255)
                          gATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATg                                                                                                                       <  5:14514/149‑1 (MQ=255)
                           aTCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAaataa                                                                                                                         >  7:164026/1‑146 (MQ=255)
                           aTCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAaataa                                                                                                                         <  8:164026/146‑1 (MQ=255)
                                       aataatGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCACCt                                                                                                          >  7:287059/1‑149 (MQ=255)
                                         taatGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTgg                                                                                                        <  1:130408/149‑1 (MQ=255)
                                         taatGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATCGCCCCATGCAACTg                                                                                                         >  8:331228/1‑148 (MQ=255)
                                          aatGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTgg                                                                                                        >  8:93603/1‑148 (MQ=255)
                                           atGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGa                                                                                                       <  2:89850/148‑1 (MQ=255)
                                            tGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAg                                                                                                     >  3:87497/1‑149 (MQ=255)
                                                aGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTa                                                                                                                                     >  8:241458/1‑113 (MQ=255)
                                                aGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTa                                                                                                                                     <  7:241458/113‑1 (MQ=255)
                                                     tGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACTAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAATGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTg                                                                                            <  2:144607/149‑1 (MQ=255)
                                                     tGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACTAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAATGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTg                                                                                            <  2:149720/149‑1 (MQ=255)
                                                       taGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAATACTGAGAGTATCGTTACATGGGAAATAATGGCGCCATGCAACTGGAAGTCCCCCGTGa                                                                                           <  4:87092/147‑1 (MQ=255)
                                                       gAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGAt                                                                                          <  6:308914/149‑1 (MQ=255)
                                                         gCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTa                                                                                                                                           >  6:296714/1‑98 (MQ=255)
                                                               agacTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGc                                                                                  <  7:301176/149‑1 (MQ=255)
                                                               agacTGGTCAGCGGAACAGAAATACAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTAGTCGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCTACTGTAAGTCCCCCGTGATCAGTATGc                                                                                  <  4:209218/149‑1 (MQ=255)
                                                                gacTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGcc                                                                                 >  7:124109/1‑149 (MQ=255)
                                                                gacTGCTCAGCGGAACAGAAATCCAGCCGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGAGCGGCATTCGACGGACAAAACGGAAAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGcc                                                                                 >  8:345648/1‑149 (MQ=255)
                                                                              aCAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAAc                                                                   <  7:328679/149‑1 (MQ=255)
                                                                                 gAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCCTGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAAccct                                                                >  3:81829/1‑147 (MQ=255)
                                                                                 gAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCAt                                                                <  6:71102/149‑1 (MQ=255)
                                                                                      ccAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAgg                                                            >  7:250856/1‑148 (MQ=255)
                                                                                       cAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAAACATGAGGaa                                                          >  5:324199/1‑149 (MQ=255)
                                                                                        aGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGaa                                                          >  7:301132/1‑148 (MQ=255)
                                                                                        aGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAAt                                                         <  8:124109/149‑1 (MQ=255)
                                                                                            gAAATTTTGTTCAACTGCAGGCCGTAGCGTCGGTGCGGCATTCGACGGACAATACGGAGAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCATGATCAGTATGCCGGCGCTGtt                                                                        <  6:69591/130‑1 (MQ=255)
                                                                                                    gTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGa                                              <  4:87497/148‑1 (MQ=255)
                                                                                                     ttCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAg                                            >  7:342816/1‑149 (MQ=255)
                                                                                                       cAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAACACGGACAGTATCGTTCCATGGGTAATAATGGCCCCAAGCAACTGGAAGTCCCCCGTGATATGTATGCCGGAGCGGTTGAAACCATGAGGAATAAGATCCGCGAAgg                                           >  5:138414/1‑148 (MQ=255)
                                                                                                        aaCTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGt                                          >  7:168709/1‑148 (MQ=255)
                                                                                                          cTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTa                                         >  1:270817/1‑147 (MQ=255)
                                                                                                           tGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAg                                      >  5:305395/1‑149 (MQ=255)
                                                                                                            gCAGGCCGTAGCGTCGGTGCGTCATTCGACGGACATCACGGACAGTATCGTTACATGGGAATTAATCGCCCGTTGCAACTGGAAGTCCCCCGTGTTCAGTATGCGGGCTCTGTTGAAAGCATGAGGAATAAGTTCCGCGAAGGTTAAg                                      <  6:337521/148‑1 (MQ=255)
                                                                                                             cAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAg                                      <  8:250856/147‑1 (MQ=255)
                                                                                                             cAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGAATCGTTACATGGTAAATAATGGCCCCATGCAACTGGTAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAg                                      <  1:20793/147‑1 (MQ=255)
                                                                                                                  cGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAgagga                               <  8:369580/149‑1 (MQ=255)
                                                                                                                      gCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGgagag                            <  8:168709/148‑1 (MQ=255)
                                                                                                                       cctCGGTGCGGCCTTCGACGTACAAAACGGACAGGAGCGTTACAGGGTAAATAATGGCCCCATGCAACTGGAAGTCCCCCGGGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGAt                          <  1:264194/147‑1 (MQ=255)
                                                                                                                         gcGGTGCGGCATTCGACGGACAAAACGGACAGTTTCGTTACATGGGATATAATGGCCCCATGCTTGTGGTAGTCCGCGGTGATCATTATGCCGGCGGTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGTTGa                        <  7:345648/148‑1 (MQ=255)
                                                                                                                           ggTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACg                      <  2:270817/149‑1 (MQ=255)
                                                                                                                                 gCATTCGACGGACAAAACGGACAGTATCGTTACATGTGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACAATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGAAGACCGGc                >  7:221989/1‑149 (MQ=255)
                                                                                                                                    ttCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAg             >  5:299676/1‑149 (MQ=255)
                                                                                                                                     tcgtcgGACTAAAGGGACAGTATCGTTACTTGGGAAATTATGGCCCCATGCAACTGGATGTCCGCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGGGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGa              <  7:331228/147‑1 (MQ=255)
                                                                                                                                       gACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATCCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAgcg           >  2:47637/1‑148 (MQ=255)
                                                                                                                                               aaaCGGACAGTATCGTTACATGGGAAGTAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCAt  >  7:22992/1‑149 (MQ=255)
                                                                                                                                               aaaCGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGCGATGACGCCCGGCGCGCGCAGCCTCAt  >  8:128503/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GACAGAGTCCGTCTGAAAAACGCACGGATCCTCGGTGATAATAATGCCAGAAATGGAGCAGACAGACTGGTCAGCGGAACAGAAATCCAGACGAAATATTGTTCAACTGCAGCCCGTAGCGTCGGTGCGGCATTCGACGGACAAAACGGACAGTATCGTTACATGGGAAATAATGGCCCCATGCAACTGGAAGTCCCCCGTGATCAGTATGCCGGCGCTGTTGAAACCATGAGGAATAAGATCCGCGAAGGTAAAGTAGAGGAGAGATGACGCCCGGCGAGCGCAGCCTCAT  >  NC_000913/2076042‑2076333

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: