Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,947,4811.A37.3% 113.1 / 135.4 67intergenic (+18/‑16)metZ/metWtRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base . (2/40);  new base A (1/24);  total (3/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.92e-01

GGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGT  >  NC_000913/2947331‑2947590
                                                                                                                                                       |                                                                                                              
gggtTAAAAATAAGGTTGCATGAAATATCGAGCGTAGTATATTTCGCATACAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑a                                                                                                               <  4:1363990/151‑1 (MQ=14)
  ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aa                                                                                                              <  6:942957/151‑1 (MQ=255)
  ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aa                                                                                                              <  6:444382/151‑1 (MQ=255)
  ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCGCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aa                                                                                                              <  2:1671127/151‑1 (MQ=255)
   gTAAAACAAAGGTTGCATGATAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAAGCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aa‑a                                                                                                            <  2:780235/151‑1 (MQ=255)
   gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aa‑a                                                                                                            <  8:1485010/151‑1 (MQ=255)
   gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aa‑a                                                                                                            <  1:445316/151‑1 (MQ=255)
   gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aa‑a                                                                                                            <  5:478081/151‑1 (MQ=255)
   gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aa‑a                                                                                                            <  5:1663272/151‑1 (MQ=255)
   catctcacaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aa‑aa                                                                                                           <  7:594153/143‑2 (MQ=17)
    aaaaaataaGGCGCCAGGAAAACGCGAGCGCTGTATAGACCGCCTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aa‑aag                                                                                                          <  7:1029246/150‑3 (MQ=11)
     aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑Ag                                                                                                           <  1:1030889/150‑1 (MQ=255)
     aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑Ag                                                                                                           <  3:1230312/150‑1 (MQ=255)
     aaaaaaaaGGTTGAATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aa‑aagc                                                                                                         <  4:663134/151‑4 (MQ=255)
      aaaacaaGGTTGCATGAAAACGACATCGTAGTATAGTCCGAATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aa‑aagca                                                                                                        <  6:1414229/151‑5 (MQ=14)
      aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCa                                                                                                         <  6:184682/151‑1 (MQ=255)
      aaaaaaaGGCTGCATGAAAACCCCAGCGGAGTATAGTGGGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aa‑aagca                                                                                                        <  4:349828/151‑5 (MQ=17)
       aaataaGGGTGAATGAAAAATCGAGCGGAGTATATTGCGCATACCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aa‑aagcag                                                                                                       <  3:36972/151‑6 (MQ=11)
       aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTGATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAAGCAATTAAAATTTGATGAAGT‑AA‑AGCAg                                                                                                        <  4:909356/151‑1 (MQ=255)
       aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAg                                                                                                        <  8:658582/151‑1 (MQ=255)
        ataaaGGTTGCATGAAAAAGCGATCGGAGTATATTTCGCATCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGt                                                                                                       <  6:147243/149‑1 (MQ=14)
        aaaaaGGTTGCATGAAATCGCGAGCGGAGTATAGTGCGCATCCACGGAGGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAAGCAATTAAAATTTGATGAAGT‑AA‑AGCAGt                                                                                                       <  3:1194493/151‑1 (MQ=255)
        aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAA‑TTGAAGAAAT‑aa‑aa                                                                                                           >  4:596899/1‑145 (MQ=255)
         acaaGGTTGAATGAAAAAGCGAGCGGTGTATACTCCGCCTCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTa                                                                                                      <  8:1019660/149‑1 (MQ=11)
         acaaGGGTGCATGAAATCGCGAGCGGAGTATAGTGCCCATACACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTa                                                                                                      <  3:1350298/149‑1 (MQ=14)
         aaacgggggCATGAAAAAGCTAGCGGAGTATACTGCGCCTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTa                                                                                                      <  8:671395/143‑1 (MQ=11)
         aaaccgctgaaTGAAAACGCGAGCGGTGTATAGTGCGCCTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTa                                                                                                      <  7:942383/141‑1 (MQ=14)
         aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTa                                                                                                      <  2:1324251/151‑1 (MQ=255)
         aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTa                                                                                                      <  4:1129154/151‑1 (MQ=255)
         aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCATCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTa                                                                                                      <  1:1552731/151‑1 (MQ=255)
         aaaaGGTTCAAGGAAAAACCGAGCGCAGTATAGTGCGCATCCCGTGACGCGGGGTGGAGCAGCGTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTa                                                                                                      <  7:1141690/151‑1 (MQ=11)
          aacgcatgaaTGAAAACGCCAGAGGAGTATACTCCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTAc                                                                                                     <  5:745821/142‑1 (MQ=11)
          aaaGGTTGCATGAAAACGCTAGCGGTGTATAGTGGGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTAc                                                                                                     <  7:1140022/151‑1 (MQ=14)
          aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTAc                                                                                                     <  6:1004093/151‑1 (MQ=255)
          aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTAc                                                                                                     <  8:1391663/151‑1 (MQ=255)
           acggcggCATGAAAAAGCGAGCGGAGTATAGAGCGCCTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTACg                                                                                                    <  3:334349/145‑1 (MQ=11)
           aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACg                                                                                                    <  3:1313713/151‑1 (MQ=255)
           aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGAAGCAGTAAA‑AGCAGTACg                                                                                                    >  5:675749/1‑151 (MQ=17)
           aaGGTTGAATGAAAAAGCTTTCGGTGTATATTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTACg                                                                                                    <  4:1001918/151‑1 (MQ=11)
            aGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACgg                                                                                                   <  6:1034437/151‑1 (MQ=255)
            aGGTTGCATGAAAACGAGATTGGTGTATATTGCGCTTCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTACgg                                                                                                   <  6:1283732/151‑1 (MQ=11)
            aGGTTGAATGAATACTCTAGGGTAGTATAGTGCGATTCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTACgg                                                                                                   <  1:677035/151‑1 (MQ=11)
             ggTTGCATGAAATCGCGTGCGTAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTACGGt                                                                                                  <  4:1508555/151‑1 (MQ=14)
             ggTTGCATGAAAACGCGTGGGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTACGGt                                                                                                  <  7:1336663/151‑1 (MQ=14)
             ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGt                                                                                                  <  7:550195/151‑1 (MQ=255)
             ggTTGCAGGAAAACTCCAGCGTAGTATAGTCCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTACGGt                                                                                                  <  4:1304660/151‑1 (MQ=14)
             ggTTGAATGAAAACTCTATCGTAGTATAGTGCGCTTCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTACGGt                                                                                                  <  1:797676/151‑1 (MQ=11)
              gTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGTg                                                                                                 <  6:1181454/151‑1 (MQ=255)
              gTTGAATGAATACGCGATGGGAGTATATTCCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTACGGTg                                                                                                 <  5:115132/151‑1 (MQ=11)
              gTTGAATGAAAACGCGAGCGTAGTATAGTTCGATTCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTACGGTg                                                                                                 <  6:434347/151‑1 (MQ=11)
               gtgaatcaaaaCAAGAGCGCAGTATTTTTCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTACGGTGa                                                                                                <  2:1392769/144‑1 (MQ=11)
                tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGTGa                                                                                                <  5:447827/150‑1 (MQ=255)
                ggaaTGAATACGCGTGCGGAGTAGAGTGCGAATCCACGGACGCGGGGTGGAGCAGCCGGGTAGCTCGGCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTACGGTGa                                                                                                <  4:1435488/147‑1 (MQ=11)
                 gaaTGAAAACTCGATCGGAGTATATTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTACGGTGAcg                                                                                              <  2:1437302/149‑1 (MQ=14)
                 gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAAATTTTCAGGATGA‑AA‑AGCAGTACGGTGAcg                                                                                              >  3:1428889/1‑151 (MQ=18)
                 gCATGAAAACCCTAGCGTTGTATACTTCGCTTACAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTACGGTGAcg                                                                                              <  1:1014255/151‑1 (MQ=11)
                 gCAGGAAAACGCGAGGGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTACGGTGAc                                                                                               <  7:1441672/150‑1 (MQ=14)
                  cATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGTGAcgc                                                                                             <  3:1125768/151‑1 (MQ=255)
                  cATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGTGAcgc                                                                                             <  6:717179/151‑1 (MQ=255)
                   aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGTGAcgcg                                                                                            <  5:259121/151‑1 (MQ=255)
                   aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGTGAcgcg                                                                                            <  3:755436/151‑1 (MQ=255)
                   aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAA‑AGCAGTACGGTGAcgcg                                                                                            <  3:1015069/151‑1 (MQ=255)
                    tGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGTCGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACTGTGAcgcg                                                                                            <  1:1068793/150‑1 (MQ=255)
                      aaaaCTCGAGCTGATTATATTGCTCATTCAAGGACGCGGGGTGGATCATCCTTTTAGCTAGTAGGGCTCATNANCCTNATGNNNTAGGTTTAAATCAGGCACCCTCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGTGACGCgggg                                                                                         <  2:91935/151‑1 (MQ=12)
                       aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGTGACGCGGGGt                                                                                        <  2:18448/151‑1 (MQ=255)
                       aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGTGACGCGGGGt                                                                                        <  4:892280/151‑1 (MQ=255)
                       aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGGCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGTGACGCGGGGt                                                                                        <  5:1601863/151‑1 (MQ=255)
                        aaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AC‑AGCAGTACGGTGACGCGGGGTg                                                                                       <  2:848779/151‑1 (MQ=255)
                        aaCGCGAGCGGAGTATAGTGCGCAGCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGTGACGCGGGGTg                                                                                       <  4:148788/151‑1 (MQ=255)
                            cgAGCGGAGTATAGTGAGCATCCACGGACGCTTGGTGGAGCAGCCTGGTAGCTCGTCGGTCTCATAACCCGAATGTCGTCTGTTCAAATCCTGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCATTACGGTGACGCGGGGTGGagc                                                                                   <  2:687331/151‑1 (MQ=18)
                             gAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGTGACGCGGGGTGGagca                                                                                  <  8:1483273/151‑1 (MQ=255)
                              cgCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGGCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCGCCGCAACCAATTAAAATTTGATGAAGT‑AA‑CGCAGTACGGTGACGCGGGGTGGagcag                                                                                 <  1:978133/150‑1 (MQ=255)
                              aggggAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCGTGGTAGCTCGTGGGGCTCATAAGCCGAAGGTCGGCGGTTCATATCCGGCCCCGGCAACCAATTAAAAATTGATGAAGT‑TA‑TGCAGTACGGTGACGCGGGGTGGagcag                                                                                 <  7:1041606/148‑1 (MQ=18)
                              aGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGTCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGTGACGCGGGGTGGagca                                                                                  <  2:546395/150‑1 (MQ=255)
                              aGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCGATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGTGACGCGGGGTGGagcag                                                                                 <  5:1513867/151‑1 (MQ=255)
                                 ggtgTATAGTGCGCATCCTCGGACGCTGGTTGGAGCAGCGTGTTAGCTGTGCTGTCTCATAAGCCTAAGGTGGTCGGTTGAAATCCTGCCCCCGCAATCAATTAAAATTTGAAGAATT‑TA‑AGCATTACCGTGACGCGGGGTGGAGCAGCCt                                                                              <  1:799572/148‑1 (MQ=12)
                                 ggAGTATAGTGAGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTGCAAATCCCGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGTGACGCGGGGTGGAGCAGCCt                                                                              <  3:59725/151‑1 (MQ=255)
                                 gcaGTATAGTGATGATCCACCGACGCGGGGATGAGCAGCCTGTTATCGCGTCTGGATCATAAACCGAAGGTCTTCGGTTCAAATAACCAAGCCGGAAACAATTAAAATTAGATGTATT‑TA‑CGCATTACTCTGACGCGGGGTGGAGCAGCCt                                                                              <  8:1572103/149‑1 (MQ=9)
                                                               tGGAGTAGCCTTGTAGCTCGGCTGTATCATAACCCGAATGTCGTCTGTTGAAAACCTTCCGCCGCAAGGAATTACACATTGTTGAATT‑CA‑TGCATTCCTGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAgg                                                 <  1:1489941/150‑1 (MQ=2)
                                                                                                           cgGTTCAAATACGGCCGGCGGAACCAATTAAACAGAGATGAAGT‑AT‑CCCACTCCCGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGaa    <  5:67516/151‑1 (MQ=14)
                                                                                                             gtgcgaaacagcccccgGCCACCAATTACAAAGGGATGAAGA‑GA‑CGCAGTACGGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGGCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAgc  <  5:550552/139‑2 (MQ=14)
                                                                                                             gTTCCAAACCGCCCCGCGCAACCAATAACAATGCGCTGAAGT‑AACCGCA‑TACCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGt  <  5:550547/151‑1 (MQ=14)
                                                                                                                             cGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATc                                  >  7:1303961/1‑103 (MQ=255)
                                                                                                                                                       |                                                                                                              
GGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AA‑AGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGT  >  NC_000913/2947331‑2947590

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