Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,351,243 A→G 5.3% Y595H (TAT→CAT)  arcB ← aerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,351,2430AG 5.3% 3010.2 / 94.5 908Y595H (TAT→CAT) arcBaerobic respiration control sensor histidine protein kinase, cognate to two‑component response regulators ArcA and RssB
Reads supporting (aligned to +/- strand):  ref base A (413/446);  new base G (25/24);  total (438/470)
Fisher's exact test for biased strand distribution p-value = 7.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 929 total aligned reads displayed.

TCCCAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATGGCGACATCAACGCTGTTACCTAATTTTTCC  >  NC_000913/3351095‑3351391
                                                                                                                                                    |                                                                                                                                                    
 cccAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATaa                                                                                                                                                     <  3:6932869/149‑1 (MQ=255)
 cccAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTCAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATaa                                                                                                                                                     <  4:5827528/149‑1 (MQ=255)
  ccAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAAc                                                                                                                                                    >  2:4446066/1‑149 (MQ=255)
   cAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATGACg                                                                                                                                                   >  3:1074659/1‑149 (MQ=255)
    aGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGt                                                                                                                                                  >  1:6630459/1‑149 (MQ=255)
      aaTTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGttt                                                                                                                                                >  1:1687559/1‑149 (MQ=255)
      aaTTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGttt                                                                                                                                                <  2:2599916/149‑1 (MQ=255)
       aTTTCTTGATCATCGCTGTAAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTCCTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACGAGCGGCGGTAAATCCTCGCGCGGATAACGGTTc                                                                                                                                               <  2:2059066/149‑1 (MQ=255)
       aTTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTc                                                                                                                                               <  1:3728833/149‑1 (MQ=255)
            ttGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATGACGTTTCGTCAg                                                                                                                                          >  1:523961/1‑149 (MQ=255)
            ttGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAg                                                                                                                                          <  1:6157469/149‑1 (MQ=255)
               atcatcGCGGTTAGCACCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTCCCCGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCCCGCGGATAACGTTTCGTCAGTTc                                                                                                                                       <  4:6519238/149‑1 (MQ=255)
                tcatcGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCa                                                                                                                                      >  1:2102286/1‑149 (MQ=255)
                tcatcGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCa                                                                                                                                      >  4:1173373/1‑149 (MQ=255)
                 catcGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCAc                                                                                                                                     <  3:5517222/149‑1 (MQ=255)
                  atcGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCTCCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACg                                                                                                                                    >  2:4651115/1‑149 (MQ=255)
                  atcGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACg                                                                                                                                    >  1:1197411/1‑149 (MQ=255)
                   tcGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACga                                                                                                                                   >  1:5987513/1‑149 (MQ=255)
                   tcGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACga                                                                                                                                   >  4:6803151/1‑149 (MQ=255)
                     gcgGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACgaga                                                                                                                                 <  3:1359177/149‑1 (MQ=255)
                         ttAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGata                                                                                                                             <  3:1631313/149‑1 (MQ=255)
                             cgcCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCa                                                                                                                         <  1:787147/149‑1 (MQ=255)
                                  gAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCTCGAGAGATCTCCAGCCCg                                                                                                                    <  4:3268275/149‑1 (MQ=255)
                                    aCAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACga                                                                                                                                   <  1:4558950/132‑1 (MQ=255)
                                        aaaGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATGACGTTTCGTCAGTTCActgtag                                                                                                                                >  3:2477920/1‑126 (MQ=255)
                                         aaGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCCTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCatat                                                                                                             <  2:6255463/149‑1 (MQ=255)
                                          aGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATc                                                                                                            <  2:3199148/149‑1 (MQ=255)
                                             ggCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTgg                                                                                                         <  2:6039306/149‑1 (MQ=255)
                                              gCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTc                                                                                                                                               <  4:6248396/110‑1 (MQ=255)
                                                  aCTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACt                                                                                                    >  2:3903175/1‑149 (MQ=255)
                                                   cTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTg                                                                                                   >  2:6855068/1‑149 (MQ=255)
                                                    tCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACga                                                                                                                                   >  3:2170788/1‑116 (MQ=255)
                                                    tCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGa                                                                                                  >  4:799652/1‑149 (MQ=255)
                                                    tCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGa                                                                                                  <  1:4986960/149‑1 (MQ=255)
                                                      aGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAat                                                                                                >  4:1624058/1‑149 (MQ=255)
                                                       gCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATGACGTTTCGTCAGTTCACGAGAGATATCCAg                                                                                                                        >  2:2904654/1‑124 (MQ=255)
                                                       gCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAata                                                                                               >  2:5222591/1‑149 (MQ=255)
                                                            tcatcCATTCCAGCATTGAGGTACTCTTGTTTGTCGTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGATCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAg                                                                                          <  3:4868877/149‑1 (MQ=255)
                                                             catcCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTc                                                                                                                 <  1:3924853/125‑1 (MQ=255)
                                                              atcCATTCCATCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGTCGGTAAATCCTCGCGCGGAAAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATAGCTGGCAACTGAATATCCAGCa                                                                                        >  3:4390942/1‑149 (MQ=255)
                                                               tcCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGcc                                                                                                                      <  1:1923351/118‑1 (MQ=255)
                                                               tcCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCaa                                                                                       >  1:100125/1‑149 (MQ=255)
                                                                ccATTCCAGCATTGAGGTACTCTTGTTTAACTTTCAGGACGTTACCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAac                                                                                      >  2:412543/1‑149 (MQ=255)
                                                                  aTTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAacac                                                                                    <  4:3067478/149‑1 (MQ=255)
                                                                   ttCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACAcc                                                                                   >  1:2617925/1‑149 (MQ=255)
                                                                   ttCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACAcc                                                                                   <  1:3223877/149‑1 (MQ=255)
                                                                     ccAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAg                                                                                 >  4:7000562/1‑149 (MQ=255)
                                                                      cAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATGACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAgg                                                                                >  2:5817928/1‑149 (MQ=255)
                                                                      cAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAgg                                                                                >  3:3961041/1‑149 (MQ=255)
                                                                               ggTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCg                                                                       >  1:7199179/1‑149 (MQ=255)
                                                                                gTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGc                                                                      >  4:4191818/1‑149 (MQ=255)
                                                                                gTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGc                                                                      >  4:1778552/1‑149 (MQ=255)
                                                                                  aatctTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGccc                                                                    <  2:3478324/147‑1 (MQ=255)
                                                                                  aCTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGccc                                                                    <  3:872015/149‑1 (MQ=255)
                                                                                    tctTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACga                                                                                                                                   >  2:1507679/1‑84 (MQ=255)
                                                                                    tctTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGat                                                                                                                              <  4:4564370/89‑1 (MQ=255)
                                                                                    tctTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCTGCAACACCAGGTCGTATTCGCCCgg                                                                  >  3:4667833/1‑149 (MQ=255)
                                                                                    tctTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCgg                                                                  >  4:2034348/1‑149 (MQ=255)
                                                                                    tctTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCgg                                                                  >  4:5505189/1‑149 (MQ=255)
                                                                                    tctTGTTTGTCTTTCAACACGTTAGCTGTTAAGGCCACCAGCGGCGGTAAATGCTCGAGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCCGCCGGGTCGTATCTGGCAACTGAATATCCAGCAACACCCGGTCGTATTCGCCCgg                                                                  >  4:4994831/1‑149 (MQ=255)
                                                                                      ttgtttgtCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGtt                                                                >  3:2806645/1‑149 (MQ=255)
                                                                                         tttgtCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGaa                                                                                                 >  4:4806964/1‑113 (MQ=255)
                                                                                         tttgtCTTTCAGCACGGTAGCGGGTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGTGAGATATCCAGCCCGGTCATAACTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTaa                                                             >  4:4001974/1‑149 (MQ=255)
                                                                                          ttgtCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTaaa                                                            <  3:5016930/149‑1 (MQ=255)
                                                                                            gtCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCa                                                                                                                                           >  4:2884800/1‑68 (MQ=255)
                                                                                             tCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATGACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACAt                                                         >  1:6948565/1‑149 (MQ=255)
                                                                                               tttCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACAtct                                                       >  2:4965981/1‑149 (MQ=255)
                                                                                               tttCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACAtct                                                       >  2:2259061/1‑149 (MQ=255)
                                                                                               tttCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACAtct                                                       >  3:3453888/1‑149 (MQ=255)
                                                                                                ttCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATGACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACAtctc                                                      <  3:324823/149‑1 (MQ=255)
                                                                                                 tCAGCACGGTAGCAGTTAAGGCCACCAGCGCCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTACACGAGAGATATCCAGACCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCGCGGTTTAAACATCTcc                                                     >  4:3654198/1‑149 (MQ=255)
                                                                                                   aGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAg                                                   >  3:3768385/1‑149 (MQ=255)
                                                                                                    gCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAgc                                                  <  1:2434513/149‑1 (MQ=255)
                                                                                                              ggTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGc                                        >  1:4136443/1‑149 (MQ=255)
                                                                                                                   aGGCCACCAGCGGCGGTAAATCCTCGCGCGGATGACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTc                                   <  2:6699633/149‑1 (MQ=255)
                                                                                                                   aGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTc                                   >  4:1034441/1‑149 (MQ=255)
                                                                                                                    gtccaccaGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCa                                  >  4:5349590/3‑149 (MQ=255)
                                                                                                                    ggCCACCAGCGGCGGTAAATCCTCGCGCGGATGACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCa                                  >  1:3016225/1‑149 (MQ=255)
                                                                                                                     gCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCt                                                                                                           >  3:5780334/1‑75 (MQ=255)
                                                                                                                      ccaccaGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATg                                >  3:6211150/1‑149 (MQ=255)
                                                                                                                       caccaGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATgg                               <  1:3568439/149‑1 (MQ=255)
                                                                                                                       caccaGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATgg                               >  3:1873700/1‑149 (MQ=255)
                                                                                                                        accaGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGt                                                                 <  2:1338995/114‑1 (MQ=255)
                                                                                                                          caGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATGGCGa                            <  1:510536/149‑1 (MQ=255)
                                                                                                                           aGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATGGCGAc                           <  3:2962159/149‑1 (MQ=255)
                                                                                                                             cggcggTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATGGCGACAt                         <  1:1422370/149‑1 (MQ=255)
                                                                                                                             cggcggTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATGGCGACAt                         <  1:6330029/149‑1 (MQ=255)
                                                                                                                              ggcggTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATGGCGACATc                        >  4:237856/1‑149 (MQ=255)
                                                                                                                                cggTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATGGCGACATCaa                      <  2:7042897/149‑1 (MQ=255)
                                                                                                                                 ggTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATGGCGACATCAAc                     <  4:269570/149‑1 (MQ=255)
                                                                                                                                 ggTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATGGCGACATCAAc                     <  2:5380623/149‑1 (MQ=255)
                                                                                                                                 ggTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATGGCGACATCAAc                     <  2:1110380/149‑1 (MQ=255)
                                                                                                                                       tCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATGGCGACATCAACGCTGtt               <  3:6803151/149‑1 (MQ=255)
                                                                                                                                       tCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATGGCGACATCAACGCTGtt               <  2:1841483/149‑1 (MQ=255)
                                                                                                                                        ccTCGCGCGGATGACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATGGCGACATCAACGCTGTTa              <  1:563202/149‑1 (MQ=255)
                                                                                                                                             cgcgGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATGGCGACATCAACGCTGTTACCTaa         <  1:5362077/149‑1 (MQ=255)
                                                                                                                                                ggATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATGGCGACATCAACGCTGTTACCTAAttt      >  2:7287977/1‑149 (MQ=255)
                                                                                                                                                 gATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGcc                                                                     <  2:3650717/85‑1 (MQ=255)
                                                                                                                                                   tAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATGGTGACATCAACGCTGTTACCTAATTTTTc   <  1:4459927/149‑1 (MQ=255)
                                                                                                                                                   tAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATGGCGACATCAACGCTGTTACCTAATTTTTc   <  3:6403820/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCCCAGAATTTCTTGATCATCGCGGTTAGCGCCGGAACAGAAAGCGGCTTACTCAGCACATCATCCATTCCAGCATTGAGGTACTCTTGTTTGTCTTTCAGCACGTTAGCGGTTAAGGCCACCAGCGGCGGTAAATCCTCGCGCGGATAACGTTTCGTCAGTTCACGAGAGATATCCAGCCCGGTCATATCTGGCAACTGAATATCCAGCAACACCAGGTCGTATTCGCCCGGTTTAAACATCTCCAGCGCCGCCTTGCCGGTCATGGCGACATCAACGCTGTTACCTAATTTTTCC  >  NC_000913/3351095‑3351391

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