Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 167,291 G→A 5.9% intergenic (+27/‑193) mrcB → / → fhuA fused glycosyl transferase and transpeptidase/ferrichrome outer membrane transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913167,2910GA 5.9% 329.6 / 7.3 102intergenic (+27/‑193)mrcB/fhuAfused glycosyl transferase and transpeptidase/ferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base G (51/45);  new base A (2/4);  total (53/49)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.25e-01

Only 100 of 104 total aligned reads displayed.

AATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAA  >  NC_000913/167145‑167437
                                                                                                                                                  |                                                                                                                                                  
aaTCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAg                                                                                                                                                  >  4:144051/1‑149 (MQ=255)
aaTCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAg                                                                                                                                                  >  2:169360/1‑149 (MQ=255)
 aTCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGt                                                                                                                                                 <  2:638440/149‑1 (MQ=255)
  tCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGtt                                                                                                                                                >  3:667461/1‑149 (MQ=255)
  tCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGtt                                                                                                                                                <  1:202382/149‑1 (MQ=255)
  tCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGtt                                                                                                                                                >  4:445811/1‑149 (MQ=255)
  tCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGtt                                                                                                                                                >  2:146066/1‑149 (MQ=255)
  tCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGtt                                                                                                                                                >  2:113791/1‑149 (MQ=255)
   ccGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCGCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTg                                                                                                                                               <  2:718558/149‑1 (MQ=255)
    cGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGc                                                                                                                                              <  2:421854/149‑1 (MQ=255)
     gTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGcc                                                                                                                                             >  1:747648/1‑149 (MQ=255)
      tttGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCa                                                                                                                                            >  1:138402/1‑149 (MQ=255)
         gATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAACGACGGCGACGGAGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATcc                                                                                                                                         <  2:637926/149‑1 (MQ=255)
          aTCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCg                                                                                                                                        >  1:483808/1‑149 (MQ=255)
           tCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAACTACCGGATGGCGAGTTGCCATCCgg                                                                                                                                       >  4:259852/1‑149 (MQ=255)
           tCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCgg                                                                                                                                       >  3:84929/1‑149 (MQ=255)
               tcttctCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATCCGGTatt                                                                                                                                   <  2:443386/149‑3 (MQ=255)
                     cAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACATCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACa                                                                                                                             <  1:241595/149‑1 (MQ=255)
                          gcagcaACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCa                                                                                                                        >  3:49095/1‑149 (MQ=255)
                            agcaACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATc                                                                                                                      <  1:265619/149‑1 (MQ=255)
                               aaCAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTaa                                                                                                                   <  4:42557/149‑1 (MQ=255)
                                 cAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGa                                                                                                                 <  2:663093/149‑1 (MQ=255)
                                  aGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGat                                                                                                                >  4:131579/1‑149 (MQ=255)
                                    ccGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGatat                                                                                                              >  4:619972/1‑149 (MQ=255)
                                       cAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTaa                                                                                                           >  3:196553/1‑149 (MQ=255)
                                       cAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTaa                                                                                                           >  2:663523/1‑149 (MQ=255)
                                           aGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCct                                                                                                       >  4:861663/1‑149 (MQ=255)
                                              aaCCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTaaaa                                                                                                                                  >  4:596874/1‑119 (MQ=255)
                                              aaCCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTaaaa                                                                                                                                  <  3:596874/119‑1 (MQ=255)
                                                ccTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCctttctt                                                                                                  <  1:699419/149‑1 (MQ=255)
                                                  tGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCtttt                                                                                                >  4:806479/1‑149 (MQ=255)
                                                     tCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGa                                                                                                                 >  4:260433/1‑129 (MQ=255)
                                                     tCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGa                                                                                                                 <  3:260433/129‑1 (MQ=255)
                                                          aaGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgt                                                                                        >  2:465279/1‑149 (MQ=255)
                                                           agagCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgtt                                                                                       <  3:16675/149‑1 (MQ=255)
                                                           agagCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgtt                                                                                       <  1:682902/149‑1 (MQ=255)
                                                           agagCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgtt                                                                                       <  4:305620/149‑1 (MQ=255)
                                                            gagCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgttg                                                                                      <  1:247490/149‑1 (MQ=255)
                                                              gCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgttgtt                                                                                    <  3:667850/149‑1 (MQ=255)
                                                              gCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgttgtt                                                                                    >  4:196828/1‑149 (MQ=255)
                                                              gCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgttgtt                                                                                    <  1:674235/149‑1 (MQ=255)
                                                               cAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGttt                                                                                   <  2:631495/149‑1 (MQ=255)
                                                                aGAAAGCCAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTtta                                                                                  <  4:124276/149‑1 (MQ=255)
                                                                aGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTtta                                                                                  >  4:102868/1‑149 (MQ=255)
                                                                aGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTtta                                                                                  <  4:538334/149‑1 (MQ=255)
                                                                 gAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTttat                                                                                 <  4:849538/149‑1 (MQ=255)
                                                                  aaaGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTttatt                                                                                >  2:850155/1‑149 (MQ=255)
                                                                  aaaGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTntatt                                                                                <  3:861870/149‑1 (MQ=255)
                                                                    aGACAGCGACTGTGTAGCCGGTTGGATCAAGGATATGTTTGGTACTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTaa                                                                              <  2:267359/149‑1 (MQ=255)
                                                                    aGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTaa                                                                              <  4:427742/149‑1 (MQ=255)
                                                                    aGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTaa                                                                              >  3:682111/1‑149 (MQ=255)
                                                                    aGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAACTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTaa                                                                              >  1:504649/1‑149 (MQ=255)
                                                                         gCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCtt                                                                         <  3:861663/149‑1 (MQ=255)
                                                                            aCGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAg                                                                      <  1:850155/149‑1 (MQ=255)
                                                                                tgtAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAAc                                                                  >  4:19401/1‑149 (MQ=255)
                                                                                  tAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGc                                                                <  1:169360/149‑1 (MQ=255)
                                                                                  tAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGc                                                                <  3:144051/149‑1 (MQ=255)
                                                                                    gCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAg                                                                                                                  >  4:810363/1‑97 (MQ=255)
                                                                                    gCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAg                                                                                                                  <  3:810363/97‑1 (MQ=255)
                                                                                    gCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTc                                                              >  3:555615/1‑149 (MQ=255)
                                                                                     ccGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCa                                                             <  2:434440/149‑1 (MQ=255)
                                                                                       ggTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGa                                                           <  3:727086/149‑1 (MQ=255)
                                                                                            gATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCg                                                      <  4:183593/149‑1 (MQ=255)
                                                                                             aTCAATGATATGTTTGGTAGTAATTAGCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGt                                                     >  4:561364/1‑149 (MQ=255)
                                                                                              tCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTa                                                    >  2:429559/1‑149 (MQ=255)
                                                                                              tCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTa                                                    >  4:664509/1‑149 (MQ=255)
                                                                                              tCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTa                                                    >  4:530338/1‑149 (MQ=255)
                                                                                                 aGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCg                                                 >  4:325411/1‑149 (MQ=255)
                                                                                                  ggATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATCCGGTa                                                                                                                                     >  2:615282/1‑64 (MQ=37)
                                                                                                  ggATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATCCGGTa                                                                                                                                     <  1:615282/64‑1 (MQ=37)
                                                                                                   gATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCt                                               >  1:307089/1‑149 (MQ=255)
                                                                                                     tatGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTg                                             <  3:70959/149‑1 (MQ=255)
                                                                                                      atGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGc                                            <  4:667461/149‑1 (MQ=255)
                                                                                                       tGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCTGGTAAAATAACCTCCAATCTAAGATATTAACCATTTCTTTTCATCTGGTTGATTATTAACCCTTCAGTAAAGCTCAGATTGCGTACCGCTTGCg                                           >  2:681532/1‑149 (MQ=255)
                                                                                                         tttGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGaa                                         >  2:392486/1‑149 (MQ=255)
                                                                                                          ttGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAAc                                        <  4:682111/149‑1 (MQ=255)
                                                                                                             gtagtaATTAACATCTAAGCGTGAAATACCGGATGGCGATTTCCCATCCGGTAAAATAACACCCCATCTAAGATATTAATCCTTTCTTTTCATCTGGTTGTTTATTAACCCATCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCg                                     >  2:244254/1‑149 (MQ=255)
                                                                                                              tagtaATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGc                                    >  4:79511/1‑149 (MQ=255)
                                                                                                              tagtaATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGc                                    >  4:604380/1‑149 (MQ=255)
                                                                                                                 taATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAATATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAAAGCTCAGATTGCGTACCGCTTGCGAACCCGCCAg                                 >  2:234906/1‑149 (MQ=255)
                                                                                                                  aaTTAACATCTAAGCGTGAAATACCGGATGGCGAGTTTCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCCTACCGCTTGCGAACCCGCCAGc                                <  2:816603/149‑1 (MQ=255)
                                                                                                                   aTTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCg                               >  1:89872/1‑149 (MQ=255)
                                                                                                                    ttAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGt                              >  4:283258/1‑149 (MQ=255)
                                                                                                                      aaCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGttt                            <  3:772362/149‑1 (MQ=255)
                                                                                                                         aTCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGa                         >  2:806274/1‑149 (MQ=255)
                                                                                                                          tCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGaa                        <  2:335163/149‑1 (MQ=255)
                                                                                                                            tAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAata                      >  3:644355/1‑149 (MQ=255)
                                                                                                                            tAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACACTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAata                      >  1:202651/1‑149 (MQ=255)
                                                                                                                             aaGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAatat                     >  4:336885/1‑149 (MQ=255)
                                                                                                                                 gTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATAttatc                 <  3:79511/149‑1 (MQ=255)
                                                                                                                                    aaaTACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTATCGCTTGCGAACCCGCCAGCGTTTCGAATAttatctta              >  3:10571/1‑149 (MQ=255)
                                                                                                                                      aTACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATAttatcttatc            >  3:840702/1‑149 (MQ=255)
                                                                                                                                         ccGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCttt         <  3:445811/149‑1 (MQ=255)
                                                                                                                                          cGGATGGCAACTCGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTta        >  4:117273/1‑149 (MQ=255)
                                                                                                                                           ggATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTtat       <  2:138402/149‑1 (MQ=255)
                                                                                                                                           ggATGGCAACTCGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTtat       <  3:117273/149‑1 (MQ=255)
                                                                                                                                             aTGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTataa     <  1:392486/149‑1 (MQ=255)
                                                                                                                                                gCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTAtaataa  >  3:797117/1‑149 (MQ=255)
                                                                                                                                                  gAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAAc                                        <  4:23064/109‑1 (MQ=255)
                                                                                                                                                  gAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAAc                                        >  3:23064/1‑109 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAA  >  NC_000913/167145‑167437

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