Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 485,845 (G)5→6 coding (85/648 nt) acrR → transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913485,8401.G98.3% 212.4 / ‑3.1 60coding (80/648 nt)acrRtranscriptional repressor
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (39/20);  total (39/21)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

TCGTCGTGCTATGGTACATACATTCACAAATGTATGTAAATCTAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑A‑GGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAACTGTCAGAATCCAATATTGGTGAACTAGAGC  >  NC_000913/485695‑485983
                                                                                                                                                   |                                                                                                                                               
tcgtcgTGCTATGGTACATACATTCACAAATGTATGTAAATCTAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑A‑ggggg                                                                                                                                            >  2:634023/1‑151 (MQ=255)
tcgtcgTGCTATGGTACATACATTCACAAATGTATGTAAATCTAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑A‑ggggg                                                                                                                                            >  2:365348/1‑151 (MQ=255)
       gCTATGGTACATACATTCACAAATGTATGTAAATCTAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTatcat                                                                                                                                      >  2:547614/1‑151 (MQ=255)
           tGGTACATACATTCACAAATGTATGTAAATCTAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCAc                                                                                                                                  >  2:336998/1‑151 (MQ=255)
               acatacatTCACAAATGTATGTAAATCTAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTc                                                                                                                               <  2:174738/150‑1 (MQ=255)
               acatacatTCACAAATGTATGTAAATCTAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCg                                                                                                                              >  2:331351/1‑151 (MQ=255)
                catacatTCACAAATGTATGTAAATCTAACGCCTGTAAATTCACGAACATATTGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGc                                                                                                                             >  1:19867/1‑151 (MQ=255)
                     atTCACAAATGTATGTAAATCTAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGc                                                                                                                        <  2:608914/151‑1 (MQ=255)
                      ttCACAACTGTATGTAAATCTATCGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGCGG‑GGGGGTATCATCCACCTCGCTGGGc                                                                                                                        <  2:654292/150‑1 (MQ=255)
                      ttCACAAATGTATGTAAATCTAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCg                                                                                                                       >  1:85671/1‑151 (MQ=255)
                       tCACAAATGTATGTAAATCTAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGAGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCg                                                                                                                       <  2:691469/150‑1 (MQ=255)
                          caAATGTATGTAAATCTAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCg                                                                                                                       >  1:459702/1‑147 (MQ=255)
                                 atgtaAATCTAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCaaa                                                                                                             >  2:523248/1‑150 (MQ=255)
                                   gtaAATCTAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAagc                                                                                                          >  2:757201/1‑151 (MQ=255)
                                        tCTAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTg                                                                                                     >  2:257200/1‑151 (MQ=255)
                                          tAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGc                                                                                                   >  1:73302/1‑151 (MQ=255)
                                          tAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGc                                                                                                   >  1:575674/1‑151 (MQ=255)
                                           aaCGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCg                                                                                                  >  1:317638/1‑151 (MQ=255)
                                                      aTTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGt                                                                                       <  1:547577/151‑1 (MQ=255)
                                                      aTTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGt                                                                                       >  2:691344/1‑151 (MQ=255)
                                                              aCATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCt                                                                                >  2:244399/1‑150 (MQ=255)
                                                              aCATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCt                                                                                >  2:780615/1‑150 (MQ=255)
                                                               cATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTAc                                                                              >  1:134210/1‑151 (MQ=255)
                                                                atatGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACt                                                                             >  2:55684/1‑151 (MQ=255)
                                                                atatGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACt                                                                             >  1:142835/1‑151 (MQ=255)
                                                                  atGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTgg                                                                           <  1:468996/151‑1 (MQ=255)
                                                                  atGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTgg                                                                           >  1:736693/1‑151 (MQ=255)
                                                                      cACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCAtt                                                                       >  1:717775/1‑151 (MQ=255)
                                                                       aCGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCAttt                                                                      >  1:112760/1‑151 (MQ=255)
                                                                       aCGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCAttt                                                                      <  2:670460/151‑1 (MQ=255)
                                                                          aaaaaCCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTaa                                                                   <  2:376019/151‑1 (MQ=255)
                                                                           aaaaCCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTaaa                                                                  <  1:783243/151‑1 (MQ=255)
                                                                           aaaaCCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTaa                                                                   >  1:269887/1‑150 (MQ=255)
                                                                              aCCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGAc                                                               >  2:338569/1‑151 (MQ=255)
                                                                                   aCAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTc                                                          <  1:677132/151‑1 (MQ=255)
                                                                                     aagaagCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCgg                                                        >  2:542152/1‑151 (MQ=255)
                                                                                       gaagCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGAt                                                      >  2:742638/1‑151 (MQ=255)
                                                                                         agCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGAttt                                                    <  1:563795/151‑1 (MQ=255)
                                                                                              aaGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCa                                               <  2:98248/151‑1 (MQ=255)
                                                                                                 aaaCGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTg                                            >  1:749365/1‑151 (MQ=255)
                                                                                                       gcCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCt                                      <  2:405857/151‑1 (MQ=255)
                                                                                                       gcCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCt                                      >  1:399290/1‑151 (MQ=255)
                                                                                                        ccAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTg                                     >  2:594008/1‑151 (MQ=255)
                                                                                                          aaCACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTggg                                   >  1:120460/1‑151 (MQ=255)
                                                                                                            cacaTCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGaa                                 <  2:194318/151‑1 (MQ=255)
                                                                                                             acaTCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAAc                                >  1:567829/1‑151 (MQ=255)
                                                                                                              caTCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAACt                               <  1:277454/151‑1 (MQ=255)
                                                                                                                    cGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAACTGTc                            >  2:560027/1‑148 (MQ=255)
                                                                                                                    cGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAACTGTc                            >  2:410518/1‑148 (MQ=255)
                                                                                                                    cGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAACTGTCa                           <  2:181283/149‑1 (MQ=255)
                                                                                                                    cGATGTGGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAACTGTCAg                          <  2:317662/150‑1 (MQ=255)
                                                                                                                        gtgGCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAACTGTCAGAATc                      <  1:732522/150‑1 (MQ=255)
                                                                                                                           gCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGGGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAACTGTCAGAATCCAat                  >  1:267922/1‑151 (MQ=255)
                                                                                                                           gCTCTACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAACTGTCAGAATCCAat                  <  2:409755/151‑1 (MQ=255)
                                                                                                                              ctACGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCAAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAACTGTCAGAATCCAatat                >  1:412440/1‑150 (MQ=255)
                                                                                                                                aCGTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAACTGTCAGAATCCAATATTg              <  2:63162/150‑1 (MQ=255)
                                                                                                                                  gTCTTTTCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAACTGTCAGAATCCAATATTGGt            >  2:90110/1‑150 (MQ=255)
                                                                                                                                      tttCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAACTGTCAGAATCCAATATTGGTGaa         <  1:310184/149‑1 (MQ=255)
                                                                                                                                      tttCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAACTGTCAGAATCCAATATTGGTGAACt       >  1:32916/1‑151 (MQ=255)
                                                                                                                                       ttCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAACTGTCAGAATCCAATATTGGTGAACt       >  2:589223/1‑150 (MQ=255)
                                                                                                                                       ttCTCACAGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAACTGTCAGAATCCAATATTGGTGAACt       >  2:268693/1‑150 (MQ=255)
                                                                                                                                           cacaGC‑AGGGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAACTGTCAGAATCCAATATTGGTGAACTAGAGc  >  1:585779/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
TCGTCGTGCTATGGTACATACATTCACAAATGTATGTAAATCTAACGCCTGTAAATTCACGAACATATGGCACGAAAAACCAAACAAGAAGCGCAAGAAACGCGCCAACACATCCTCGATGTGGCTCTACGTCTTTTCTCACAGC‑A‑GGGGGTATCATCCACCTCGCTGGGCGAGATTGCAAAAGCAGCTGGCGTTACGCGCGGTGCAATCTACTGGCATTTTAAAGACAAGTCGGATTTGTTCAGTGAGATCTGGGAACTGTCAGAATCCAATATTGGTGAACTAGAGC  >  NC_000913/485695‑485983

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