Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,208,443 A→G 5.1% Q104Q (CAA→CAG tfaP → e14 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,208,4430AG 5.1% 246.6 / 3.0 78Q104Q (CAA→CAGtfaPe14 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (36/38);  new base G (2/2);  total (38/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 = 6.38e-01

AGAGGTTATCCGTACTCATGATGATGAAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGT  >  NC_000913/1208296‑1208590
                                                                                                                                                   |                                                                                                                                                   
agagGTTATCCATACTCATGATGATGAAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACaa                                                                                                                                                     >  6:219916/1‑148 (MQ=255)
 gagGTTATCCGTACTCATGATGATGAAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGACGATGAGCTGAAACAata                                                                                                                                                   >  2:330822/1‑149 (MQ=255)
  agGTTATCCGTACTCATGATGATGAAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAata                                                                                                                                                   >  1:17376/1‑148 (MQ=255)
   ggTTATCCGTACTCATGATGATGAAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAAtata                                                                                                                                                 <  1:392484/149‑1 (MQ=255)
    gTTATCCGTACTCATGATGATGAAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGCCAGTGGCAGGGTATGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATaa                                                                                                                                                >  2:286867/1‑149 (MQ=255)
        tCCGTACTCATGATGATGAAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAAttt                                                                                                                                             <  8:389021/148‑1 (MQ=255)
         gcgacttcatcatcaTGAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTtgtg                                                                                                                                          >  6:54751/7‑149 (MQ=255)
          cGTACTCATGATGATGAAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTtgt                                                                                                                                           <  1:164156/148‑1 (MQ=255)
          cGTACTCATGATGATGAAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTtgt                                                                                                                                           <  7:373244/148‑1 (MQ=255)
           gTACTCATGATGATGAAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTgg                                                                                                                                         >  3:295852/1‑149 (MQ=255)
            tACTCATGATGATGAAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGc                                                                                                                                        >  2:78910/1‑149 (MQ=255)
               tCATGATGATGAAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTg                                                                                                                                      >  2:347565/1‑148 (MQ=255)
                 atgatgatgAAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGa                                                                                                                                    <  3:257018/148‑1 (MQ=255)
                     tgatgaAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTTATGCAACTGATTTTATCAATAGCAGACAGTGTCAGGGTAAGGCTGCATTGGTAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATc                                                                                                                               <  5:184145/149‑1 (MQ=255)
                     tgatgaAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATc                                                                                                                               <  2:59110/149‑1 (MQ=255)
                      gatgaAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCt                                                                                                                              >  5:76609/1‑149 (MQ=255)
                        tgaAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTgg                                                                                                                            <  1:78910/149‑1 (MQ=255)
                         gaAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGTCTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTgg                                                                                                                            >  2:287498/1‑148 (MQ=255)
                               gCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGTCTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACt                                                                                                                     <  1:287498/149‑1 (MQ=255)
                               gCAGATGCGAACTTTAAAACGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATCGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACGTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACt                                                                                                                     <  4:337252/149‑1 (MQ=255)
                                 aGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTgg                                                                                                                   <  4:295852/149‑1 (MQ=255)
                                 aGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTgg                                                                                                                   <  5:13909/149‑1 (MQ=255)
                                   aTGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCGGGATTATCTGGAAGCACTGGaa                                                                                                                 <  5:305195/149‑1 (MQ=255)
                                        aaCTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGt                                                                                                            >  8:325246/1‑149 (MQ=255)
                                        aaCATTCAAAAGCAGATGCTTATCTCTGATGCAACAGAATTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGt                                                                                                            <  5:309011/149‑1 (MQ=255)
                                         aCTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGtt                                                                                                           <  1:339894/149‑1 (MQ=255)
                                             tCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGa                                                                                                         <  3:370999/147‑1 (MQ=255)
                                               aaaaGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACa                                                                                                     >  7:207117/1‑149 (MQ=255)
                                               aaaaGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACa                                                                                                     >  7:186719/1‑149 (MQ=255)
                                                 aaGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGTTTGATACATc                                                                                                   >  6:87108/1‑149 (MQ=255)
                                                 aaGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATc                                                                                                   >  6:142329/1‑149 (MQ=255)
                                                   gCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCa                                                                                                 <  8:135760/149‑1 (MQ=255)
                                                      gATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTg                                                                                              >  6:32384/1‑149 (MQ=255)
                                                           ttATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTgcgc                                                                                           >  2:313081/1‑147 (MQ=255)
                                                             aTCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCa                                                                                         <  3:32091/147‑1 (MQ=255)
                                                              tctcTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGa                                                                                       <  1:195258/148‑1 (MQ=255)
                                                                tctGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGatat                                                                                    >  4:267476/1‑149 (MQ=255)
                                                                tctGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGatat                                                                                    >  5:95211/1‑149 (MQ=255)
                                                                 ctGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATAtt                                                                                   >  5:341874/1‑149 (MQ=255)
                                                                    aTGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATAtt                                                                                   <  5:54751/146‑1 (MQ=255)
                                                                    aTGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGaa                                                                                <  8:186719/149‑1 (MQ=255)
                                                                       cAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGaa                                                                                                                 >  7:48937/1‑113 (MQ=255)
                                                                       cAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGaa                                                                                                                 <  8:48937/113‑1 (MQ=255)
                                                                       cAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATg                                                                              <  7:158160/148‑1 (MQ=255)
                                                                       cAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATg                                                                              <  5:219916/148‑1 (MQ=255)
                                                                        aaCTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGc                                                                            >  5:42166/1‑149 (MQ=255)
                                                                                  aTCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTc                                                                   <  1:368592/148‑1 (MQ=255)
                                                                                     aaTAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGCTATTGAATGGCCTACGCCTCCgg                                                                >  5:313331/1‑148 (MQ=255)
                                                                                           aGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCa                                                         <  8:207117/149‑1 (MQ=255)
                                                                                            gACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAg                                                        <  6:42166/149‑1 (MQ=255)
                                                                                                 tGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGc                                                      >  3:328456/1‑146 (MQ=255)
                                                                                                   gCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGa                                                  <  1:347565/148‑1 (MQ=255)
                                                                                                    cAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATg                                                <  6:76609/149‑1 (MQ=255)
                                                                                                       ggTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGTAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGACAGATGACa                                             >  2:166227/1‑149 (MQ=255)
                                                                                                       ggTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACa                                             >  4:64622/1‑149 (MQ=255)
                                                                                                         tAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATc                                           <  7:288342/149‑1 (MQ=255)
                                                                                                               tGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGgcgc                                     >  5:65205/1‑149 (MQ=255)
                                                                                                               tGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGgcgc                                     >  1:120061/1‑149 (MQ=255)
                                                                                                                 cATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGgcgc                                     <  7:325246/147‑1 (MQ=255)
                                                                                                                 cATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGgcgc                                     <  3:267476/147‑1 (MQ=255)
                                                                                                                  aTTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGCCTGGATTATCTGGAAGCACTGGAACAGGTTGATACATCCAGTGCGCCAGATATTGAATGGACTACGCCTCCGGCAGTTCACGCCAGATGACATCCGGCGCGGt                                  >  8:202627/1‑149 (MQ=255)
                                                                                                                      ggAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTg                              <  8:61867/149‑1 (MQ=255)
                                                                                                                       gAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCTGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTgg                             >  7:222430/1‑149 (MQ=255)
                                                                                                                            cTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCTGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCt                        <  8:222430/149‑1 (MQ=255)
                                                                                                                            cTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCt                        <  5:32384/149‑1 (MQ=255)
                                                                                                                            cTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCt                        >  8:272498/1‑149 (MQ=255)
                                                                                                                             ttAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCTGTGCTggattaca                       >  5:303094/1‑143 (MQ=255)
                                                                                                                             ttAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGAGCGGTGCTGGTATCTg                       >  7:369036/1‑149 (MQ=255)
                                                                                                                               aaaGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGtt                     >  6:293166/1‑149 (MQ=255)
                                                                                                                               aaaGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGtt                     >  8:278255/1‑149 (MQ=255)
                                                                                                                                   aagaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCg                 >  2:233794/1‑149 (MQ=255)
                                                                                                                                    agaTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGt                <  2:17376/149‑1 (MQ=255)
                                                                                                                                     gaTGAGCTGAAACAGTATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGa                                                                                       >  2:149635/1‑77 (MQ=255)
                                                                                                                                     gaTGAGCTGAAACAGTATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGa                                                                                       >  2:144056/1‑77 (MQ=255)
                                                                                                                                     gaTGAGCTGAAACAGTATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGa                                                                                       <  1:149635/77‑1 (MQ=255)
                                                                                                                                     gaTGAGCTGAAACAGTATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGa                                                                                       <  1:144056/77‑1 (MQ=255)
                                                                                                                                      aTGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCa              >  4:104934/1‑149 (MQ=255)
                                                                                                                                       tGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCAc             <  5:245727/149‑1 (MQ=255)
                                                                                                                                          gCTGAAACAATATAATTTGTGGCTGGCTTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCGGTTGCCGTCACcgc          <  6:341874/149‑1 (MQ=255)
                                                                                                                                            tGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGt        >  8:183403/1‑149 (MQ=255)
                                                                                                                                                aCAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAAt    >  3:139570/1‑149 (MQ=255)
                                                                                                                                                   atatAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGt  >  3:68712/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
AGAGGTTATCCGTACTCATGATGATGAAGTCGCAGATGCGAACTTTCAAAAGCAGATGCTTATCTCTGATGCAACTGATTTTATCAATAGCAGACAGTGGCAGGGTAAGGCTGCATTGGGAAGACTTAAAGAAGATGAGCTGAAACAATATAATTTGTGGCTGGATTATCTGGAAGCACTGGAACTGGTTGATACATCCAGTGCGCCAGATATTGAATGGCCTACGCCTCCGGCAGTTCAGGCCAGATGACATCCGGCGCGGTGCTGGTATCTGTTGCCGTCACCGCGTCAATGT  >  NC_000913/1208296‑1208590

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