Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,850,593 C→T 6.2% G45E (GGG→GAG)  uhpA ← response regulator in two‑component regulatory system wtih UhpB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,850,5930CT 6.2% 316.6 / 7.2 96G45E (GGG→GAG) uhpAresponse regulator in two‑component regulatory system wtih UhpB
Reads supporting (aligned to +/- strand):  ref base C (50/40);  new base T (3/3);  total (53/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.91e-01

CTCAACCAGCGCAGGACTGTCGTGAACGGAGAGCATAATCGTCGCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCC  >  NC_000913/3850463‑3850736
                                                                                                                                  |                                                                                                                                               
cTCAACCAGCGCAGGACTGTCGTGAACGGAGAGCATAATCGTCGCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAgc                                                                                                                                >  5:313867/1‑148 (MQ=255)
      cAGCGCAGGACTGTCGTGAACGGAGAGCATAATCGTCGCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACCCCTGCACACCGCGCCCCg                                                                                                                                              >  3:128483/1‑128 (MQ=255)
      cAGCGCAGGACTGTCGTGAACGGAGAGCATAATCGTCGCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCg                                                                                                                                              <  4:128483/128‑1 (MQ=255)
        gcgcAGGACTGTCGTGAACGGAGAGCATAATCGTCGCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTcgcg                                                                                                                        >  2:321263/1‑148 (MQ=255)
          gcAGGACTGTCGTGAACGGAGAGCATAATCGTCGCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGcc                                                                                                                                        <  4:298771/130‑1 (MQ=255)
          gcAGGACTGTCGTGAACGGAGAGCATAATCGTCGCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGcc                                                                                                                                        >  3:298771/1‑130 (MQ=255)
               aCTGTCGTGAACGGAGAGCATAATCGTCGCCATACCTTTCGGCGGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCgca                                                                                                                 >  3:221899/1‑146 (MQ=255)
                  gTCGTGAACGGAGAGCATAATCGTCGCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAAcc                                                                                                               >  5:115799/1‑147 (MQ=255)
                       gAACGGAGAGCATAATAGTCGCCATACCTTACGGCAGCTGGAATAGAAGCTCAAGACCGGATATATCGGGCATGGACATATCGCAAATACACACCTCCACACCGCCCCCCGACAGCCCCGCCAGCGCCTCTCGCCCCGAACCAACCTc                                                                                                         >  2:266880/1‑148 (MQ=255)
                           ggAGAGCATAATCGTCGCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCa                                                                                                     <  1:266880/148‑1 (MQ=255)
                           ggAGAGCATAATCGTCGCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCa                                                                                                     <  1:321263/148‑1 (MQ=255)
                               agCATAATCGTCGCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTa                                                                                                 <  6:313867/148‑1 (MQ=255)
                                   tAATCGTCGCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTg                                                                                             <  1:327160/148‑1 (MQ=255)
                                   tAATCGTCGCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTg                                                                                             <  7:212868/148‑1 (MQ=255)
                                    aaTCGTCGCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAAc                                                                                                           >  1:163074/1‑133 (MQ=255)
                                    aaTCGTCGCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAAc                                                                                                           <  2:163074/133‑1 (MQ=255)
                                    aaTCGTCGCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTg                                                                                             >  7:179767/1‑147 (MQ=255)
                                       cgtcgCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTc                                                                                                         >  3:248066/1‑132 (MQ=255)
                                       cgtcgCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTc                                                                                                         <  4:248066/132‑1 (MQ=255)
                                                aCCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACt                                                                                                          <  2:292435/122‑1 (MQ=255)
                                                aCCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACt                                                                                                          >  1:292435/1‑122 (MQ=255)
                                                aCCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCa                                                                                                     <  4:307687/127‑1 (MQ=255)
                                                aCCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCa                                                                                                     >  3:307687/1‑127 (MQ=255)
                                                aCCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTc                                                                                 >  7:266001/1‑147 (MQ=255)
                                                 ccTTTCGGCAGCTGGCTTAGCAGCTCCAGCCCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCa                                                                               >  5:299087/1‑148 (MQ=255)
                                                 ccTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTcc                                                                                <  8:179767/147‑1 (MQ=255)
                                                        gCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCcag                                                                        >  8:136442/1‑148 (MQ=255)
                                                         cAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGcccc                                                                          >  6:338245/1‑145 (MQ=255)
                                                         cAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGcccc                                                                          <  5:338245/145‑1 (MQ=255)
                                                         cAGCTGGCTTAGCAGCTCCAGACCGGAGAAATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTcc                                                                                <  8:149503/139‑1 (MQ=255)
                                                         cAGCTGGCTTAGCAGCTCCAGACCGGAGAAATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTcc                                                                                >  7:149503/1‑139 (MQ=255)
                                                                     cagcTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGc                                                                                            >  7:323200/1‑115 (MQ=255)
                                                                     cagcTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGc                                                                                            <  8:323200/115‑1 (MQ=255)
                                                                         tCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCg                                                        >  2:217249/1‑147 (MQ=255)
                                                                                  gagaTATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCg                                                                                                                                     <  7:9614/61‑1 (MQ=255)
                                                                                  gagaTATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCg                                                                                                                                     >  8:9614/1‑61 (MQ=255)
                                                                                   agaTATCGGGCAGGGAGATATCGCAAATACACCCCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGCACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGAt                                             >  5:342313/1‑148 (MQ=255)
                                                                                     atatCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGAt                                             >  3:32343/1‑146 (MQ=255)
                                                                                      tatCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCa                                                                                                     >  3:41413/1‑89 (MQ=255)
                                                                                      tatCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCa                                                                                                     <  4:41413/89‑1 (MQ=255)
                                                                                      tatCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAg                                          >  5:371689/1‑148 (MQ=255)
                                                                                         cGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAg                                                     >  6:166138/1‑134 (MQ=255)
                                                                                         cGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAg                                                     <  5:166138/134‑1 (MQ=255)
                                                                                          gggCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGa                                      <  6:299087/148‑1 (MQ=255)
                                                                                           ggCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGAt                                     <  6:115799/148‑1 (MQ=255)
                                                                                                ggggaTATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGAt                                             >  4:15419/4‑135 (MQ=255)
                                                                                                ggggaTATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGAt                                             <  3:15419/132‑1 (MQ=255)
                                                                                                ggAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGcccc                                                                          >  8:7228/1‑106 (MQ=255)
                                                                                                ggAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGcccc                                                                          <  7:7228/106‑1 (MQ=255)
                                                                                                ggAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGc                                                          >  3:146377/1‑122 (MQ=255)
                                                                                                ggAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGc                                                          <  4:146377/122‑1 (MQ=255)
                                                                                                ggAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACg                                               <  5:250478/133‑1 (MQ=255)
                                                                                                ggAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACg                                               >  6:250478/1‑133 (MQ=255)
                                                                                                ggAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTc                                 >  4:217338/1‑147 (MQ=255)
                                                                                                ggAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTc                                 >  8:325394/1‑147 (MQ=255)
                                                                                                ggAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTc                                 >  2:11897/1‑147 (MQ=255)
                                                                                                 gagaTATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGa                                                      <  6:185448/125‑1 (MQ=255)
                                                                                                 gagaTATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCt                                <  8:266001/147‑1 (MQ=255)
                                                                                                 gagaTATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAAACAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGa                                                      >  5:185448/1‑125 (MQ=255)
                                                                                                   gaTATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCtata                             <  1:11897/148‑1 (MQ=255)
                                                                                                     tatCGCAAATACACACCTGCACACCTCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCCAACTCGGCAACTCCCTGCCAATCAGGTTCCCGCCCCCGCAGCGGCGCAACGCCGGCGCGGACGAGGCGGTGCGCGTCtatacg                           >  3:34816/1‑146 (MQ=255)
                                                                                                     tatCGCAAATACACACCTGCACACCGCGCTCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAACCTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCtat                              <  3:328929/145‑1 (MQ=255)
                                                                                                     tatCGCAAATACACACCTGCACACCGCGCTCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGTTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCtat                              <  1:305407/145‑1 (MQ=255)
                                                                                                     tatCGCAAATACACACCTGCACACCGCGCTCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCtat                              >  4:328929/1‑145 (MQ=255)
                                                                                                     tatCGCAAATACACACCTGCACACCGCGCTCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCtat                              <  7:21831/145‑1 (MQ=255)
                                                                                                     tatCGCAAATACACACCTGCACACCGCGCTCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCtat                              >  8:21831/1‑145 (MQ=255)
                                                                                                     tatCGCAAATACACACCTGCACACCGCGCTCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCtat                              >  2:305407/1‑145 (MQ=255)
                                                                                                        cGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAggg                         >  7:75768/1‑147 (MQ=255)
                                                                                                        cGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAggg                         >  8:154146/1‑147 (MQ=255)
                                                                                                        cGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAggg                         >  5:136746/1‑147 (MQ=255)
                                                                                                        cGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAggg                         >  5:320782/1‑147 (MQ=255)
                                                                                                           aaaTACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAAc                     >  7:390317/1‑148 (MQ=255)
                                                                                                           aaaTACACACCTGCAAACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGGCGATGAGGTGATCGTCTATAAGGGCAAc                     >  7:288776/1‑148 (MQ=255)
                                                                                                             aTACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCt                                <  2:372444/135‑1 (MQ=255)
                                                                                                             aTACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCt                                >  1:372444/1‑135 (MQ=255)
                                                                                                             aTACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACgg                   >  1:272397/1‑148 (MQ=255)
                                                                                                              tACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCtata                             <  4:260370/137‑1 (MQ=255)
                                                                                                              tACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCtata                             >  3:260370/1‑137 (MQ=255)
                                                                                                                cacacCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGACCCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGCAGCGGACGATGAGGTGATCGTCGATAAGGGCCACGGTGa                >  2:263427/1‑148 (MQ=255)
                                                                                                                cacacCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGa                <  8:75768/148‑1 (MQ=255)
                                                                                                                 acacCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCt                                >  6:364219/1‑131 (MQ=255)
                                                                                                                 acacCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCt                                <  5:364219/131‑1 (MQ=255)
                                                                                                                  cacCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATc              <  7:325394/148‑1 (MQ=255)
                                                                                                                  cacCTGCACACCGCGCCCCGGCAGACCAGCCAGCCACTCTCTCACCGAACCAAACTCGGCAAATACCTGCACATCCGGTTACAGCCCCAGCCGCTGCGCAACGCAGGAGCGGACGATGAGGGGATCGGCTATAAGGGCAACGGTGATc              <  4:34816/148‑1 (MQ=255)
                                                                                                                     cTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCACACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATg           >  1:225737/1‑148 (MQ=255)
                                                                                                                      tGCACACCGCTCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCATAGCCGGAGCGGACg                                               <  5:52910/111‑1 (MQ=255)
                                                                                                                      tGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATgg          >  4:333770/1‑148 (MQ=255)
                                                                                                                       gCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGt         >  4:54175/1‑148 (MQ=255)
                                                                                                                         acacCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCt       <  4:32343/148‑1 (MQ=255)
                                                                                                                         acacCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCt       <  6:136746/148‑1 (MQ=255)
                                                                                                                           acCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTg     <  6:342313/148‑1 (MQ=255)
                                                                                                                           acCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCACCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTg     >  8:133008/1‑148 (MQ=255)
                                                                                                                            ccGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGTTCGTCTATAAGGGCAACGGTGATCATGGTCTTGt    >  1:99523/1‑148 (MQ=255)
                                                                                                                              gcgcCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGGCGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTcc  <  8:288776/148‑1 (MQ=255)
                                                                                                                              gcgcCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTcc  <  8:390317/148‑1 (MQ=255)
                                                                                                                              gcgcCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGGGTTCGTCTATAAGGGCAACGGTGATCATGGTCTTGTcc  <  2:99523/148‑1 (MQ=255)
                                                                                                                                  cccGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAg                           <  6:15351/119‑1 (MQ=255)
                                                                                                                                  cccGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAg                           >  5:15351/1‑119 (MQ=255)
                                                                                                                                  |                                                                                                                                               
CTCAACCAGCGCAGGACTGTCGTGAACGGAGAGCATAATCGTCGCCATACCTTTCGGCAGCTGGCTTAGCAGCTCCAGACCGGAGATATCGGGCATGGAGATATCGCAAATACACACCTGCACACCGCGCCCCGGCAGCCCCGCCAGCGCCTCGCGCCCCGAACCAAACTCGGCAACTACCTGCAAATCAGGTTCCAGCCCCAGCAGCTGCGCAAAGCCGGAGCGGACGATGAGGTGATCGTCTATAAGGGCAACGGTGATCATGGTCTTGTCC  >  NC_000913/3850463‑3850736

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: