Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,850,719:1 +T 11.3% intergenic (+26/‑141) sppA → / → ansA protease IV, a signal peptide peptidase/L‑asparaginase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,850,7191.T11.3% 239.8 / 30.5 71intergenic (+26/‑141)sppA/ansAprotease IV, a signal peptide peptidase/L‑asparaginase 1
Reads supporting (aligned to +/- strand):  ref base . (36/27);  new base T (4/4);  total (40/31)
Fisher's exact test for biased strand distribution p-value = 7.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.42e-01

CATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAA  >  NC_000913/1850570‑1850865
                                                                                                                                                      |                                                                                                                                                  
cATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGAtc                                                                                                                                                     >  1:977885/1‑150 (MQ=255)
   gTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGAtcgtc                                                                                                                                                 >  1:480662/1‑150 (MQ=255)
    ttACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGAtcgtcg                                                                                                                                                <  1:1006875/150‑1 (MQ=255)
     tACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCgg                                                                                                                                               <  2:1189310/150‑1 (MQ=255)
      aCCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGc                                                                                                                                              <  1:1028016/150‑1 (MQ=255)
         tGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcac                                                                                                                                           <  1:220493/150‑1 (MQ=255)
           cACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcactc                                                                                                                                         <  1:1234763/150‑3 (MQ=255)
           cACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcactc                                                                                                                                         <  1:113211/150‑3 (MQ=255)
                  ggCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTAttttt                                                                                                                                  >  2:956923/1‑150 (MQ=255)
                    ccTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcac                                                                                                                                           <  2:169091/139‑1 (MQ=255)
                    ccTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACGATCTGTCGGCcac                                                                                                                                           >  1:169091/1‑139 (MQ=255)
                         gTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGTCCACATTATTTTTTACGTcg                                                                                                                           >  1:992278/1‑150 (MQ=255)
                          tAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTcga                                                                                                                          <  2:1246256/150‑1 (MQ=255)
                               tctACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGGAt                                                                                                                     >  2:134172/1‑150 (MQ=255)
                                   ccGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCctc                                                                                                                 <  1:772305/150‑1 (MQ=255)
                                    cGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCctct                                                                                                                <  2:1167216/150‑1 (MQ=255)
                                       tAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTcc                                                                                                             >  2:479317/1‑150 (MQ=255)
                                           agtgaaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCTCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCt                                                                                                         <  2:831696/150‑1 (MQ=255)
                                           agtgaaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCt                                                                                                         <  1:344066/150‑1 (MQ=255)
                                            gtgaaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATc                                                                                                                    <  2:622206/138‑1 (MQ=255)
                                            gtgaaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATc                                                                                                                    >  1:622206/1‑138 (MQ=255)
                                            gtgaaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTg                                                                                                        >  1:40685/1‑150 (MQ=255)
                                               aaagtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGttt                                                                                                     <  2:40685/150‑1 (MQ=255)
                                                  gtgaCAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGc                                                                                                  >  2:958892/1‑150 (MQ=255)
                                                             gCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTcc                                                                                       <  2:1389907/150‑1 (MQ=255)
                                                              ccGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGCCAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCt                                                                                      >  2:967816/1‑150 (MQ=255)
                                                              ccGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCt                                                                                      <  1:1352294/150‑1 (MQ=255)
                                                                gcctttAATGAACCACAAACCCGTTATGCGTTTTGCCTGACCTGAGCAAACATGAGTTAAGTCATCTACAGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCCCCGAATCCTCTTCCCGCTGTGTCGCACCATATTGCCTTa                                                                                    >  2:1012253/4‑150 (MQ=255)
                                                                     tAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGc                                                                               <  2:480662/150‑1 (MQ=255)
                                                                       aTGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCAt                                                                             >  1:1276717/1‑150 (MQ=255)
                                                                        tGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCGTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATa                                                                            >  1:166757/1‑150 (MQ=255)
                                                                          aCCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGt                                                                          >  2:779632/1‑150 (MQ=255)
                                                                           cccACAAATCCGTTATGCGTTTTGCCTGACCTTCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTa                                                                         <  2:219902/150‑1 (MQ=255)
                                                                           cccACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTa                                                                         >  2:123215/1‑150 (MQ=255)
                                                                             cacaAAACTGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTAACGTCGACGAATCCTATTACCGCTGTTTAGACCCATATTTCCTTATACGAATAGTATc                                                                       >  2:933112/1‑150 (MQ=255)
                                                                              acaAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCa                                                                      <  1:123215/150‑1 (MQ=255)
                                                                               aaaaaCCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAg                                                                     >  1:6176/2‑150 (MQ=255)
                                                                                    ccGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgcg                                                                <  1:479317/150‑1 (MQ=255)
                                                                                     cGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTgcgc                                                               <  2:977885/150‑1 (MQ=255)
                                                                                      gTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCt                                                              >  1:610587/1‑150 (MQ=255)
                                                                                             gTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTg                                                       <  2:1276717/150‑1 (MQ=255)
                                                                                                ttGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCt                                                    <  1:99746/150‑1 (MQ=255)
                                                                                                    cTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACg                                                <  1:956923/150‑1 (MQ=255)
                                                                                                      gACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGta                                              >  1:976565/1‑150 (MQ=255)
                                                                                                        ccTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCCCATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGtata                                            >  2:1275588/1‑150 (MQ=255)
                                                                                                        ccTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGtata                                            >  2:154726/1‑150 (MQ=255)
                                                                                                        ccTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTAACAGGTGCGCTCCCCCTGCCTCACGtata                                            >  2:7979/1‑150 (MQ=255)
                                                                                                          tGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTACTTATCCGCATAGTATAAGGTGCGCACCCCCTGCCTCACGtatata                                          >  2:625885/1‑150 (MQ=255)
                                                                                                            cgcCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCATCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTAAATACt                                        >  1:880833/1‑150 (MQ=255)
                                                                                                              ccAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGACGGCCACATTATTTTTTAAGTCGAAGAATCCTCTTCCCGCTGTTTAGCACCAAATTTCCATATCCGCAAAGTATCAGGTGAGCTCCCACTGACTCACGTATATAAttt                                      >  2:659365/1‑150 (MQ=255)
                                                                                                                       gTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTc                             <  1:1107311/150‑1 (MQ=255)
                                                                                                                            gTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGattt                        >  2:942110/1‑148 (MQ=255)
                                                                                                                            gTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGatat                        >  1:617744/1‑150 (MQ=255)
                                                                                                                               ttGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcat                     <  1:154726/150‑1 (MQ=255)
                                                                                                                               ttGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcat                     >  1:800718/1‑150 (MQ=255)
                                                                                                                                 gTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcattc                   >  1:851922/1‑150 (MQ=255)
                                                                                                                                 aaataatgTGGCCGACGATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcattc                   >  1:814065/8‑150 (MQ=255)
                                                                                                                                  taatGAGTCGCCGAAAGATCGTCGGACACATTATTTTTTACTTAGACTAATACACTAAAAAATTTATCGACCCATATTTCCTTATCCGCATAGTATAAGGTGACCTCCCCCTGCCACACGTATATACTTTTGCTCTTTCGATAtcattca                  >  2:201193/4‑150 (MQ=255)
                                                                                                                                  tACTGAGTGGGCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTGTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTTCTCTTTCGATAtcattca                  <  2:367444/150‑1 (MQ=255)
                                                                                                                                  tACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGc                                   <  2:109034/133‑1 (MQ=255)
                                                                                                                                  tACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGc                                   >  1:109034/1‑133 (MQ=255)
                                                                                                                                  tACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcattca                  <  1:625885/150‑1 (MQ=255)
                                                                                                                                   aCTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATAtcattcat                 <  2:1134450/150‑1 (MQ=255)
                                                                                                                                     tGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatat               >  2:276709/1‑150 (MQ=255)
                                                                                                                                      gAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatatc              >  2:444307/1‑150 (MQ=255)
                                                                                                                                      atgTGGCCGACGATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatatc              >  1:843871/3‑150 (MQ=255)
                                                                                                                                        gTGGCCGACGATCTGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatatcaa            >  1:595218/1‑150 (MQ=255)
                                                                                                                                           gCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatatcaatat         >  2:512403/1‑150 (MQ=255)
                                                                                                                                           gCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATACGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCatatcaatat         >  2:55785/1‑150 (MQ=255)
                                                                                                                                              gACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCAt      >  1:645721/1‑150 (MQ=255)
                                                                                                                                                cAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGc    <  1:709335/150‑1 (MQ=255)
                                                                                                                                                aaGATCGTCGGACACATTATTTTTTACGTCGTCGAATGCTCTTCCCGCTGTTTCGGGCCATATTTCCTTATCCGCATAGTATCATGTGCGCTCCCCCTTCCTCACGTATATAATTTTGCTCTTTCGATATCATTCATATCAATATCATGc    <  2:994621/149‑1 (MQ=255)
                                                                                                                                                  gATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCaa  >  2:870513/1‑150 (MQ=255)
                                                                                                                                                  gATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCaa  >  1:327126/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
CATGTTACCTGCACCGCTGGCCTCGGTAGCCTCTACCGTTAAAAGTGAAAGTGACAAGCTGGCCGCGTTTAATGACCCACAAAACCGTTATGCGTTTTGCCTGACCTGCGCCAACATGCGTTAAGTCTTGTACTGAGTGGCCGACAGATCGTCGGCCACATTATTTTTTACGTCGACGAATCCTCTTCCCGCTGTTTCGCCCCATATTTCCTTATCCGCATAGTATCAGGTGCGCTCCCCCTGCCTCACGTATATACTTTTGCTCTTTCGATATCATTCATATCAATATCATGCAA  >  NC_000913/1850570‑1850865

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