Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,038,118 C→G 5.6% L127V (CTG→GTG)  cbdA → cytochrome bd‑II oxidase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,038,1180CG 5.6% 324.5 / 9.0 90L127V (CTG→GTG) cbdAcytochrome bd‑II oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base C (44/41);  new base G (3/2);  total (47/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.42e-01

TACTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGCTGGTC  >  NC_000913/1037973‑1038264
                                                                                                                                                 |                                                                                                                                                  
tACTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACctgc                                                                                                                                                 >  3:301082/1‑149 (MQ=255)
  cTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACctgctg                                                                                                                                               >  2:643859/1‑149 (MQ=255)
     aCTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTg                                                                                                                                            >  4:629289/1‑149 (MQ=255)
        ggTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGtattt                                                                                                                                         <  1:120838/149‑5 (MQ=255)
         gTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGt                                                                                                                                        <  1:167418/149‑1 (MQ=255)
         gTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGt                                                                                                                                        <  4:301082/149‑1 (MQ=255)
         gTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGt                                                                                                                                        <  2:54482/149‑1 (MQ=255)
         gTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGt                                                                                                                                        <  4:344947/149‑1 (MQ=255)
            aTTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACg                                                                                                                                         >  1:434378/1‑145 (MQ=255)
            aTTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGc                                                                                                                                     >  3:585414/1‑149 (MQ=255)
            aTTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCCGGTGACg                                                                                                                                         <  2:434378/145‑1 (MQ=255)
               cTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCtgg                                                                                                                                  <  4:293051/149‑1 (MQ=255)
                 aTTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGTtggtg                                                                                                                                >  2:429102/1‑149 (MQ=255)
                 aTTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCtggtg                                                                                                                                <  1:599194/149‑1 (MQ=255)
                  ttCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATAGCAGCACCTGCTGGTGACGTGGCtggtgg                                                                                                                               >  4:513326/1‑149 (MQ=255)
                  ttCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAcc                                                                                                                                                    <  2:266859/128‑1 (MQ=255)
                  ttCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAcc                                                                                                                                                    >  1:266859/1‑128 (MQ=255)
                   tCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTTGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGc                                                                                                                              <  2:109555/149‑1 (MQ=255)
                   tCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGc                                                                                                                              <  4:649800/149‑1 (MQ=255)
                     cAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGt                                                                                                                            >  4:608974/1‑149 (MQ=255)
                     cAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGt                                                                                                                            >  2:725460/1‑149 (MQ=255)
                      aaCTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGtt                                                                                                                           <  2:137009/149‑1 (MQ=255)
                      aaCTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGtt                                                                                                                           >  2:45854/1‑149 (MQ=255)
                         tATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGTTGGTGGCGTTCgg                                                                                                                        <  1:429102/149‑1 (MQ=255)
                         tATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCgg                                                                                                                        >  3:570146/1‑149 (MQ=255)
                             tgGGCGATGTTTTTGGCGCACCGCTGGCAATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCa                                                                                                                    >  4:305725/1‑149 (MQ=255)
                                  gATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtct                                                                                                               >  3:283877/1‑149 (MQ=255)
                                        tttGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGc                                                                                                         <  1:266510/149‑1 (MQ=255)
                                                  cGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGata                                                                                               <  3:313922/149‑1 (MQ=255)
                                                      ggCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGa                                                                                           <  2:163979/149‑1 (MQ=255)
                                                         gATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATg                                                                                        >  4:567441/1‑149 (MQ=255)
                                                           tGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGcc                                                                                      <  4:708621/149‑1 (MQ=255)
                                                            ggAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCa                                                                                     >  2:612294/1‑149 (MQ=255)
                                                             gAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCaa                                                                                    <  1:143218/149‑1 (MQ=255)
                                                             gAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCaa                                                                                    >  2:143218/1‑149 (MQ=255)
                                                              aaGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAAc                                                                                   >  3:300993/1‑149 (MQ=255)
                                                              aaGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAAc                                                                                   >  3:207355/1‑149 (MQ=255)
                                                                  aTTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGtt                                                                               <  3:386946/149‑1 (MQ=255)
                                                                  aTTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGtt                                                                               <  4:283877/149‑1 (MQ=255)
                                                                  aTTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGtt                                                                               >  3:586291/1‑149 (MQ=255)
                                                                  aTTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGtt                                                                               <  1:643859/149‑1 (MQ=255)
                                                                     aaTGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGa                                                                            >  1:332614/1‑149 (MQ=255)
                                                                     aaTGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGa                                                                            <  2:140868/149‑1 (MQ=255)
                                                                      aTGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGAt                                                                           >  3:527900/1‑149 (MQ=255)
                                                                      aTGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGAt                                                                           <  3:770082/149‑1 (MQ=255)
                                                                       tGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATg                                                                          <  2:721293/149‑1 (MQ=255)
                                                                        ggCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGc                                                                         <  3:629289/149‑1 (MQ=255)
                                                                         gCCTTCTTCCTTGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCa                                                                        >  3:633210/1‑149 (MQ=255)
                                                                           cttcttCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAAt                                                                      >  2:420877/1‑149 (MQ=255)
                                                                             tcttcCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATAc                                                                    <  4:207355/149‑1 (MQ=255)
                                                                                 ccTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGa                                                                <  2:170513/149‑1 (MQ=255)
                                                                                   tCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGAcc                                                              >  1:722072/1‑149 (MQ=255)
                                                                                   tCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGAcc                                                              >  2:558734/1‑149 (MQ=255)
                                                                                        tCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGAt                                                                           >  4:440400/1‑131 (MQ=255)
                                                                                        tCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGAt                                                                           <  3:440400/131‑1 (MQ=255)
                                                                                           aCCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGt                                                                                >  1:551278/1‑123 (MQ=255)
                                                                                           aCCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGt                                                                                <  2:551278/123‑1 (MQ=255)
                                                                                           aCCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCa                                                      >  2:328873/1‑149 (MQ=255)
                                                                                              tttGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCAttt                                                   >  1:334277/1‑149 (MQ=255)
                                                                                                    gggCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGatat                                             >  3:437154/1‑149 (MQ=255)
                                                                                                     ggCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATc                                            <  4:680484/149‑1 (MQ=255)
                                                                                                          ttcttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGacac                                       >  3:119931/1‑149 (MQ=255)
                                                                                                            cttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGa                                                                            <  1:158206/110‑1 (MQ=255)
                                                                                                            cttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGa                                                                            >  2:158206/1‑110 (MQ=255)
                                                                                                             ttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCt                                    >  4:727580/1‑149 (MQ=255)
                                                                                                             ttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCt                                    >  1:458406/1‑149 (MQ=255)
                                                                                                               cttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGc                                  >  3:738308/1‑149 (MQ=255)
                                                                                                                ttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCg                                 >  1:572782/1‑149 (MQ=255)
                                                                                                                    gCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATg                             <  4:633210/149‑1 (MQ=255)
                                                                                                                      tGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGga                           <  1:612294/149‑1 (MQ=255)
                                                                                                                        gCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGgaga                         <  2:579530/149‑1 (MQ=255)
                                                                                                                          aaCGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATg                       >  3:140241/1‑149 (MQ=255)
                                                                                                                            cGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGAc                     <  1:328873/149‑1 (MQ=255)
                                                                                                                             ggCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGAcc                    <  4:300993/149‑1 (MQ=255)
                                                                                                                               cTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAg                  <  4:437154/149‑1 (MQ=255)
                                                                                                                                 gAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCt                >  2:229393/1‑149 (MQ=255)
                                                                                                                                 gAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCt                >  4:684099/1‑149 (MQ=255)
                                                                                                                                  aataaataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCtt               <  1:229393/149‑1 (MQ=255)
                                                                                                                                     aaataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAg            <  4:140241/149‑1 (MQ=255)
                                                                                                                                      aataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGc           >  4:434737/1‑149 (MQ=255)
                                                                                                                                      aataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGc           <  4:527900/149‑1 (MQ=255)
                                                                                                                                       tcaCCAGCAGGTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCg          >  2:649891/3‑149 (MQ=255)
                                                                                                                                       tcaCCAGCAGGTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCg          >  4:167097/3‑149 (MQ=255)
                                                                                                                                       ataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCg          >  2:404803/1‑149 (MQ=255)
                                                                                                                                        taCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGa         >  4:71666/1‑149 (MQ=255)
                                                                                                                                        taCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGa         >  3:340261/1‑149 (MQ=255)
                                                                                                                                         aCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAg        >  2:718143/1‑149 (MQ=255)
                                                                                                                                          ccAGCAGGTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGat                                               <  3:71998/109‑1 (MQ=255)
                                                                                                                                          ccAGCAGGTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGat                                               >  4:71998/1‑109 (MQ=255)
                                                                                                                                          ccAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGc       <  4:738308/149‑1 (MQ=255)
                                                                                                                                            aGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGCTg     >  1:569591/1‑149 (MQ=255)
                                                                                                                                             gCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGCTgg    >  1:77276/1‑149 (MQ=255)
                                                                                                                                               aCCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGCTGGTc  <  1:530310/149‑1 (MQ=255)
                                                                                                                                               aCCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGCTGGTc  <  3:567441/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
TACTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGCTGGTC  >  NC_000913/1037973‑1038264

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