Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,916,216 G→A G387E (GGG→GAG)  barA → sensor histidine kinase BarA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,916,2160GA97.1% 196.9 / ‑2.2 70G387E (GGG→GAG) barAsensor histidine kinase BarA
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (33/34);  total (35/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.64e-01

TATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGGGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCGGGATACCGGCATTGGTATTCCTGAACGCGATCAATCGCGCTTATTCCAG  >  NC_000913/2915989‑2916454
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               
tATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATAtt                                                                                                                                                                                                                           <  1:18987/250‑1 (MQ=255)
    ggCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAg                                                                                                                                                                                                                      <  1:337817/251‑1 (MQ=255)
          aCCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGCATACGATTGACCGTACGGCAAATAATTCGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGTTTCTGGAAAGTATGCCCTTCCCAATACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATACTTCTCCCGATAAAGAGTTAAGACTGAAGCTCACTATTAAAAACGACGt                                                                                                                                                                                                                >  2:93540/1‑251 (MQ=255)
            ccgccagaagcagAAAACAGAATTAACACCAACTCAGAGCGATCACCTGAATACGATTGAACGTGCGGCAAATAATTTGCTGGCAATTAGTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGGTAGAACTGACGCTCAATATTAAAAGCGACGTGc                                                                                                                                                                                                              <  2:186252/239‑1 (MQ=255)
            ccGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTAAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGCAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGGGtt                                                                                                                                                                                                                                              >  2:290678/1‑219 (MQ=255)
            ccGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTAAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGGGtt                                                                                                                                                                                                                                              <  1:290678/219‑1 (MQ=255)
                  gactgtGAAAACAGAATTTACTCCAACGCAGCGCGATCACCGGAATACGATTGAACGTTCGGCAAATAAGTTGCTGGCAATTATTAATGATGTTCTCGTCTTCTCGAAACTGGAAGCATGTAAGCTGATTCTTGAAAGTATTCCATTCCCACTACGCAGCACGCTGTATGAAGTCGTTTCTCTGCTGGTACAGTCTTCTCACGATTAAGAGTTAGAACTGACGCTCAATATTTAAAGCGACGTgtatgata                                                                                                                                                                                                        <  2:656454/246‑8 (MQ=255)
                          aaaCAGAATTAACACAAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCACATAATTGGCTGGCAAGTATTAATGATGGTCGCGACTTCTCTAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACTCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATc                                                                                                                                                                                                <  2:620481/251‑1 (MQ=255)
                             cAGAATTAACACCATCGCAGCGCGATCAACTGAATACGATTGAACGTCCGGCAAAGAATTTGCTGGTAATTATTAATGATGTTCTCGACTTCTCGAAACTTGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGGCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTATAACTGACGCTCAATATTAAAAGCGACGTGCCTGATACCGTGATCGGc                                                                                                                                                                                             <  2:42901/251‑1 (MQ=255)
                                   tAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTGCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGAccc                                                                                                                                                                                        <  1:636535/250‑1 (MQ=255)
                                          agcGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGTCTTCTCGAAACTGGAAGCAGGTAAGCTTATGCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGTAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGAt                                                                                                                                                                                <  2:489455/249‑1 (MQ=255)
                                              cAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCGGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGTTTCTGGAAAGTATTCCATTCCCACTCCGCAGCACGCTGCATGCAGTCGTTACTCTGCTGGCACATTCTTCTCACGCTAAAGTGTTATAACTGACGCTCAATATTACATGCGACGTGCCTGATACCGTGATCGGCGACCCGCTGCGATTAc                                                                                                                                                                             >  2:171126/1‑250 (MQ=255)
                                               agctagATCACCTGAATACGATTGAAAGTTCGGCAAATAATTTGATGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGATGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACTCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACTATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAg                                                                                                                                                                           <  2:490348/246‑1 (MQ=255)
                                                 cgcgACCACCTGAATACGATTGAACAGTCGGCATATAATATGCTGTCACTTATTAATGATGTTCTCGACTTCTCGAAACTAGAAGCAGGTACGCTGATTCTGGAAAGAATTCCATTCCCACTACGCAGCACGCTGGATGAGGTCGTTACTCTGCTGGCACATACTTCTCACGATAAAAAGTTAGGACTGACGCTCAAAATTAAAAGCGACGTGCCTGATAACGTGaacagagaccagccccgcttacagc                                                                                                                                                                          >  2:565324/1‑226 (MQ=255)
                                                   cgATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGAGGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGGATTCCATTCCCACTACGCAGCACGCGGGATGAAGTCGTTACTCTGCTGGCGCATTCTTCTCACGATAAAGAGTTAGAACTGACTCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCaaa                                                                                                                                                                       <  1:442816/251‑1 (MQ=255)
                                                        aCCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAatcatc                                                                                                                                                                  >  1:76984/1‑251 (MQ=255)
                                                             aaTACGATTGAACCTTAGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGATAGTATTCCATTCCCACTACGAAGAACGCTGGATGAAGTCGTTACTCAGCTGGCACATTCTTCTCACGCTAAAGAGTTAGAACTGCCGCTCAATATTAATCGCGAAGGGCCTGCTACAGTACTCGGCGCCCCGCTGCGATTACAGCAAGTCCTCACTa                                                                                                                                                              >  2:392309/1‑250 (MQ=255)
                                                              aTACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAAc                                                                                                                                                            >  1:121622/1‑251 (MQ=255)
                                                                 cGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCtg                                                                                                                                                         >  1:324702/1‑251 (MQ=255)
                                                                  gATTGAACGGTCGGCAAATAATTTGCTGGCAAGTATTAATGATGTTCGCTACTTCTCGAAACTTGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACTCGCTGCTATTACAGCAAATCATCACTAACCtgg                                                                                                                                                        <  2:93092/251‑1 (MQ=255)
                                                                   aTTGAACGTTCGGCAAATAATTTGCTGGAAATTATTCATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCtggg                                                                                                                                                       >  2:229616/1‑250 (MQ=255)
                                                                        aCGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTgggg                                                                                                                                                   >  1:173973/1‑250 (MQ=255)
                                                                            tCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATg                                                                                                                                               <  1:386274/250‑1 (MQ=255)
                                                                            tCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATg                                                                                                                                               <  1:481371/250‑1 (MQ=255)
                                                                             cGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGc                                                                                                                                              >  1:258473/1‑250 (MQ=255)
                                                                             cGGCAAATAAGTTGCTGGCAATTTGTAATGATGGTCTCTACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTACCCTGGTGGGGAATGCa                                                                                                                                             <  2:274869/251‑1 (MQ=255)
                                                                                      atTTGCTGGCAATTATTCATAATGTTCTCGACTTATCGAAACTGGAAGCAGGTAAGCTGATTCTGGAACGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAACGAGTTAGAACTGACGCTCAATAGTAAAAGCGACGTGCCTGATAACGTGATCGGCGCCCCGCTGCGCTTACAACCAATCATCACTAACCTGGTGGGGCATGCAATTAAATTc                                                                                                                                    >  2:591847/1‑251 (MQ=255)
                                                                                      atTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTc                                                                                                                                    >  1:119723/1‑251 (MQ=255)
                                                                                        ttGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCa                                                                                                                                   <  1:275962/250‑1 (MQ=255)
                                                                                          gCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCg                                                                                                                                <  2:154816/251‑1 (MQ=255)
                                                                                          gCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGGCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCg                                                                                                                                <  2:274824/251‑1 (MQ=255)
                                                                                            tGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAata                                                                                                                                                                                                                             <  2:509328/156‑1 (MQ=255)
                                                                                            tGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAata                                                                                                                                                                                                                             >  1:509328/1‑156 (MQ=255)
                                                                                                 attattAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATgg                                                                                                                         <  1:101002/251‑1 (MQ=255)
                                                                                                   tattaATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGc                                                                                                                        >  1:136425/1‑250 (MQ=255)
                                                                                                   tattaATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTAGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGATGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGc                                                                                                                        >  1:336817/1‑250 (MQ=255)
                                                                                                                ctcGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGACCTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGACAACATTTATATTCt                                                                                                         >  2:358759/1‑251 (MQ=255)
                                                                                                                  cGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTgg                                                                                                       >  1:165744/1‑251 (MQ=255)
                                                                                                                  cGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTgg                                                                                                       <  1:7165/251‑1 (MQ=255)
                                                                                                                    aCTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAATGATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTa                                                                                                     <  2:119723/251‑1 (MQ=255)
                                                                                                                       tctcGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGa                                                                                                   >  1:370513/1‑250 (MQ=255)
                                                                                                                            aaactgtcaacaagTAAGCTGATTCTGGAAAGTATTACATTCCCACTACGCAGCAAGCTGGAGGAAGTCGTTACTCGGATGGCACATTCTTCTCCCGATAAAGAGTTAGAACTGACAGTCAATATTAAAAGCGACGTGCCTGATAACGCGATCGGCGCCCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGCATGCAATTAAATTCACCGTGGATGGCCAACTTGATATTTTGGTAGaaatatg                                                                                             >  2:361152/14‑247 (MQ=255)
                                                                                                                            aaaCTGGAAGCAGGTAAGCTGATTCTGGTAAGTATTCCTTTCCCACTACGCAGCTCGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACg                                                                                             <  2:538670/251‑1 (MQ=255)
                                                                                                                              aCTGGAAGCAGGTAAGCTGAGTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTg                                                                                           <  2:136425/251‑1 (MQ=255)
                                                                                                                               cTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTAACACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGTATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAAAGTgc                                                                                          >  1:604974/1‑251 (MQ=255)
                                                                                                                                 ggAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTgc                                                                                          <  1:444116/249‑1 (MQ=255)
                                                                                                                                  gAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGACACATTGGTATTCTGGTAGaaaaccgatcgat                                                                                       >  2:322517/1‑242 (MQ=255)
                                                                                                                                            aGCTGAATATGGAAAGTATTCCATTCCCACTACGCAGCACGCTAGATGAAGTCGTTACTATGCTGGCACATACTTCTCACGATACAGAGTTAGAACTGACGCTCAATATTAAAAGCGCCGAGCCTGATAACGTGATCGGCGACCCGATACGATTATAGCAAATCATACCTAACCTGGTGGGGAATGCGTTTACATTAACCGAGACGGGCAACATTAATATTCTGAGAGAAAAACGTGCGCTCGGTAATAcc                                                                             >  2:417404/1‑251 (MQ=255)
                                                                                                                                                   tCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGc                                                                      >  1:27398/1‑251 (MQ=255)
                                                                                                                                                         aaGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTg                                                                >  1:504820/1‑251 (MQ=255)
                                                                                                                                                            tacgccattccCACTACGCAGCACGCTGGATGAAGTCGTTACTATGCAGTCACATTCTTCTCACGCTAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGCCGTGCCTGATAACGTGATCGGCGACCCGCTTAGATTACAGCAAATCATCACTAACCTGGTGGGGCATGCAATTCAATTCACCGAGAAAGGCAACATTCCTCTTCTGTGAGACAAAAGTTCGCTGAGTACTAGCAAAATGCGGATTGatg                                                             >  2:372013/5‑249 (MQ=255)
                                                                                                                                                            tATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGa                                                               >  1:485192/1‑249 (MQ=255)
                                                                                                                                                                ccgttccCACTACGCAGCACGCTGGATGAAGTTGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCa                                                         <  2:216419/248‑1 (MQ=255)
                                                                                                                                                                  attccCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGc                                                          >  1:615519/1‑248 (MQ=255)
                                                                                                                                                                   ttccCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATAcccaa                                                                          >  2:668584/1‑228 (MQ=255)
                                                                                                                                                                   ttccCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCaaa                                                                          <  1:668584/231‑1 (MQ=255)
                                                                                                                                                                       cACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACCTTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGCCCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGGTATTCTGTTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCgg                                                  >  2:285508/1‑251 (MQ=255)
                                                                                                                                                                        aCTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTc                                                                                                          >  2:513970/1‑194 (MQ=255)
                                                                                                                                                                        aCTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAAGGGCAACATTGATATTc                                                                                                          <  1:513970/194‑1 (MQ=255)
                                                                                                                                                                           aCGCAGCACGCTGGATGAAGTCGGTACGCTGCTGGCACATTCTTCTCACGATAAAGAGTGAGAACTGACGCTCAAGATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCGGGATa                                              <  1:392309/251‑1 (MQ=255)
                                                                                                                                                                              cagcaCGCTGGTTGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCGGGATACCg                                           >  1:282436/1‑251 (MQ=255)
                                                                                                                                                                                                  ttacccttcaggtacaTTATTCTCACGATAAAGAGTTAGAACGTACGCTCAATATTAAAAGCGACGTGCCGGATAACGTGATCGTCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAAGTCACCGAGAATGGCAACATGGATATTCTGGTAGAAAAACGTGCGCCGAGTAATACCAAAGTGCAGATTGAAGTTCAGATTCTGGATTCCGGCATTGGTATTCTTGAAc                         <  2:452649/236‑1 (MQ=255)
                                                                                                                                                                                                        tgctgGCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCAGGATACCGGCATTGGTATTCCTGAACGCGAtaa                  >  2:229774/1‑248 (MQ=255)
                                                                                                                                                                                                            ggCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCGGGATACCGGCATTGGTATTCCTGAACGCGATCAATcg              <  1:136080/250‑1 (MQ=255)
                                                                                                                                                                                                             gCACATTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCGGGATACCGGCATTGGTATTCCTGAACGCGATCAATcgcg            <  2:258473/251‑1 (MQ=255)
                                                                                                                                                                                                                 aTTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCATAGTGCAGATTTAAGTGCAGATTCGGGATACCGGCATTGGTCTTCCTGAACGCGATCAATCGCGCtt         >  2:585922/1‑250 (MQ=255)
                                                                                                                                                                                                                 aTTCTTCTCACGATAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCGGGATACCGGCATTGGTATTCCTGAACGCGATCAATCGCGCTTa        <  1:3086/251‑1 (MQ=255)
                                                                                                                                                                                                                         cACGATAAAGAGTGAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAAGCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGAGTCGGGATACCGGCATTGGTATTCCTGAACGCGATCAATCGCGCTTATTCCAg  <  1:229616/249‑1 (MQ=255)
                                                                                                                                                                                                                              tAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGatat                                                                                                            <  2:629702/138‑1 (MQ=255)
                                                                                                                                                                                                                              tAAAGAGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGatat                                                                                                            >  1:629702/1‑138 (MQ=255)
                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                               
TATTGGCTTTACCCGCCTGACGCTGAAAACAGAATTAACACCAACGCAGCGCGATCACCTGAATACGATTGAACGTTCGGCAAATAATTTGCTGGCAATTATTAATGATGTTCTCGACTTCTCGAAACTGGAAGCAGGTAAGCTGATTCTGGAAAGTATTCCATTCCCACTACGCAGCACGCTGGATGAAGTCGTTACTCTGCTGGCACATTCTTCTCACGATAAAGGGTTAGAACTGACGCTCAATATTAAAAGCGACGTGCCTGATAACGTGATCGGCGACCCGCTGCGATTACAGCAAATCATCACTAACCTGGTGGGGAATGCAATTAAATTCACCGAGAATGGCAACATTGATATTCTGGTAGAAAAACGTGCGCTGAGTAATACCAAAGTGCAGATTGAAGTGCAGATTCGGGATACCGGCATTGGTATTCCTGAACGCGATCAATCGCGCTTATTCCAG  >  NC_000913/2915989‑2916454

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