Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 167,293 G→C 5.2% intergenic (+29/‑191) mrcB → / → fhuA peptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB/ferrichrome outer membrane transporter/phage receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913167,2930GC 5.2% 421.1 / 8.9 124intergenic (+29/‑191)mrcB/fhuApeptidoglycan glycosyltransferase/peptidoglycan DD‑transpeptidase MrcB/ferrichrome outer membrane transporter/phage receptor
Reads supporting (aligned to +/- strand):  ref base G (60/56);  new base C (4/2);  total (66/58)
Fisher's exact test for biased strand distribution p-value = 6.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.88e-01

Only 100 of 125 total aligned reads displayed.

CAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCAT  >  NC_000913/167144‑167441
                                                                                                                                                     |                                                                                                                                                    
   tCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTg                                                                                                                                                   <  2:138120/150‑1 (MQ=255)
      gtgtGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCGCAGCAAGAGCAGAAAGACAGCGCCGGTGTAGCCGGTTGGATCAAGGATATGTGTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCa                                                                                                                                                <  1:767443/147‑1 (MQ=255)
       tttGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCAt                                                                                                                                               <  1:731541/150‑1 (MQ=255)
          gATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCg                                                                                                                                            >  2:495077/1‑150 (MQ=255)
           aTCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCgg                                                                                                                                           >  2:1536040/1‑150 (MQ=255)
           aTCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCgg                                                                                                                                           >  1:1283436/1‑150 (MQ=255)
            tCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATCCGGt                                                                                                                                          <  2:1586594/150‑1 (MQ=255)
               gTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTaaa                                                                                                                                       <  1:651865/150‑1 (MQ=255)
                    ctcAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAAc                                                                                                                                  <  2:54780/150‑1 (MQ=255)
                    ctcAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAAc                                                                                                                                  >  1:956742/1‑150 (MQ=255)
                      cAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACAt                                                                                                                                >  2:1520286/1‑150 (MQ=255)
                        gCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATcc                                                                                                                              >  1:270067/1‑150 (MQ=255)
                        gCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGACATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATcc                                                                                                                              <  1:1098577/150‑1 (MQ=255)
                           gcagcaACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTTGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCAt                                                                                                                           >  2:1030204/1‑150 (MQ=255)
                           gcagcaACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCAt                                                                                                                           >  2:573295/1‑150 (MQ=255)
                            cagcaACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATc                                                                                                                          <  2:328499/150‑1 (MQ=255)
                            cagcaACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATc                                                                                                                          <  2:47928/150‑1 (MQ=255)
                              gcaACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTa                                                                                                                        <  2:1018909/150‑1 (MQ=255)
                                aaCAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAg                                                                                                                      <  2:1426175/150‑1 (MQ=255)
                                 aCAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGa                                                                                                                     <  2:375834/150‑1 (MQ=255)
                                  cAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGat                                                                                                                    >  2:617987/1‑150 (MQ=255)
                                   aGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGata                                                                                                                   <  1:410296/150‑1 (MQ=255)
                                    gCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGatat                                                                                                                  >  2:1495912/1‑150 (MQ=255)
                                    gCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGatat                                                                                                                  >  1:1146450/1‑150 (MQ=255)
                                     ccGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATAtt                                                                                                                 >  1:1083162/1‑150 (MQ=255)
                                        cAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAAc                                                                                                              >  2:1510080/1‑150 (MQ=255)
                                         aaCAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAAcc                                                                                                             <  1:671453/150‑1 (MQ=255)
                                          aCAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAAccc                                                                                                            >  1:1223784/1‑150 (MQ=255)
                                           cAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATCCGGt                                                                                                                                          >  2:1567702/1‑119 (MQ=255)
                                           cAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGAAGGCAACTCGCCATCCGGt                                                                                                                                          <  1:1567702/119‑1 (MQ=255)
                                              cAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCctttc                                                                                                        <  2:56060/150‑1 (MQ=255)
                                              cAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCctttc                                                                                                        >  2:1267990/1‑150 (MQ=255)
                                               aaCCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCctttct                                                                                                       >  2:448291/1‑150 (MQ=255)
                                               aaCCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCctttct                                                                                                       <  2:270067/150‑1 (MQ=255)
                                                 ccTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCctttcttt                                                                                                     <  1:1520286/150‑1 (MQ=255)
                                                  cTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCtttt                                                                                                    <  1:1536040/150‑1 (MQ=255)
                                                   tGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTc                                                                                                   >  2:913693/1‑150 (MQ=255)
                                                      tCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATc                                                                                                <  2:263470/150‑1 (MQ=255)
                                                      tCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATc                                                                                                >  1:1054518/1‑150 (MQ=255)
                                                      tCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATc                                                                                                >  2:751383/1‑150 (MQ=255)
                                                      tCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATc                                                                                                <  1:1402968/150‑1 (MQ=255)
                                                         gCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTgg                                                                                             >  1:1157127/1‑150 (MQ=255)
                                                         gCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTgg                                                                                             >  1:1300208/1‑150 (MQ=255)
                                                           aaGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgtt                                                                                           >  1:1286763/1‑150 (MQ=255)
                                                               gCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGttt                                                                                       >  1:341571/1‑150 (MQ=255)
                                                               gCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGttt                                                                                       <  1:230485/150‑1 (MQ=255)
                                                               gCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGttt                                                                                       <  1:224446/150‑1 (MQ=255)
                                                                cAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTtta                                                                                      <  2:1525891/150‑1 (MQ=255)
                                                                cAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTtta                                                                                      <  1:1437457/150‑1 (MQ=255)
                                                                  gAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATTTGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTttatt                                                                                    >  1:529266/1‑150 (MQ=255)
                                                                  gAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATTTGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTttatt                                                                                    >  1:523752/1‑150 (MQ=255)
                                                                  gAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTttatt                                                                                    >  1:1549277/1‑150 (MQ=255)
                                                                   aaaGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTttatta                                                                                   >  2:91508/1‑150 (MQ=255)
                                                                    aaGACAGCGACGTTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTaa                                                                                  >  2:267123/1‑150 (MQ=255)
                                                                    aaGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTaa                                                                                  >  1:802334/1‑150 (MQ=255)
                                                                    aaGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTaa                                                                                  >  2:1335732/1‑150 (MQ=255)
                                                                     aGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAAc                                                                                 <  2:1074964/150‑1 (MQ=255)
                                                                         aGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCtt                                                                             >  2:131214/1‑150 (MQ=255)
                                                                          gCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTc                                                                            <  2:977207/150‑1 (MQ=255)
                                                                          gCGACGGTGTAGCCGGTGGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTc                                                                            <  2:734056/150‑1 (MQ=255)
                                                                           cGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCa                                                                           <  1:951359/150‑1 (MQ=255)
                                                                             aCGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAgg                                                                         <  2:1223784/150‑1 (MQ=255)
                                                                                       cGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGa                                                               >  2:209906/1‑150 (MQ=255)
                                                                                          ttGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTg                                                            <  1:913693/150‑1 (MQ=255)
                                                                                             gATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGt                                                         >  2:1333150/1‑150 (MQ=255)
                                                                                             gATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGt                                                         >  2:1438989/1‑150 (MQ=255)
                                                                                              aTCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTa                                                        >  2:170694/1‑150 (MQ=255)
                                                                                                cAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTAcc                                                      >  2:430364/1‑150 (MQ=255)
                                                                                                 aaGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCg                                                     >  1:916730/1‑150 (MQ=255)
                                                                                                     atatGCTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTg                                                 <  1:1537346/150‑1 (MQ=255)
                                                                                                      tatGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGc                                                <  2:776632/150‑1 (MQ=255)
                                                                                                         gTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGaa                                             >  1:1601921/1‑150 (MQ=255)
                                                                                                          tttGGTAGTAATTTACATCTAATCGTTAAATACCGGATTGCGAGTTGCCATCCGGTAAAATAACATTGCATCTAAGATATTAACCCTTTCTTTTGATCTGGTTGTTTTTTATCCCTTCAGGAACGCTGAGATTGCGTACCGCTTGGGAAc                                            <  2:333261/150‑1 (MQ=255)
                                                                                                            tGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAAccc                                          >  1:1565504/1‑150 (MQ=255)
                                                                                                             ggTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCg                                         >  1:634053/1‑150 (MQ=255)
                                                                                                               tagtaATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGcc                                       >  1:1218733/1‑150 (MQ=255)
                                                                                                               tagtaATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGcc                                       >  1:498549/1‑150 (MQ=255)
                                                                                                               tagtaATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGcc                                       >  1:281098/1‑150 (MQ=255)
                                                                                                                agtaATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCa                                      >  1:1294955/1‑150 (MQ=255)
                                                                                                                     ttAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGtt                                 <  1:751383/150‑1 (MQ=255)
                                                                                                                     ttAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGtt                                 >  2:939877/1‑150 (MQ=255)
                                                                                                                     ttAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGtt                                 <  1:253570/150‑1 (MQ=255)
                                                                                                                      tAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGttt                                >  2:1089144/1‑150 (MQ=255)
                                                                                                                       aaCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTc                               <  2:1470985/150‑1 (MQ=255)
                                                                                                                        aCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCg                              <  2:523752/150‑1 (MQ=255)
                                                                                                                        aCATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTATTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCg                              <  1:1510080/150‑1 (MQ=255)
                                                                                                                          aTCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGaa                            <  1:1030204/150‑1 (MQ=255)
                                                                                                                               aGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAAtatta                       <  2:1218733/150‑1 (MQ=255)
                                                                                                                                  gTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATAttatct                    <  1:131214/150‑1 (MQ=255)
                                                                                                                                    gataTACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATAttatctta                  >  1:1340745/4‑150 (MQ=255)
                                                                                                                                    gAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATAttatctta                  >  2:291090/1‑150 (MQ=255)
                                                                                                                                       aTACCCGATGGCGAATTGCCATCCAGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTATTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATtatcttatct               >  2:1434504/1‑150 (MQ=255)
                                                                                                                                        tACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACTCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTatcttatctt              <  2:1541816/150‑1 (MQ=255)
                                                                                                                                        tACCGGATGGCAACTCGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTatcttatctt              >  2:1538860/1‑150 (MQ=255)
                                                                                                                                          ccGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTta            >  1:701356/1‑150 (MQ=255)
                                                                                                                                          ccGGATGGCAACTCGCCATCCGGTAAAATAACATCCCATCTAAGATATTACCCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTta            >  1:1556854/1‑150 (MQ=255)
                                                                                                                                           cGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAGCCCGCCAGCGTTTCGAATATTATCTTATCTTtat           >  1:500333/1‑150 (MQ=255)
                                                                                                                                               tGGCGAGTTGCCATCCGGTAAAATAACATCCAATCTAAGATATTAACCCTTTCTTTTCATCTGCTTGTTTATTAACCCTTCAGGAACGCTCAGATTGAGTACCGCTTGCGAACCAGCCAGCGTTTCGAATATTATATTATCTTTataata       >  2:795942/1‑150 (MQ=255)
                                                                                                                                                 gCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATaataat     >  2:1060605/1‑150 (MQ=255)
                                                                                                                                                    aGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCAt  <  2:1549277/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCAT  >  NC_000913/167144‑167441

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