Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,038,118 C→G 5.5% 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.5% 647.2 / 22.0 183L127V (CTG→GTG) cbdAcytochrome bd‑II oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base C (88/84);  new base G (5/5);  total (94/89)
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.24e-01

Only 100 of 193 total aligned reads displayed.

TACTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGCTGGTCTT  >  NC_000913/1037973‑1038266
                                                                                                                                                 |                                                                                                                                                    
tACTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACctgc                                                                                                                                                   <  1:1039485/149‑1 (MQ=255)
tACTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACctgc                                                                                                                                                   >  2:613530/1‑149 (MQ=255)
  cTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACctgctg                                                                                                                                                 >  3:19903/1‑149 (MQ=255)
  cTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACctgctg                                                                                                                                                 >  1:1016647/1‑149 (MQ=255)
  cTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACctgctg                                                                                                                                                 >  4:4005/1‑149 (MQ=255)
   tAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTgg                                                                                                                                                <  2:1263131/149‑1 (MQ=255)
      cTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGtat                                                                                                                                             <  1:861254/149‑3 (MQ=255)
       tGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGAc                                                                                                                                            <  4:350863/149‑1 (MQ=255)
       tGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCAGTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGAc                                                                                                                                            >  2:1304519/1‑149 (MQ=255)
         gTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGtattta                                                                                                                                          <  3:989555/149‑6 (MQ=255)
         gTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGt                                                                                                                                          <  2:847700/149‑1 (MQ=255)
         gTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGt                                                                                                                                          >  2:60141/1‑149 (MQ=255)
             ttCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCt                                                                                                                                      >  3:385651/1‑149 (MQ=255)
                  ttCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCtggtgg                                                                                                                                 >  3:886404/1‑149 (MQ=255)
                    ccAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGtat                                                                                                                                             <  4:1100681/135‑3 (MQ=255)
                     cAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAGTGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGt                                                                                                                              >  2:553995/1‑149 (MQ=255)
                      aaCTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGta                                                                                                                                              >  1:1323058/1‑131 (MQ=255)
                      aaCTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCAGGTGCTGGta                                                                                                                                              <  2:1323058/132‑2 (MQ=255)
                        cTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCg                                                                                                                           <  2:398835/149‑1 (MQ=255)
                         tATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCgg                                                                                                                          >  1:286460/1‑149 (MQ=255)
                          aTGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGt                                                                                                                         <  3:423366/149‑1 (MQ=255)
                           tgtgGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGtt                                                                                                                        <  3:287347/149‑1 (MQ=255)
                               ggCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCaaa                                                                                                                    >  4:1300481/1‑149 (MQ=255)
                               ggCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCaaa                                                                                                                    <  2:1303491/149‑1 (MQ=255)
                               ggCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCaaa                                                                                                                    >  4:458848/1‑149 (MQ=255)
                                gCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAGTGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAt                                                                                                                   <  1:553995/149‑1 (MQ=255)
                                gCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAt                                                                                                                   >  1:1205662/1‑149 (MQ=255)
                                  gATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCtggtg                                                                                                                                  >  2:748931/1‑132 (MQ=255)
                                  gATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCtggtg                                                                                                                                  <  1:748931/132‑1 (MQ=255)
                                  gATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtct                                                                                                                 <  2:286460/149‑1 (MQ=255)
                                  gATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtct                                                                                                                 >  2:1391175/1‑149 (MQ=255)
                                     atTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAAtctctc                                                                                                              <  3:702866/149‑1 (MQ=255)
                                      tttttGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATctctct                                                                                                             <  3:458848/149‑1 (MQ=255)
                                             cgcACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGTTTGTTCTTCTTCGGCTGGTAACGGTGGAATAAATACCAGCCCCTGCTGGTGACGTGGCTGGTGGCGTGCGGTTCAACTCTCTCTGCGTtgt                                                                                                      <  2:665232/149‑1 (MQ=38)
                                             cgcACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTtgt                                                                                                      <  1:638426/149‑1 (MQ=255)
                                             cgcACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTtgt                                                                                                      <  4:558890/149‑1 (MQ=255)
                                              gcACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTtgtg                                                                                                     >  4:200711/1‑149 (MQ=255)
                                               cACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAGTCTCTCTGCGTTGTgg                                                                                                    <  1:1250142/149‑1 (MQ=255)
                                               cACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTgg                                                                                                    >  2:192505/1‑149 (MQ=255)
                                                aCCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGa                                                                                                   >  4:61632/1‑149 (MQ=255)
                                                aCCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGa                                                                                                   >  4:73665/1‑149 (MQ=255)
                                                  cGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGata                                                                                                 <  3:1300481/149‑1 (MQ=255)
                                                       gCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGaa                                                                                            <  1:1294494/149‑1 (MQ=255)
                                                        cGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATGCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAAt                                                                                           <  2:1472595/149‑1 (MQ=255)
                                                        cGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAAt                                                                                           <  2:680652/149‑1 (MQ=255)
                                                        cGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAAt                                                                                           <  2:1016647/149‑1 (MQ=255)
                                                                gCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATg                                                                                          <  4:664009/142‑1 (MQ=255)
                                                                gCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATg                                                                                          >  3:664009/1‑142 (MQ=255)
                                                                gCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACgg                                                                                   <  4:554042/149‑1 (MQ=255)
                                                                  aTTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGtt                                                                                 <  1:187057/149‑1 (MQ=255)
                                                                      aTGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGAt                                                                             >  1:1405756/1‑149 (MQ=255)
                                                                          ccTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCaa                                                                         >  2:4941/1‑149 (MQ=255)
                                                                                  cTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGa                                                                                                   <  3:510862/115‑1 (MQ=255)
                                                                                  cTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGAc                                                                 <  1:368147/149‑1 (MQ=255)
                                                                                   tCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGtt                                                                                 <  4:169938/132‑1 (MQ=255)
                                                                                   tCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGAcc                                                                >  2:1284607/1‑149 (MQ=255)
                                                                                    cGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCGGACCg                                                               <  4:1351977/149‑1 (MQ=255)
                                                                                    cGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGCTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCg                                                               >  3:1212732/1‑149 (MQ=255)
                                                                                      aaTCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGt                                                             <  3:61667/149‑1 (MQ=255)
                                                                                           aCCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCa                                                        >  2:869648/1‑149 (MQ=255)
                                                                                             cTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCAtt                                                      >  1:414432/1‑149 (MQ=255)
                                                                                             cTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCAtt                                                      >  4:52395/1‑149 (MQ=255)
                                                                                              tttGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCAttt                                                     <  4:1323597/149‑1 (MQ=255)
                                                                                                  tCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCCCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGat                                                 <  2:701791/149‑1 (MQ=255)
                                                                                                      gCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCg                                             >  2:969724/1‑149 (MQ=255)
                                                                                                          ttcttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGacac                                         >  3:551267/1‑149 (MQ=255)
                                                                                                          ttcttcttcGGCTGGCAACGGCTGAAGAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTgc                                                           <  3:928836/131‑1 (MQ=255)
                                                                                                             ttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCt                                      >  4:277650/1‑149 (MQ=255)
                                                                                                             ttcttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCt                                      >  1:385989/1‑149 (MQ=255)
                                                                                                               cttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGc                                    <  3:50689/149‑1 (MQ=255)
                                                                                                               cttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGc                                    <  4:705097/149‑1 (MQ=255)
                                                                                                               cttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGc                                    >  2:709479/1‑149 (MQ=255)
                                                                                                               cttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGc                                    >  4:1087056/1‑149 (MQ=255)
                                                                                                                ttcGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCg                                   >  4:1439690/1‑149 (MQ=255)
                                                                                                                     cTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATgg                              >  3:258873/1‑149 (MQ=255)
                                                                                                                      tGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGga                             <  3:1087056/149‑1 (MQ=255)
                                                                                                                      tGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGga                             <  1:1391175/149‑1 (MQ=255)
                                                                                                                      tGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGga                             <  1:402084/149‑1 (MQ=255)
                                                                                                                        gCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGgaga                           <  1:969724/149‑1 (MQ=255)
                                                                                                                         cAACGGCTGAATAAATACCAGCACCTGTTGGTGACGTGGCTGGTGGCGTTCGGTTAAAATCTCTCTGCGTTGTGAATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGAt                          <  1:1304519/149‑1 (MQ=255)
                                                                                                                          aaCGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATg                         >  2:1156566/1‑149 (MQ=255)
                                                                                                                           aCGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGa                        <  1:307454/149‑1 (MQ=255)
                                                                                                                            cGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGAc                       <  3:1382434/149‑1 (MQ=255)
                                                                                                                              gCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGa                        <  3:794291/146‑1 (MQ=255)
                                                                                                                              gCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGa                        >  4:794291/1‑146 (MQ=255)
                                                                                                                                tGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGc                   >  1:1444402/1‑149 (MQ=255)
                                                                                                                                 gAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCt                  >  4:1076394/1‑149 (MQ=255)
                                                                                                                                  aataaataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCtt                 >  3:1033008/1‑149 (MQ=255)
                                                                                                                                     aaataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAg              >  1:270222/1‑149 (MQ=255)
                                                                                                                                      aataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGAt                                                                             >  3:708943/1‑85 (MQ=255)
                                                                                                                                      aataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGc             >  4:1148462/1‑149 (MQ=255)
                                                                                                                                       tcaCCAGCAGGTGCTGGTGACGTGGGTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCg            >  4:915452/3‑149 (MQ=255)
                                                                                                                                       tcaCCAGCAGGTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGAt                                                                             >  2:30572/3‑84 (MQ=255)
                                                                                                                                       ataCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCg            <  4:258873/149‑1 (MQ=255)
                                                                                                                                         aCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATctctct                                                                                                             <  4:289226/50‑1 (MQ=255)
                                                                                                                                         aCCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATctctct                                                                                                             >  3:289226/1‑50 (MQ=255)
                                                                                                                                          ccAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGc         <  2:1205662/149‑1 (MQ=255)
                                                                                                                                              cACCTGCTGGTGAAGTGGCTGGTGGCGTTCGGTTAAAATCTCTCTCCGTTTTGGATAATTCATGCCAACGGTTGGATGCAAAACCCGACCGGTCCGCATTTTGATATCGACACACTGCGTATGGAAATGACCAGATTGAGCGGGCTGGt     <  4:447594/149‑1 (MQ=255)
                                                                                                                                                ccTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGCTGGTCt   <  3:1125936/149‑1 (MQ=255)
                                                                                                                                                ccTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGCTGGTCt   <  3:52395/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
TACTAACTGGTCATTCTATTCCAACTATGTGGGCGATATTTTTGGCGCACCGCTGGCGATGGAAGCATTAATGGCCTTCTTCCTCGAATCCACCTTTGTCGGGCTGTTCTTCTTCGGCTGGCAACGGCTGAATAAATACCAGCACCTGCTGGTGACGTGGCTGGTGGCGTTCGGTTCAAATCTCTCTGCGTTGTGGATATTGAATGCCAACGGTTGGATGCAATACCCGACCGGTGCGCATTTTGATATCGACACCCTGCGTATGGAGATGACCAGCTTCAGCGAGCTGGTCTT  >  NC_000913/1037973‑1038266

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