Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,016,800 A→G 5.7% K247K (AAA→AAG ygfK → putative Fe‑S subunit oxidoreductase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,016,8000AG 5.7% 214.1 / 2.9 70K247K (AAA→AAGygfKputative Fe‑S subunit oxidoreductase subunit
Reads supporting (aligned to +/- strand):  ref base A (32/34);  new base G (2/2);  total (34/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 = 6.77e-01

GCCACGGATTGCAGGAAAAACGCGAAAGCTTGCAAGCCTTACCCGCTCGCATCCCCACCAGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATGC  >  NC_000913/3016652‑3016948
                                                                                                                                                    |                                                                                                                                                    
gCCACGGATTGCAGGAAAAACGCGAAAGCTTGCAAGCCTTACCCGCTCGCATCCCCACCAGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCTAAGCCATTTGCCACTACATGCTGGAAGaaaaa                                                                                                                                                      <  3:163678/149‑1 (MQ=255)
 ccACGGATTGCAGGAAAAACGCGAAAGCTTGCAAGCCTTACCCGCTCGCATCCCCACCAGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAg                                                                                                                                                     <  2:62272/149‑1 (MQ=255)
  cACGGATTGCAGGAAAAACGCGAAAGCTTGCAAGCCTTACCCGCTCGCATCCCCACCAGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAgg                                                                                                                                                    >  4:121285/1‑149 (MQ=255)
             ggAAAAACGCGAAAGCTTGCAAGCCTTACCCGCTCGCATCCCCACCAGTATGGTGCATGCCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCAATTTGCAGCTACATGCTGGAAGAAAAAGGGCTCAACACCt                                                                                                                                         >  4:170230/1‑149 (MQ=255)
                  aaCGCGAAAGCTTGCAAGCCTTACCCGCTCGCATCCCCACAAGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAACTCGAAGCCATTTGCCGCTCCATGCTGGAAGAAAAAGGGCTCCACACCTTtgt                                                                                                                                     >  5:339488/1‑148 (MQ=255)
                  aaCGCGAAAGCTTGCAAGACTTACGCGCTCGCATCCCCACCAGTATTGTGCATGGCGTCACGCTCTCCACCATGCACGGCTGTCGTCCGCATGAAATCGAAGCGATTTGCCGCTACATGCTGTAAGAAAAATGGCTCAACACCTTtgtg                                                                                                                                    <  7:82629/149‑1 (MQ=255)
                     gcgAAAGCTTGCAAGCCTTACCCGCTCGCATCCCCCCCAGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAACTCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGaaa                                                                                                                                 <  6:290347/149‑1 (MQ=255)
                       gAAAGCTTGCAAGCCTTACCCGCTCGCATCCCCACCAGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTc                                                                                                                                                >  2:271486/1‑132 (MQ=255)
                       gAAAGCTTGCAAGCCTTACCCGCTCGCATCCCCACCAGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTc                                                                                                                                                <  1:271486/132‑1 (MQ=255)
                          aGCTTGCAAGCCTTACCCGCTCGCATCCCCACCAGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTaa                                                                                                                            >  2:114790/1‑149 (MQ=255)
                               gCAAGCCTTACCCGCTCGCATCCCCACCAGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGa                                                                                                                       >  4:317493/1‑149 (MQ=255)
                               gCAAGCCTTACCCGCTCGCATCCCCACCAGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGa                                                                                                                       >  7:223331/1‑149 (MQ=255)
                                cAAGCCTTACCCGCTCGCATCCCCACCAGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGAc                                                                                                                      <  7:62743/149‑1 (MQ=255)
                                        aCCCGCTCGCATCCCCACCAGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACt                                                                                                                <  4:223679/147‑1 (MQ=255)
                                        aCCCGCTCGCATCCCCACCAGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTg                                                                                                               <  2:229644/148‑1 (MQ=255)
                                        aCCCGCTCGCATCCCCACCAGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTg                                                                                                               <  7:115516/148‑1 (MQ=255)
                                        aCCCGCTCGCATCCCCACCAGTATGGGGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTg                                                                                                               <  8:329218/148‑1 (MQ=255)
                                         cccGCTCGCATCCCCACCAGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTggg                                                                                                             <  5:166/149‑1 (MQ=255)
                                                   tCCCCACCAGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCgt                                                                                                   >  2:1729/1‑149 (MQ=255)
                                                      ccaccaGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGATCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGtt                                                                                                >  6:330137/1‑149 (MQ=255)
                                                       caccaGTATGGTGAATGGCGTCAGCCTCTCCACCATGCACGGCTGTACTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTc                                                                                               <  2:271497/149‑1 (MQ=255)
                                                           aGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAGGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTAcgcgc                                                                                                     >  3:257089/1‑139 (MQ=255)
                                                           aGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAGGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTAcgcgc                                                                                                     <  4:257089/139‑1 (MQ=255)
                                                            gTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTgag                                                                                          <  8:223331/149‑1 (MQ=255)
                                                               tGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATc                                                                                       <  4:96671/149‑1 (MQ=255)
                                                                ggTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAACAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATcc                                                                                      >  1:53007/1‑149 (MQ=255)
                                                                 gTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTGTGTGAAACTTAACCCGACAGTCCTGGGGTACGCGCGTGTTAGTGAGATCCt                                                                                     >  3:227408/1‑149 (MQ=255)
                                                                   gCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTg                                                                                                               >  8:174184/1‑121 (MQ=255)
                                                                   gCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTg                                                                                                               <  7:174184/121‑1 (MQ=255)
                                                                    cATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTTCATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCg                                                                                   <  3:170230/148‑1 (MQ=255)
                                                                        gCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGt                                                                               <  8:123734/148‑1 (MQ=255)
                                                                            cACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCg                                                                          >  5:260272/1‑149 (MQ=255)
                                                                              cccTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTACCCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGCTCCTCGATGTCTGCGGt                                                                        >  5:322177/1‑149 (MQ=255)
                                                                                ctctCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGttt                                                                      <  6:260272/149‑1 (MQ=255)
                                                                                ctctCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGttt                                                                      >  7:81344/1‑149 (MQ=255)
                                                                                 tctcCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAGGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTc                                                                     >  8:202556/1‑149 (MQ=255)
                                                                                 tctcCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTc                                                                     >  8:106343/1‑149 (MQ=255)
                                                                                 tctcCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTc                                                                     >  6:159522/1‑149 (MQ=255)
                                                                                   tcCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTc                                                                     <  4:3028/147‑1 (MQ=255)
                                                                                   tcCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCgg                                                                   >  5:207694/1‑149 (MQ=255)
                                                                                   tcCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCgg                                                                   >  1:175645/1‑149 (MQ=255)
                                                                                     ctccaTGCACGGATGTCCTACGCATGAAATCGAAGCTATTTGCCGCTACATGCTGTAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGt                                                                  <  7:6484/146‑1 (MQ=255)
                                                                                     cacctTGCACGGCTGTCCTCCGCAGGAAAGCGCAGACATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACGTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGGGGTCGGGAGATCCTCGATGTCTGCGGGTTCGGt                                                                  <  4:94934/148‑1 (MQ=255)
                                                                                     caccaTGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCg                                                                    <  5:330137/146‑1 (MQ=255)
                                                                                     caccaTGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGtt                                                                 <  5:122446/149‑1 (MQ=255)
                                                                                        caTGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACa                                                              >  8:47904/1‑149 (MQ=255)
                                                                                          tGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACa                                                              >  1:316845/1‑147 (MQ=255)
                                                                                            cACGGCTGTCCGCCGCATGAAATCGAAGCCTTTTGCCGCTACATGCGGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTGTCGGTTACATa                                                            <  2:64358/147‑1 (MQ=255)
                                                                                                   gTCCTCCGCATGAAATCCAAGCAATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTACCCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGAAGTCTGAGGTTTCGGTTACATAGGCTTaa                                                     >  2:186538/1‑147 (MQ=255)
                                                                                                    tcctccGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAaag                                                  <  1:1729/149‑1 (MQ=255)
                                                                                                     cctccGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAaaga                                                 <  7:109220/149‑1 (MQ=255)
                                                                                                         cGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAagaaga                                              <  3:317493/148‑1 (MQ=255)
                                                                                                           cATGAAATCGAAGCCATTTGCCGCTACATGCTGGAATAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCTACCTTACTGTGGTACGCGCGTGGTCGTGAGATCCTCGATGTCTGCGGGTTCGGTTACATAGGCTTAAAAGAAGAGTc                                           <  6:339488/149‑1 (MQ=255)
                                                                                                            aTGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCa                                          >  8:125467/1‑149 (MQ=255)
                                                                                                                    gAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCa                                  <  5:298314/149‑1 (MQ=255)
                                                                                                                    gAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGACTGCGGTTTAGGTTAAATAGGCTTAAATGAAGAGTCATTTGAGCa                                  >  4:183932/1‑149 (MQ=255)
                                                                                                                     aaGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCAc                                 <  5:268219/149‑1 (MQ=255)
                                                                                                                          aTTTGCCGCTACATGCTGGATGAAAAAGGGATCAACACCTTTGTGAAACTTAACACTAGCTTAATGGTGTTCGCGCGTGTTCGTGAGATGCTCGATGTCTGCGGTTTCGGTTACATAGTCTTAAAAGAAGAGTCATTTGATCACGAGCt                            >  4:70629/1‑149 (MQ=255)
                                                                                                                          aTTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTACCCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTGTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCt                            >  6:76300/1‑149 (MQ=255)
                                                                                                                          aTTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCt                            >  2:44598/1‑149 (MQ=255)
                                                                                                                           tttGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGt                                                                                             <  4:114222/83‑1 (MQ=255)
                                                                                                                           tttGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGt                                                                                             >  3:114222/1‑83 (MQ=255)
                                                                                                                           tttGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCt                            >  7:166246/1‑148 (MQ=255)
                                                                                                                               ccGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGc                       <  6:207694/149‑1 (MQ=255)
                                                                                                                                 gCTACATGCTGGAAGAAAAAGGGCAAAACACCTTTGTGAAACATAACCCGACAGTACTGGGGTACGCGCTTGTTCGTGAGATCCTAGATGTCTGCGGTTTCTGTTACATATTCTTAAATCAATATTCATTTTATAACGACCTCAATCt                      >  8:76782/1‑148 (MQ=255)
                                                                                                                                         cTGGAAGAAAAGGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAg              <  7:202556/148‑1 (MQ=255)
                                                                                                                                         cTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAg              <  8:81344/148‑1 (MQ=255)
                                                                                                                                            gaagaaAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACt          <  8:166246/149‑1 (MQ=255)
                                                                                                                                            gaagaaAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACt          >  6:143687/1‑149 (MQ=255)
                                                                                                                                              agaaAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGtt                                                                 >  3:216780/1‑92 (MQ=255)
                                                                                                                                              agaaAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGtt                                                                 <  4:216780/92‑1 (MQ=255)
                                                                                                                                                 aaaaGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGaaa     >  3:17424/1‑149 (MQ=255)
                                                                                                                                                 aaaaGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGaaa     >  4:211161/1‑149 (MQ=255)
                                                                                                                                                 aaaaGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGaaa     >  5:11181/1‑149 (MQ=255)
                                                                                                                                                   aaGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATg   >  1:283742/1‑149 (MQ=255)
                                                                                                                                                    aGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATGc  <  1:114790/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GCCACGGATTGCAGGAAAAACGCGAAAGCTTGCAAGCCTTACCCGCTCGCATCCCCACCAGTATGGTGCATGGCGTCACCCTCTCCACCATGCACGGCTGTCCTCCGCATGAAATCGAAGCCATTTGCCGCTACATGCTGGAAGAAAAAGGGCTCAACACCTTTGTGAAACTTAACCCGACCTTACTGGGGTACGCGCGTGTTCGTGAGATCCTCGATGTCTGCGGTTTCGGTTACATAGGCTTAAAAGAAGAGTCATTTGATCACGACCTCAAGCTGACGCAAGCACTGGAAATGC  >  NC_000913/3016652‑3016948

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: