Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 654,563 C→G 9.5% intergenic (+21/+20) citB → / ← dcuC response regulator in two‑component regulatory system with CitA/anaerobic C4‑dicarboxylate transport

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913654,5630CG 9.5% 326.5 / 21.0 95intergenic (+21/+20)citB/dcuCresponse regulator in two‑component regulatory system with CitA/anaerobic C4‑dicarboxylate transport
Reads supporting (aligned to +/- strand):  ref base C (41/45);  new base G (3/6);  total (44/51)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

AGGCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCAT  >  NC_000913/654415‑654711
                                                                                                                                                    |                                                                                                                                                    
aGGCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAAc                                                                                                                                                      <  1:616513/149‑1 (MQ=255)
  gCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACtt                                                                                                                                                    <  4:1089408/149‑1 (MQ=255)
    aCTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGc                                                                                                                                                  <  1:927489/149‑1 (MQ=255)
       gACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGTAAACTTGCCTc                                                                                                                                               <  1:1084798/149‑1 (MQ=255)
           aTCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAATGCCGGAGGCAAGTTTCCTCCGGc                                                                                                                                           <  3:1295320/149‑1 (MQ=255)
            tCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCa                                                                                                                                          >  2:726629/1‑149 (MQ=255)
             cAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCAt                                                                                                                                         <  1:722104/149‑1 (MQ=255)
                ccGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTAc                                                                                                                                      >  4:1299120/1‑149 (MQ=255)
                  gCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTa                                                                                                                                    >  1:1425126/1‑149 (MQ=255)
                  gCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTa                                                                                                                                    >  2:535195/1‑149 (MQ=255)
                    accacTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAtt                                                                                                                                  <  4:1248140/149‑1 (MQ=255)
                      cacTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACtt                                                                                                                                                    >  2:956626/1‑129 (MQ=255)
                      cacTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACtt                                                                                                                                                    <  1:956626/129‑1 (MQ=255)
                               gCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCtgt                                                                                                                       >  3:124603/1‑149 (MQ=255)
                                 gTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGa                                                                                                                     <  4:124603/149‑1 (MQ=255)
                                     tCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCgg                                                                                                                                            >  4:138607/1‑122 (MQ=255)
                                     tCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCgg                                                                                                                                            <  3:138607/122‑1 (MQ=255)
                                     tCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgc                                                                                                                 <  2:1425126/149‑1 (MQ=255)
                                       ttGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGtt                                                                                                                                                    <  2:422737/112‑1 (MQ=255)
                                       ttGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGtt                                                                                                                                                    <  1:317500/112‑1 (MQ=255)
                                       ttGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCa                                                                                                                                          <  3:1096172/122‑1 (MQ=255)
                                       ttGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCa                                                                                                                                          >  4:1096172/1‑122 (MQ=255)
                                       ttGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctg                                                                                                               >  2:168648/1‑149 (MQ=255)
                                       ttGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctg                                                                                                               >  2:891294/1‑149 (MQ=255)
                                        tGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctgc                                                                                                              >  2:935924/1‑149 (MQ=255)
                                         gAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACtt                                                                                                                                                    <  3:1259355/110‑1 (MQ=255)
                                         gAATATTGCACCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACtt                                                                                                                                                    >  4:1259355/1‑110 (MQ=255)
                                            tatTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctgc                                                                                                              >  3:1415069/1‑145 (MQ=255)
                                            tatTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctgc                                                                                                              <  4:1415069/145‑1 (MQ=255)
                                              ttGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACtt                                                                                                                                                    <  4:373265/105‑1 (MQ=255)
                                              ttGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACtt                                                                                                                                                    >  3:373265/1‑105 (MQ=255)
                                                  gcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCt                                                                                                    <  3:1299120/149‑1 (MQ=255)
                                                   ccAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTg                                                                                                   <  1:935924/149‑1 (MQ=255)
                                                    cAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACATTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTgg                                                                                                  <  1:891294/149‑1 (MQ=255)
                                                         gccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCa                                                                                             <  1:535195/149‑1 (MQ=255)
                                                                gATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACtt                                                                                                                                                    >  2:51429/1‑87 (MQ=255)
                                                                gATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGAGAAAGCCGGAGGAAACtt                                                                                                                                                    <  1:51429/87‑1 (MQ=255)
                                                                        cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAAc                                                                                                                                                      <  1:892512/77‑1 (MQ=255)
                                                                        cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAAc                                                                                                                                                      >  2:892512/1‑77 (MQ=255)
                                                                        cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAtt                                                                                                                                  >  2:152977/1‑97 (MQ=255)
                                                                        cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAtt                                                                                                                                  <  1:152977/97‑1 (MQ=255)
                                                                        cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCAcc                                                                                              >  4:488723/1‑133 (MQ=255)
                                                                        cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCAcc                                                                                              <  3:488723/133‑1 (MQ=255)
                                                                        cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCaa                                                                              <  1:726629/149‑1 (MQ=255)
                                                                        cTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCaa                                                                              <  1:168648/149‑1 (MQ=255)
                                                                         tGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAAc                                                                             >  4:38863/1‑149 (MQ=255)
                                                                           aaaTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTa                                                                                                                                       >  1:529734/1‑89 (MQ=255)
                                                                           aaaTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTa                                                                                                                                       <  2:529734/89‑1 (MQ=255)
                                                                            aaTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCt                                                                                                    <  4:443875/123‑1 (MQ=255)
                                                                            aaTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCt                                                                                                    >  3:443875/1‑123 (MQ=255)
                                                                              ttGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGcc                                                                                                          <  4:137670/115‑1 (MQ=255)
                                                                              ttGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGcc                                                                                                          >  3:137670/1‑115 (MQ=255)
                                                                                   cACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTAc                                                                                                                               >  2:1347091/1‑89 (MQ=255)
                                                                                   cACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTAc                                                                                                                               <  1:1347091/89‑1 (MQ=255)
                                                                                    aCGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTgg                                                                                                  <  4:735885/117‑1 (MQ=255)
                                                                                    aCGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTgg                                                                                                  >  3:735885/1‑117 (MQ=255)
                                                                                    aCGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCa                                                                  <  3:38863/149‑1 (MQ=255)
                                                                                       gCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCTGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGa                                                                           >  4:1220641/1‑137 (MQ=255)
                                                                                       gCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCTGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGa                                                                           <  3:1220641/137‑1 (MQ=255)
                                                                                          aaGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCgg                                                                                                                                            >  3:716543/1‑69 (MQ=255)
                                                                                          aaGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCgg                                                                                                                                            <  4:716543/69‑1 (MQ=255)
                                                                                             ttGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAt                                                                                                                                   <  1:446638/75‑1 (MQ=255)
                                                                                             ttGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAt                                                                                                                                   >  2:446638/1‑75 (MQ=255)
                                                                                             ttGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaaccaa                                                         >  1:46182/1‑149 (MQ=255)
                                                                                               ggCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACcaaca                                                       >  4:409293/1‑149 (MQ=255)
                                                                                               ggCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACcaaca                                                       >  2:675320/1‑149 (MQ=255)
                                                                                                   gACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAAc                                                                             >  1:481521/1‑123 (MQ=255)
                                                                                                   gACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAAc                                                                             <  2:481521/123‑1 (MQ=255)
                                                                                                    aCCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCAc                                                  <  1:675320/149‑1 (MQ=255)
                                                                                                      cacaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTgg                                                                                                  >  3:1425260/1‑99 (MQ=255)
                                                                                                      cacaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTgg                                                                                                  <  4:1425260/99‑1 (MQ=255)
                                                                                                      cacaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCa                                                                  <  1:285923/131‑1 (MQ=255)
                                                                                                      cacaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCa                                                                  >  2:285923/1‑131 (MQ=255)
                                                                                                            gCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACtt                                                                                                                             <  2:810444/66‑1 (MQ=255)
                                                                                                            gCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACtt                                                                                                                             >  1:810444/1‑66 (MQ=255)
                                                                                                              atatACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTg                                                                                    <  4:564455/105‑1 (MQ=255)
                                                                                                              atatACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTg                                                                                    >  3:564455/1‑105 (MQ=255)
                                                                                                                   ccACAGTGGGTGACATGATAAAGCCGGAGTAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTcc                                                                                                     >  1:1147818/1‑83 (MQ=255)
                                                                                                                   ccACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTcc                                                                                                     <  2:1147818/83‑1 (MQ=255)
                                                                                                                         tGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTAc                             <  2:46182/149‑1 (MQ=255)
                                                                                                                                cATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGt                                      <  2:967873/133‑1 (MQ=255)
                                                                                                                                cATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGt                                      >  1:967873/1‑133 (MQ=255)
                                                                                                                                 aTGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctct                                                                                     >  2:32610/1‑85 (MQ=255)
                                                                                                                                 aTGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctct                                                                                     <  1:32610/85‑1 (MQ=255)
                                                                                                                                     tAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACgg                 <  1:1366529/149‑1 (MQ=255)
                                                                                                                                      aaaGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTa                                                                                  <  3:1305536/83‑1 (MQ=255)
                                                                                                                                      aaaGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTa                                                                                  >  4:1305536/1‑83 (MQ=255)
                                                                                                                                        aGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCga              <  3:409293/149‑1 (MQ=255)
                                                                                                                                         gCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCa                                                                                          >  2:785716/1‑72 (MQ=38)
                                                                                                                                         gCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCa                                                                                          <  1:785716/72‑1 (MQ=38)
                                                                                                                                         gCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAAc                                                                       >  4:888955/1‑91 (MQ=255)
                                                                                                                                         gCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAAc                                                                       <  3:888955/91‑1 (MQ=255)
                                                                                                                                          ccGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTa                                             <  1:903846/116‑1 (MQ=255)
                                                                                                                                          ccGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTa                                             >  2:903846/1‑116 (MQ=255)
                                                                                                                                              aggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCtt        >  3:308729/1‑149 (MQ=255)
                                                                                                                                                gaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCg      >  2:1194208/1‑149 (MQ=255)
                                                                                                                                                 aaaCTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctc                                                                                      >  2:1424545/1‑68 (MQ=255)
                                                                                                                                                 aaaCTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctc                                                                                      <  1:1424545/68‑1 (MQ=255)
                                                                                                                                                    cTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCAt  >  4:291088/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
AGGCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCTTCGCCAT  >  NC_000913/654415‑654711

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