Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,233,809 G→A 5.7% S112S (AGC→AGT higA ← antitoxinof the HigB‑HigA toxin‑antitoxin system

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,233,8090GA 5.7% 330.6 / 4.5 103S112S (AGC→AGThigAantitoxinof the HigB‑HigA toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base G (54/43);  new base A (3/3);  total (57/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 105 total aligned reads displayed.

CGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCT  >  NC_000913/3233663‑3233956
                                                                                                                                                  |                                                                                                                                                   
cGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATAGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCg                                                                                                                                                     >  7:47463/1‑147 (MQ=255)
   cTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTc                                                                                                                                                  >  6:46798/1‑147 (MQ=255)
          gggAGATTTAAATCGTTATTTGACGCGCCGGAGGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAAcgcg                                                                                                                                          <  6:42849/147‑1 (MQ=255)
            gagaTTTAAATCGTTATTTGAAGCGCCGTATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTg                                                                                                                                        <  1:207243/148‑1 (MQ=255)
            gagaTTTAAATCGTTATTTGAAGCGCCGGATGAAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTg                                                                                                                                        <  1:187846/148‑1 (MQ=255)
                      tCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGa                                                                                                                              <  2:341928/148‑1 (MQ=255)
                           aTTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTa                                                                                                                          <  5:46798/147‑1 (MQ=255)
                             ttGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTg                                                                                                                       >  7:200620/1‑148 (MQ=255)
                                 aGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCaa                                                                                                                   >  4:80691/1‑148 (MQ=255)
                                     ccGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTc                                                                                                               >  8:196840/1‑148 (MQ=255)
                                     ccGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTc                                                                                                               >  2:125645/1‑148 (MQ=255)
                                       ggATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCg                                                                                                             <  1:125645/148‑1 (MQ=255)
                                         aTGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCgg                                                                                                            >  5:247819/1‑147 (MQ=255)
                                            cAACGCATCCAGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGGTGCCAATTTTTTAGCGTGTTCCAGCTTTAATTTCCTCTTCCCGCTCAAAACGCGTg                                                                                                                                        <  1:309500/116‑1 (MQ=255)
                                             aaCGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTg                                                                                                                       >  8:393872/1‑132 (MQ=255)
                                             aaCGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTg                                                                                                                       <  7:393872/132‑1 (MQ=255)
                                              aCGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTa                                                                                                                          <  3:237282/128‑1 (MQ=255)
                                              aCGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTa                                                                                                                          >  4:237282/1‑128 (MQ=255)
                                              aCGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCgg                                                                                                            >  5:308436/1‑142 (MQ=255)
                                                gCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCgg                                                                                                            >  3:358015/1‑140 (MQ=255)
                                                 cATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGa                                                                                                   <  6:247819/148‑1 (MQ=255)
                                                    ccGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAg                                                                                                >  2:427542/1‑148 (MQ=255)
                                                    ccGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCACTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAg                                                                                                >  6:222857/1‑148 (MQ=255)
                                                     cGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAg                                                                                                <  6:308436/147‑1 (MQ=255)
                                                         aCGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGGTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTa                                                                                           <  1:206712/148‑1 (MQ=255)
                                                            tCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAAc                                                                                        >  3:13750/1‑148 (MQ=255)
                                                            tCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCACTCAAAACGCTTGACACCATAGATTTACTGCCAATTTCCGTCAGATCGGAAAGTGTTaat                                                                                         >  2:163807/1‑146 (MQ=255)
                                                            tCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCACTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAAc                                                                                        >  2:157732/1‑148 (MQ=255)
                                                             cgcgTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACAcc                                                                                                                                   <  3:283779/104‑1 (MQ=255)
                                                             cgcgTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTa                                                                                                                          <  6:266436/113‑1 (MQ=255)
                                                             cgcgTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAAc                                                                                        >  6:275368/1‑147 (MQ=255)
                                                             cgcgTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGGGACAcc                                                                                                                                   >  4:283779/1‑104 (MQ=255)
                                                                 ttAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGaa                                                                                                  <  7:143757/133‑1 (MQ=255)
                                                                 ttAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGaa                                                                                                  >  8:143757/1‑133 (MQ=255)
                                                                  tAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTg                                                                                                                                        >  5:67678/1‑94 (MQ=255)
                                                                  tAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTg                                                                                                                                        <  6:67678/94‑1 (MQ=255)
                                                                   aaTCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATAtt                                                                                  >  1:151314/1‑147 (MQ=255)
                                                                    aTCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGa                                                                                >  2:189543/1‑148 (MQ=255)
                                                                     tCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGAt                                                                               >  4:102239/1‑148 (MQ=255)
                                                                      cAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGAttt                                                                                                                           >  8:254861/1‑103 (MQ=255)
                                                                      cAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGAttt                                                                                                                           <  7:254861/103‑1 (MQ=255)
                                                                        aTAAACAAGGCGGGAGAAATGCCGACTCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGa                                                                                                                                       >  1:378906/1‑89 (MQ=255)
                                                                        aTAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGc                                                                                                           >  4:223131/1‑117 (MQ=255)
                                                                        aTAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGc                                                                                                           <  3:223131/117‑1 (MQ=255)
                                                                        aTAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCa                                                                            >  8:177490/1‑148 (MQ=255)
                                                                         tAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCa                                                                            >  8:333262/1‑147 (MQ=255)
                                                                          aaaCAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCatat                                                                                   <  2:225829/139‑1 (MQ=255)
                                                                          aaaCAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCatat                                                                                   >  1:225829/1‑139 (MQ=255)
                                                                           aaCAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATa                                                                          <  7:177490/147‑1 (MQ=255)
                                                                             cAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAgg                                                                       <  3:80691/148‑1 (MQ=255)
                                                                              aaGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAggg                                                                      >  7:187692/1‑148 (MQ=255)
                                                                               aGGCGGGAGAACTGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCACTCAAAACGCGGGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGt                                                                     <  1:163807/148‑1 (MQ=255)
                                                                               aGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGt                                                                     <  4:358015/148‑1 (MQ=255)
                                                                               aGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCACTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGt                                                                     <  1:157732/148‑1 (MQ=255)
                                                                                ggCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTa                                                                    >  4:326848/1‑148 (MQ=255)
                                                                                ggCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTa                                                                    >  1:201823/1‑148 (MQ=255)
                                                                                         aaTGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATc                                                                              <  1:181532/129‑1 (MQ=255)
                                                                                           tGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGc                                                         <  1:189543/148‑1 (MQ=255)
                                                                                           tGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCACATTGATCCATAAGGGTACGAATCACGGc                                                         >  5:63407/1‑148 (MQ=255)
                                                                                             ccGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCt                                                        <  1:427542/147‑1 (MQ=255)
                                                                                             ccGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCt                                                        <  3:326848/147‑1 (MQ=255)
                                                                                               gAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCtata                                                     >  4:260906/1‑148 (MQ=255)
                                                                                                 aTCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCtata                                                     >  5:15743/1‑146 (MQ=255)
                                                                                                 aTCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCtata                                                     <  6:15743/146‑1 (MQ=255)
                                                                                                     cgTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGc                                                 >  4:180907/1‑146 (MQ=255)
                                                                                                          ccAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATAtt                                                                                  >  5:17078/1‑108 (MQ=255)
                                                                                                          ccAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATAtt                                                                                  <  6:17078/108‑1 (MQ=255)
                                                                                                          ccAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCa                                           <  7:333262/147‑1 (MQ=255)
                                                                                                           cAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCgg                                                                                                    >  6:48742/1‑89 (MQ=255)
                                                                                                           cAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCgg                                                                                                    <  5:48742/89‑1 (MQ=255)
                                                                                                            aaTTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATGGATTTACTGCCAATTTCCGGCAGAACGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATgg                                        >  5:244706/1‑148 (MQ=255)
                                                                                                            aaTTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGc                                            >  3:185102/1‑144 (MQ=255)
                                                                                                            aaTTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGc                                            <  4:185102/144‑1 (MQ=255)
                                                                                                             aTTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTTAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGc                                       >  8:51281/1‑148 (MQ=255)
                                                                                                              ttttttAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAg                                                                        >  6:163759/1‑114 (MQ=255)
                                                                                                              ttttttAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAg                                                                        <  5:163759/114‑1 (MQ=255)
                                                                                                                  ttAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGa                                   <  7:304624/147‑1 (MQ=255)
                                                                                                                  ttAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAg                                  >  4:313881/1‑148 (MQ=255)
                                                                                                                     gCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGcc                                <  5:78000/147‑1 (MQ=255)
                                                                                                                     gCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCa                               <  3:260906/148‑1 (MQ=255)
                                                                                                                      cGTGTTCCAGCGTTAATTTCCTCTTCCCGCTTAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCAt                              <  7:51281/148‑1 (MQ=255)
                                                                                                                       gtgtTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATg                             <  8:200620/148‑1 (MQ=255)
                                                                                                                       gtgtTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATg                             >  1:425755/1‑148 (MQ=255)
                                                                                                                         gtTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATgg                            >  3:379779/1‑147 (MQ=255)
                                                                                                                           tCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCAt                         >  4:302290/1‑148 (MQ=255)
                                                                                                                             cAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCAtt                        >  5:426426/1‑147 (MQ=255)
                                                                                                                             cAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCAtt                        >  2:168643/1‑147 (MQ=255)
                                                                                                                               gCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGc                                       >  3:253225/1‑130 (MQ=255)
                                                                                                                               gCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGc                                       <  4:253225/130‑1 (MQ=255)
                                                                                                                                 gTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATGGATTTACTGCCAATTTCCGGCAGAACGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAAtt                   <  6:244706/148‑1 (MQ=255)
                                                                                                                                  ttAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTc                  >  8:344467/1‑148 (MQ=255)
                                                                                                                                       ttCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTaaa                                                                                         >  1:201568/1‑72 (MQ=255)
                                                                                                                                       ttCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTaaa                                                                                         <  2:201568/72‑1 (MQ=255)
                                                                                                                                       ttCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCa              >  8:192842/1‑147 (MQ=255)
                                                                                                                                          ctctTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAAtt          >  3:164461/1‑148 (MQ=255)
                                                                                                                                           tctTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTc         >  3:192905/1‑148 (MQ=255)
                                                                                                                                               cccACTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCggg      <  5:222857/147‑1 (MQ=255)
                                                                                                                                                ccGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGcc                                >  3:318006/1‑120 (MQ=255)
                                                                                                                                                ccGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGcc                                <  4:318006/120‑1 (MQ=255)
                                                                                                                                                  gCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCt  >  5:105621/1‑148 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGGCTGGCGCTGGAGATTTAAATCGTTATTTGAAGCGCCGGATGCAACGCATCCGGCACGTCGCGTTAATCAATAAACAAGGCGGGAGAAATGCCGAATCGCGTTGCCAATTTTTTAGCGTGTTCCAGCGTTAATTTCCTCTTCCCGCTCAAAACGCGTGACACCATAGATTTACTGCCAATTTCCGGCAGATCGGAAAGGGTTAAACCATATTGATCCATAAGGGTACGAATCACGGCTATACCGCCAGGCATGGCTTGAGCCATGGCATTAAATTCCGCAAATTCGGGCGCT  >  NC_000913/3233663‑3233956

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