Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913599,5760TG44.7% 139.1 / 35.5 85V288G (GTG→GGG) cusAcopper/silver efflux system, membrane component
Reads supporting (aligned to +/- strand):  ref base T (4/43);  new base G (38/0);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 5.37e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACA  >  NC_000913/599428‑599722
                                                                                                                                                    |                                                                                                                                                    
aCGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGggtgg                                                                                                                                                    <  2:54786/151‑1 (MQ=255)
 cGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggg                                                                                                                                                    >  2:920460/1‑150 (MQ=255)
   cTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggg                                                                                                                                                  >  2:1388524/1‑148 (MQ=255)
     cgacgaCTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGa                                                                                                                                                 <  2:340256/149‑1 (MQ=255)
        cgacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTcc                                                                                                                                              >  3:723238/1‑149 (MQ=255)
         gacTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATcc                                                                                                                                              >  2:1470362/1‑148 (MQ=255)
            tttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCg                                                                                                                                          >  2:1055946/1‑149 (MQ=255)
             ttAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCCGCGATc                                                                                                                                       >  2:518402/1‑151 (MQ=255)
              tAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGgggggggggaccttcgggcc                                                                                                                                      >  3:718499/1‑137 (MQ=255)
               aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATcc                                                                                                                                      <  2:1457527/150‑1 (MQ=255)
               aaTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑g                                                                                                                                    <  3:26950/151‑1 (MQ=255)
                aTCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGATCCTTCGATcccg                                                                                                                                    >  2:657780/1‑149 (MQ=255)
                  cacaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGc                                                                                                                                           >  2:543244/1‑142 (MQ=255)
                   acaTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTccctccatcc                                                                                                                                      >  3:122299/1‑138 (MQ=255)
                     aTCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCaaaa                                                                                                                              <  2:1159119/151‑1 (MQ=255)
                        gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACg                                                                                                                           <  2:415127/151‑1 (MQ=255)
                        gTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCaaaaa                                                                                                                             >  2:647622/1‑149 (MQ=255)
                          tttAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGGTCC‑GGCaaaaa                                                                                                                             >  2:170752/1‑147 (MQ=255)
                            tAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGccc                                                                                                                        <  2:1267838/150‑1 (MQ=255)
                             aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCg                                                                                                                       <  2:1014828/150‑1 (MQ=255)
                             aaaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCC‑GGCAAAAAAGccc                                                                                                                        >  3:211646/1‑149 (MQ=255)
                              aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTTCCCCGATCC‑GGCAAAAAACCCCg                                                                                                                       >  2:1447479/1‑149 (MQ=255)
                              aaaGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCAAAAAAGccc                                                                                                                        >  2:570156/1‑148 (MQ=255)
                                  ccAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGGTCC‑GGCCAAAACAcccccgaa                                                                                                                   >  3:1583949/1‑144 (MQ=255)
                                   cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGGTCCTGCGATCC‑CGCGAAAAAGCCCCCGAAGTg                                                                                                                >  3:194217/1‑151 (MQ=255)
                                   cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTg                                                                                                                <  3:678664/151‑1 (MQ=255)
                                   cAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTg                                                                                                                <  2:471034/151‑1 (MQ=255)
                                     gTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGAt                                                                                                              <  3:344825/151‑1 (MQ=255)
                                           aTGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGAGCCTGCGATCC‑GGCCAAAAAGCCCGAGAAGTGGTcgccgc                                                                                                        >  3:809468/1‑151 (MQ=255)
                                            tGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCC‑GGCAAAAAAGCCCGCGAAGGGATCgccgcc                                                                                                       >  2:996127/1‑151 (MQ=255)
                                             ggCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGccgccg                                                                                                      <  3:882363/151‑1 (MQ=255)
                                               cGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGTGGGCTTGGCATCC‑GGCAAAAACGCCCGAGGAGAGGTGCccgccggg                                                                                                    >  3:574890/1‑149 (MQ=255)
                                                  gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGGGGGGGGGGGATCCCGCGGTCC‑CGCAAAAACGCCCGCGGAGGGAGCgccgcc                                                                                                       >  3:207350/1‑145 (MQ=255)
                                                  gCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGaa                                                                                                  <  3:802296/150‑1 (MQ=255)
                                                    ccGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGa                                                                                               <  3:1403250/151‑1 (MQ=255)
                                                      gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑Ca                                                                                            <  3:212591/151‑1 (MQ=255)
                                                      gTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑Ca                                                                                            <  3:705566/151‑1 (MQ=255)
                                                       tttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGTGAAGTGGCGGGCGGGGTGGGGATCCTGCGATCC‑CGCAAAAACGCCCGAGGAGGGATCGCCGCCGGGAAgg                                                                                                >  2:893019/1‑147 (MQ=255)
                                                        ttATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCC‑GGCAAAAAAGCCCGAGGAGTGATCGCCGGCGTGAAGGa                                                                                               >  3:1401774/1‑147 (MQ=255)
                                                         tATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑Caaa                                                                                          <  2:783787/150‑1 (MQ=255)
                                                              gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACt                                                                                        <  2:1306398/147‑1 (MQ=255)
                                                              gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACt                                                                                        <  3:176960/147‑1 (MQ=255)
                                                              gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTg                                                                                       <  3:189528/148‑1 (MQ=255)
                                                              gcgcGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGa                                                                                     <  2:1289486/150‑1 (MQ=255)
                                                                    tGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCt                                                                               <  2:983711/150‑1 (MQ=255)
                                                                     gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGa                                                                             <  2:572542/151‑1 (MQ=255)
                                                                     gTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGATCC‑CGCACAAAAGCCCGCGGAGTGAGCGCCGCCGCGAAGGA‑CAAAATGGaaaaactg                                                                              >  3:13115/1‑145 (MQ=255)
                                                                        gCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGaaaa                                                                          <  2:590181/151‑1 (MQ=255)
                                                                         cGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAag                                                                         <  3:108941/151‑1 (MQ=255)
                                                                           aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAagta                                                                       <  2:938920/151‑1 (MQ=255)
                                                                           aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAagta                                                                       <  2:1360021/151‑1 (MQ=255)
                                                                           aaGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGATCCGGGCAAAAACCCCCG‑GGAGGGATCGCCGCCGGGGAGGA‑GAAAAAGGGAAAGCTgggaa                                                                          >  3:47534/1‑144 (MQ=255)
                                                                             ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAagtagt                                                                     <  3:334889/151‑1 (MQ=255)
                                                                             ggTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGA‑CAAACTGGAAACGCTGAAAagtag                                                                      >  2:744509/1‑150 (MQ=255)
                                                                              gTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGTCCTGCGCTCC‑GGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGA‑CAAACCGGAAACGCTGaaa                                                                           >  3:367148/1‑144 (MQ=255)
                                                                                 cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCCAAAACGCCCGAGCAGTGGTCGCCGCCGGGAAGGA‑CAAACTGGAAACGCTGaaaa                                                                          >  2:1339543/1‑142 (MQ=255)
                                                                                 cAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGc                                                                 >  2:849394/1‑151 (MQ=255)
                                                                                   gATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCg                                                               <  3:833518/151‑1 (MQ=255)
                                                                                    aTTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGGGCCCGGGAACC‑GGGCAAAAAACCCGGGAAGGGAACGCCGCCGGGGAGGAGCAAAGGGGAAA‑GCTGAAAaggagtgagcccg                                                              >  3:843877/1‑144 (MQ=255)
                                                                                     ttGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCgg                                                              <  2:631966/150‑1 (MQ=255)
                                                                                       ggCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAg                                                           <  2:1429140/151‑1 (MQ=255)
                                                                                        gCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAgg                                                          <  2:618219/151‑1 (MQ=255)
                                                                                         cccGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGc                                                         <  3:931356/151‑1 (MQ=255)
                                                                                           cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGt                                                       <  3:191741/151‑1 (MQ=255)
                                                                                           cGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCAGCGGACC‑GGCAAAAACGCCCGCGAAGTGATCGCCGCCCGGTAGGA‑GAAAATGGAAAAGCCGGAAAGAAAACTTCCGGCAGGCgg                                                       >  2:508858/1‑150 (MQ=255)
                                                                                            ggAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTg                                                      <  3:474411/151‑1 (MQ=255)
                                                                                                   cgcCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGGGCCTGCGATCC‑GGCAAAAACGCCCCAGAAGGGATCGGCGCCGGGGAGGA‑CAAAATGGAAAAGCTGGAAAATAGGCTGCCGGAAGGCgggg                                                     >  3:958076/1‑145 (MQ=255)
                                                                                                     ccGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGtt                                             <  3:958916/151‑1 (MQ=255)
                                                                                                     ccGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGGGGAGATAGtt                                             >  2:1306773/1‑151 (MQ=255)
                                                                                                      cGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAAAGCTGGAAAGAAGTCTGCCGGAAGGCGGGGGGATAGtt                                             >  2:780628/1‑150 (MQ=255)
                                                                                                          gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGAGCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGGAGTCTGCCGGAAGGCGTGGAGATAGTTacaa                                         <  3:1137789/150‑1 (MQ=255)
                                                                                                          gCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGGTCC‑GGCAAAAAAGCCCGAGAAGGGATCGCCGCCGGGGAGGA‑CAAACTGGGAAAGGTGAAAAGTAGTCTGCCGGGAGGGGGGGGGATAGTTaaaa                                         >  3:685053/1‑150 (MQ=255)
                                                                                                           cATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTacaa                                         >  2:1615522/1‑149 (MQ=255)
                                                                                                              tGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAAGA‑CAAACTGGAAACGCTGAAAAGGAGTCTGCCGGGAGGCGGGGAGATAGGTACAACATa                                     >  2:1219587/1‑150 (MQ=255)
                                                                                                               gCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGGAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGGAGTCTGCCGGGAGGCGGGGAGATAGTTACAACATa                                     >  2:678266/1‑149 (MQ=255)
                                                                                                                 cGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGGGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAAATGGAAACGCTGAAAAGGAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATAc                                    >  2:1560902/1‑148 (MQ=255)
                                                                                                                    aCTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGc                              <  3:1246885/151‑1 (MQ=255)
                                                                                                                        aaCGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGcagcc                          <  2:279178/151‑1 (MQ=255)
                                                                                                                             cgaaggcgaagTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCAAAAAAGCCCGAGAAGGGATCGCCGCCGTGAAGGA‑CAAAATGGAAACGCTGAAAAGTAGTCTGCCGGGAGGCGTGGGGATAGTTAAAACAAACGATCGCAGCCAGCTc                     >  3:826961/1‑151 (MQ=255)
                                                                                                                              gaaggcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCa                    <  3:1287437/151‑1 (MQ=255)
                                                                                                                              gaaggcgaagTGGCGGGCGGGGGGGGGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGGAAAGGAGTCTGCCGGAAGGCGTGGGGATAGTTACAACATACGATCGCAGCCAGCTc                     >  2:790717/1‑150 (MQ=255)
                                                                                                                                  gcgaagTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGa                <  2:17750/151‑1 (MQ=255)
                                                                                                                                     aagTGGCGGGCGGGGGGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGGGATCGCCGCCGTGAAGGA‑CAAACTGGAAAAGCTGAAAAGTAGTCTGCCGGAAGGCGGGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGAcc              >  3:1308386/1‑150 (MQ=255)
                                                                                                                                       gTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcgcg           <  3:1304492/151‑1 (MQ=255)
                                                                                                                                       gTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACcgc            <  3:1239026/150‑1 (MQ=255)
                                                                                                                                         ggcgggcgGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGGGAAGGA‑CAAACGGGAAACGCTGAAAGGAAGTCTGCCGGAAGGCGTGGGGATAGTTACAACATCCGATCGCAGCCAGCTCATTGCCCGCGcc         >  1:9225/1‑151 (MQ=255)
                                                                                                                                            gggcgggGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATc      >  3:1221686/1‑151 (MQ=255)
                                                                                                                                                ggggTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGAc   >  3:1194807/1‑150 (MQ=255)
                                                                                                                                                ggggTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACa  <  2:1513931/151‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ACGCTCGACGACTTTAATCACATCGTTTTAAAAGCCAGTGAAAATGGCGTGCCCGTTTATCTGCGCGATGTTGCGAAGGTCCAGATTGGCCCGGAGATGCGCCGGGGCATTGCCGAACTAAACGGCGAAGGCGAAGTGGCGGGCGGGGTGGTGATCCTGCGATCC‑GGCAAAAACGCCCGAGAAGTGATCGCCGCCGTGAAGGA‑CAAACTGGAAACGCTGAAAAGTAGTCTGCCGGAAGGCGTGGAGATAGTTACAACATACGATCGCAGCCAGCTCATTGACCGCGCCATCGACA  >  NC_000913/599428‑599722

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