Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 654,560 A→C 7.0% intergenic (+18/+23) 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,5600AC 7.0% 354.1 / 15.1 114intergenic (+18/+23)citB/dcuCresponse regulator in two‑component regulatory system with CitA/anaerobic C4‑dicarboxylate transport
Reads supporting (aligned to +/- strand):  ref base A (51/55);  new base C (4/4);  total (55/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.53e-01

Only 100 of 118 total aligned reads displayed.

CAGGCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCT  >  NC_000913/654414‑654704
                                                                                                                                                  |                                                                                                                                                
cAGGCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGaaa                                                                                                                                                <  4:90920/149‑1 (MQ=255)
 aGGCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAAc                                                                                                                                               <  4:63821/149‑1 (MQ=255)
 aGGCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAAc                                                                                                                                               <  2:183030/149‑1 (MQ=255)
  ggCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACtt                                                                                                                                             <  5:15158/149‑1 (MQ=255)
   gCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACtt                                                                                                                                             <  3:408866/149‑1 (MQ=255)
     aCTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGc                                                                                                                                           >  7:213021/1‑149 (MQ=255)
           cATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCgg                                                                                                                                     <  1:173825/149‑1 (MQ=255)
            aTCAGCCGCACCCCTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCgg                                                                                                                                     >  5:101181/1‑148 (MQ=255)
             tCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCa                                                                                                                                   >  7:125657/1‑149 (MQ=255)
                gCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGc                                                                                                                                    <  3:359814/145‑1 (MQ=255)
                gCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTa                                                                                                                                <  7:451183/149‑1 (MQ=255)
                gCCGAACCACTGCCAGGCGATATCTTGAATATTTCGCCACCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAATTGGCAGACCACAACGCATATACAACAGTGGGTGACATGATAAAGCCTGAGGAAACTTTCCTCCCGc                                                                                                                                    >  4:359814/1‑145 (MQ=255)
                   gCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTa                                                                                                                             <  2:423778/149‑1 (MQ=255)
                     agctcTGGCATGCGTTATCTTGAATCTTGCGCCATCAGACATCTGATTATTTCTGATATTTTTAACGGCATAGTTGGCAGACCACATCGCAAATTGCACTGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAtt                                                                                                                           <  4:423978/145‑1 (MQ=255)
                     accacTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAtt                                                                                                                           <  7:258947/149‑1 (MQ=255)
                      ccacTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCAGGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAttt                                                                                                                          <  6:205591/149‑1 (MQ=255)
                      ccacTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAtt                                                                                                                           >  7:302992/1‑148 (MQ=255)
                        acTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGATCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTAc                                                                                                                        <  1:459091/149‑1 (MQ=255)
                          tGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGTTTCCTCCGGct                                                                                                                                   <  3:21884/136‑2 (MQ=255)
                            ccAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGc                                                                                                                    >  6:242003/1‑149 (MQ=255)
                              aGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGc                                                                                                                                           >  5:360438/1‑124 (MQ=255)
                              aGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGc                                                                                                                                           <  6:360438/124‑1 (MQ=255)
                              aGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCt                                                                                                                  <  7:279917/149‑1 (MQ=255)
                              aGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCt                                                                                                                  <  7:435000/149‑1 (MQ=255)
                              aGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCt                                                                                                                  <  7:15902/149‑1 (MQ=255)
                                gCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCtgt                                                                                                                <  2:94405/149‑1 (MQ=255)
                                        ttGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgctg                                                                                                        >  2:110955/1‑149 (MQ=255)
                                          gAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCgc                                                                                                          <  6:204624/145‑1 (MQ=255)
                                             tatTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGcc                                                                                                   >  5:162832/1‑149 (MQ=255)
                                               ttGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCg                                                                                                  <  7:75057/148‑1 (MQ=255)
                                                 gcgcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGtt                                                                                                <  6:162832/148‑1 (MQ=255)
                                                   gcCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCt                                                                                             <  1:110955/149‑1 (MQ=255)
                                                    ccAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACTCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTg                                                                                            <  6:114652/149‑1 (MQ=255)
                                                       gccgccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCa                                                                                         <  6:69572/149‑1 (MQ=255)
                                                        ccgccATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCAc                                                                                        >  5:295478/1‑149 (MQ=255)
                                                             aTCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCa                                                                                   <  8:301085/149‑1 (MQ=255)
                                                                  attattGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctct                                                                              >  3:418026/1‑149 (MQ=255)
                                                                  attattGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctc                                                                               <  8:213021/148‑1 (MQ=255)
                                                                  attattGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTAGTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctct                                                                              >  2:127665/1‑149 (MQ=255)
                                                                  attattGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACAACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGctct                                                                              >  1:109120/1‑149 (MQ=255)
                                                                     attGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGt                                                                            <  2:366542/148‑1 (MQ=255)
                                                                      ttGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAg                                                                          >  5:393793/1‑149 (MQ=255)
                                                                          tGAAATTGTTCACGGCAAAGTTGGCAGACTACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAAc                                                                      >  4:328003/1‑149 (MQ=255)
                                                                          tGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAAc                                                                      >  5:179869/1‑149 (MQ=255)
                                                                          tGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAAc                                                                      >  1:247625/1‑149 (MQ=255)
                                                                          tGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAAc                                                                      >  4:319234/1‑149 (MQ=255)
                                                                           gAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGAGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAAc                                                                      >  5:340892/1‑148 (MQ=255)
                                                                             aaTTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGAGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGATGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGAt                                                                   <  6:340892/149‑1 (MQ=255)
                                                                                  ttCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGc                                                                                                                                           >  3:260967/1‑72 (MQ=255)
                                                                                  ttCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGc                                                                                                                                           <  4:260967/72‑1 (MQ=255)
                                                                                  ttCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACaa                                                              >  4:28904/1‑149 (MQ=255)
                                                                                   tCACGGCAAAGTTGGCAGCCCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAACCTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAAt                                                             >  1:308713/1‑149 (MQ=255)
                                                                                   tCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAAt                                                             >  6:287499/1‑149 (MQ=255)
                                                                                     aCGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGc                                                                                                                    <  6:108901/92‑1 (MQ=255)
                                                                                     aCGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATc                                                            <  4:418026/148‑1 (MQ=255)
                                                                                        gCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTGCCTGGCACCATCAGCTCTGTAGCAACGATAACCCTCacc                                                         >  3:287222/1‑148 (MQ=255)
                                                                                         cAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTAc                                                                                                                        >  7:120854/1‑84 (MQ=255)
                                                                                         cAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTAc                                                                                                                        <  8:120854/84‑1 (MQ=255)
                                                                                              ttGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGCAAGtt                                                                                                                                             <  8:127181/58‑1 (MQ=255)
                                                                                              ttGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaaccaa                                                  <  6:101181/149‑1 (MQ=255)
                                                                                                ggCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACcaac                                                 <  6:179869/148‑1 (MQ=255)
                                                                                                 gCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACcaacaa                                               >  2:431498/1‑149 (MQ=255)
                                                                                                       cacaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCg                                         >  6:20578/1‑149 (MQ=255)
                                                                                                       cacaACGCATATACCACAGTGGTGACATGTAAAACCCCGACGAAACTTGCCTCCGGCAATACTATTTACATGCCTGTGCCCGCTTCAGGCGGTCCGGGCCCCATCAGCTCAGTAGCAACGCTTACAATCCCCAAACCAACCAGCAcccg                                        >  3:219087/1‑147 (MQ=255)
                                                                                                        acaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCgg                                        <  1:425111/149‑1 (MQ=255)
                                                                                                        acaACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCgg                                        <  1:431498/149‑1 (MQ=255)
                                                                                                          aaCGCATATACCACAGTGGGTGACATGATAAAGACGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCCCAATCAGCTCTGTAGCAACGATAACAATCAACAACCCAACAAGCACCGGTa                                      >  6:23323/1‑149 (MQ=255)
                                                                                                              cATATACCACATTGGTTGAAATTATAAAGCCGGAGGAAACTTTCCTCCTGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGa                                  <  5:426260/149‑1 (MQ=255)
                                                                                                               atatACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACcaac                                                 <  3:97623/133‑1 (MQ=255)
                                                                                                                tataCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAgg                                >  5:399110/1‑149 (MQ=255)
                                                                                                                  taCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGt                               <  8:302992/148‑1 (MQ=255)
                                                                                                                   aCCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGc                             >  7:304941/1‑149 (MQ=255)
                                                                                                                     cacaGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCTCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGt                           >  6:208209/1‑149 (MQ=255)
                                                                                                                     cacaGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaaccaa                                                  >  6:383871/1‑126 (MQ=255)
                                                                                                                     cacaGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCCCCAAACCAACAAGCACCGGTACCGAGGTGCGt                           >  3:430710/1‑149 (MQ=255)
                                                                                                                     cacaGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGt                           >  5:449265/1‑149 (MQ=255)
                                                                                                                      acaGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGtt                          <  8:304941/149‑1 (MQ=255)
                                                                                                                       caGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGttt                         <  3:28904/149‑1 (MQ=255)
                                                                                                                         gTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTa                       <  2:247625/149‑1 (MQ=255)
                                                                                                                         gTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTa                       <  6:449265/149‑1 (MQ=255)
                                                                                                                          tGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaac                                                     <  3:81072/118‑1 (MQ=255)
                                                                                                                          tGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCaccaaac                                                     >  4:81072/1‑118 (MQ=255)
                                                                                                                            ggTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGa                    <  6:399110/149‑1 (MQ=255)
                                                                                                                              tGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACt                  <  3:319234/149‑1 (MQ=255)
                                                                                                                                aCATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGtttt                        <  1:225179/141‑1 (MQ=255)
                                                                                                                                aCATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGtttt                        >  2:225179/1‑141 (MQ=255)
                                                                                                                                aCATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACtt                 >  3:118047/1‑148 (MQ=255)
                                                                                                                                  aTGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCa               >  7:130265/1‑148 (MQ=255)
                                                                                                                                  aTGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCa               >  4:5157/1‑148 (MQ=255)
                                                                                                                                    gATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTAc                                     <  6:267876/124‑1 (MQ=255)
                                                                                                                                    gATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTAc                                     >  5:267876/1‑124 (MQ=255)
                                                                                                                                        atgCCGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGc         >  3:195558/3‑148 (MQ=255)
                                                                                                                                        aaGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACTAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGc         <  5:242003/148‑1 (MQ=255)
                                                                                                                                         tgCCGGAGGCAAGTTTCCTCCGGCTTTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGc         >  5:152449/2‑147 (MQ=255)
                                                                                                                                         aGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGc         >  1:94812/1‑147 (MQ=255)
                                                                                                                                         aGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCga       <  5:20578/149‑1 (MQ=255)
                                                                                                                                           ccGGAGGCAAGTTTCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCgaga     <  4:195558/149‑1 (MQ=255)
                                                                                                                                           ccGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGc         >  7:125130/1‑145 (MQ=255)
                                                                                                                                             ggaggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCCCCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGCCTTCAAACGGCGAGATc   >  5:387449/1‑149 (MQ=255)
                                                                                                                                              gaggaAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCt  >  5:417280/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CAGGCACTGACCATCAGCCGCACCACTGCCAGGCGTTATCTTGAATATTGCGCCAGCCGCCATCTGATTATTGCTGAAATTGTTCACGGCAAAGTTGGCAGACCACAACGCATATACCACAGTGGGTGACATGATAAAGCCGGAGGAAACTTGCCTCCGGCATTACTATTTACTTGCCTGTGACCGCTGCTGCCGTTCCTGGCACCATCAGCTCTGTAGCAACGATAACAATCACCAAACCAACAAGCACCGGTACCGAGGTGCGTTTTACGACTTCAAACGGCGAGATCT  >  NC_000913/654414‑654704

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