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

GGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGT  >  NC_000913/2947331‑2947590
                                                                                                                                                       |                                                                                                              
gggtTAAAAATAAGGTTGCATGAAATATCGAGCGTAGTATATTTCGCATACAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTa                                                                                                               <  4:1363990/151‑1 (MQ=14)
  ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaa                                                                                                              <  6:942957/151‑1 (MQ=255)
  ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaa                                                                                                              <  6:444382/151‑1 (MQ=255)
  ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCGCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaa                                                                                                              <  2:1671127/151‑1 (MQ=255)
   gTAAAACAAAGGTTGCATGATAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAAGCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaaa                                                                                                            <  2:780235/151‑1 (MQ=255)
   gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaaa                                                                                                            <  8:1485010/151‑1 (MQ=255)
   gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaaa                                                                                                            <  1:445316/151‑1 (MQ=255)
   gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaaa                                                                                                            <  5:478081/151‑1 (MQ=255)
   gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTaaa                                                                                                            <  5:1663272/151‑1 (MQ=255)
   catctcacaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaa                                                                                                           <  7:594153/143‑2 (MQ=17)
    aaaaaataaGGCGCCAGGAAAACGCGAGCGCTGTATAGACCGCCTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaag                                                                                                          <  7:1029246/150‑3 (MQ=11)
     aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAg                                                                                                           <  1:1030889/150‑1 (MQ=255)
     aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAg                                                                                                           <  3:1230312/150‑1 (MQ=255)
     aaaaaaaaGGTTGAATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaagc                                                                                                         <  4:663134/151‑4 (MQ=255)
      aaaacaaGGTTGCATGAAAACGACATCGTAGTATAGTCCGAATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaagca                                                                                                        <  6:1414229/151‑5 (MQ=14)
      aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCa                                                                                                         <  6:184682/151‑1 (MQ=255)
      aaaaaaaGGCTGCATGAAAACCCCAGCGGAGTATAGTGGGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaagca                                                                                                        <  4:349828/151‑5 (MQ=17)
       aaataaGGGTGAATGAAAAATCGAGCGGAGTATATTGCGCATACCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaagcag                                                                                                       <  3:36972/151‑6 (MQ=11)
       aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTGATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAAGCAATTAAAATTTGATGAAGTAAAGCAg                                                                                                        <  4:909356/151‑1 (MQ=255)
       aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAg                                                                                                        <  8:658582/151‑1 (MQ=255)
        ataaaGGTTGCATGAAAAAGCGATCGGAGTATATTTCGCATCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGt                                                                                                       <  6:147243/149‑1 (MQ=14)
        aaaaaGGTTGCATGAAATCGCGAGCGGAGTATAGTGCGCATCCACGGAGGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAAGCAATTAAAATTTGATGAAGTAAAGCAGt                                                                                                       <  3:1194493/151‑1 (MQ=255)
        aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAATTGAAGAAATaaaa                                                                                                           >  4:596899/1‑145 (MQ=255)
         acaaGGTTGAATGAAAAAGCGAGCGGTGTATACTCCGCCTCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTa                                                                                                      <  8:1019660/149‑1 (MQ=11)
         acaaGGGTGCATGAAATCGCGAGCGGAGTATAGTGCCCATACACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTa                                                                                                      <  3:1350298/149‑1 (MQ=14)
         aaacgggggCATGAAAAAGCTAGCGGAGTATACTGCGCCTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTa                                                                                                      <  8:671395/143‑1 (MQ=11)
         aaaccgctgaaTGAAAACGCGAGCGGTGTATAGTGCGCCTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTa                                                                                                      <  7:942383/141‑1 (MQ=14)
         aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTa                                                                                                      <  2:1324251/151‑1 (MQ=255)
         aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTa                                                                                                      <  4:1129154/151‑1 (MQ=255)
         aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCATCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTa                                                                                                      <  1:1552731/151‑1 (MQ=255)
         aaaaGGTTCAAGGAAAAACCGAGCGCAGTATAGTGCGCATCCCGTGACGCGGGGTGGAGCAGCGTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTa                                                                                                      <  7:1141690/151‑1 (MQ=11)
          aacgcatgaaTGAAAACGCCAGAGGAGTATACTCCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTAc                                                                                                     <  5:745821/142‑1 (MQ=11)
          aaaGGTTGCATGAAAACGCTAGCGGTGTATAGTGGGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTAc                                                                                                     <  7:1140022/151‑1 (MQ=14)
          aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTAc                                                                                                     <  6:1004093/151‑1 (MQ=255)
          aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTAc                                                                                                     <  8:1391663/151‑1 (MQ=255)
           acggcggCATGAAAAAGCGAGCGGAGTATAGAGCGCCTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACg                                                                                                    <  3:334349/145‑1 (MQ=11)
           aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACg                                                                                                    <  3:1313713/151‑1 (MQ=255)
           aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGAAGCAGTAAAAGCAGTACg                                                                                                    >  5:675749/1‑151 (MQ=17)
           aaGGTTGAATGAAAAAGCTTTCGGTGTATATTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACg                                                                                                    <  4:1001918/151‑1 (MQ=11)
            aGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACgg                                                                                                   <  6:1034437/151‑1 (MQ=255)
            aGGTTGCATGAAAACGAGATTGGTGTATATTGCGCTTCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACgg                                                                                                   <  6:1283732/151‑1 (MQ=11)
            aGGTTGAATGAATACTCTAGGGTAGTATAGTGCGATTCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACgg                                                                                                   <  1:677035/151‑1 (MQ=11)
             ggTTGCATGAAATCGCGTGCGTAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGt                                                                                                  <  4:1508555/151‑1 (MQ=14)
             ggTTGCATGAAAACGCGTGGGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGt                                                                                                  <  7:1336663/151‑1 (MQ=14)
             ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGt                                                                                                  <  7:550195/151‑1 (MQ=255)
             ggTTGCAGGAAAACTCCAGCGTAGTATAGTCCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGt                                                                                                  <  4:1304660/151‑1 (MQ=14)
             ggTTGAATGAAAACTCTATCGTAGTATAGTGCGCTTCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGt                                                                                                  <  1:797676/151‑1 (MQ=11)
              gTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTg                                                                                                 <  6:1181454/151‑1 (MQ=255)
              gTTGAATGAATACGCGATGGGAGTATATTCCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTg                                                                                                 <  5:115132/151‑1 (MQ=11)
              gTTGAATGAAAACGCGAGCGTAGTATAGTTCGATTCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTg                                                                                                 <  6:434347/151‑1 (MQ=11)
               gtgaatcaaaaCAAGAGCGCAGTATTTTTCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGa                                                                                                <  2:1392769/144‑1 (MQ=11)
                tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGa                                                                                                <  5:447827/150‑1 (MQ=255)
                ggaaTGAATACGCGTGCGGAGTAGAGTGCGAATCCACGGACGCGGGGTGGAGCAGCCGGGTAGCTCGGCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGa                                                                                                <  4:1435488/147‑1 (MQ=11)
                 gaaTGAAAACTCGATCGGAGTATATTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcg                                                                                              <  2:1437302/149‑1 (MQ=14)
                 gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAAATTTTCAGGATGAAAAGCAGTACGGTGAcg                                                                                              >  3:1428889/1‑151 (MQ=18)
                 gCATGAAAACCCTAGCGTTGTATACTTCGCTTACAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcg                                                                                              <  1:1014255/151‑1 (MQ=11)
                 gCAGGAAAACGCGAGGGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAc                                                                                               <  7:1441672/150‑1 (MQ=14)
                  cATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGAcgc                                                                                             <  3:1125768/151‑1 (MQ=255)
                  cATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGAcgc                                                                                             <  6:717179/151‑1 (MQ=255)
                   aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGAcgcg                                                                                            <  5:259121/151‑1 (MQ=255)
                   aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGAcgcg                                                                                            <  3:755436/151‑1 (MQ=255)
                   aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcgcg                                                                                            <  3:1015069/151‑1 (MQ=255)
                    tGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGTCGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACTGTGAcgcg                                                                                            <  1:1068793/150‑1 (MQ=255)
                      aaaaCTCGAGCTGATTATATTGCTCATTCAAGGACGCGGGGTGGATCATCCTTTTAGCTAGTAGGGCTCATNANCCTNATGNNNTAGGTTTAAATCAGGCACCCTCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCgggg                                                                                         <  2:91935/151‑1 (MQ=12)
                       aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGt                                                                                        <  2:18448/151‑1 (MQ=255)
                       aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGt                                                                                        <  4:892280/151‑1 (MQ=255)
                       aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGGCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGt                                                                                        <  5:1601863/151‑1 (MQ=255)
                        aaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTACAGCAGTACGGTGACGCGGGGTg                                                                                       <  2:848779/151‑1 (MQ=255)
                        aaCGCGAGCGGAGTATAGTGCGCAGCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTg                                                                                       <  4:148788/151‑1 (MQ=255)
                            cgAGCGGAGTATAGTGAGCATCCACGGACGCTTGGTGGAGCAGCCTGGTAGCTCGTCGGTCTCATAACCCGAATGTCGTCTGTTCAAATCCTGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCATTACGGTGACGCGGGGTGGagc                                                                                   <  2:687331/151‑1 (MQ=18)
                             gAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGagca                                                                                  <  8:1483273/151‑1 (MQ=255)
                              cgCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGGCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCGCCGCAACCAATTAAAATTTGATGAAGTAACGCAGTACGGTGACGCGGGGTGGagcag                                                                                 <  1:978133/150‑1 (MQ=255)
                              aggggAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCGTGGTAGCTCGTGGGGCTCATAAGCCGAAGGTCGGCGGTTCATATCCGGCCCCGGCAACCAATTAAAAATTGATGAAGTTATGCAGTACGGTGACGCGGGGTGGagcag                                                                                 <  7:1041606/148‑1 (MQ=18)
                              aGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGTCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGagca                                                                                  <  2:546395/150‑1 (MQ=255)
                              aGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCGATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGagcag                                                                                 <  5:1513867/151‑1 (MQ=255)
                                 ggtgTATAGTGCGCATCCTCGGACGCTGGTTGGAGCAGCGTGTTAGCTGTGCTGTCTCATAAGCCTAAGGTGGTCGGTTGAAATCCTGCCCCCGCAATCAATTAAAATTTGAAGAATTTAAGCATTACCGTGACGCGGGGTGGAGCAGCCt                                                                              <  1:799572/148‑1 (MQ=12)
                                 ggAGTATAGTGAGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTGCAAATCCCGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCt                                                                              <  3:59725/151‑1 (MQ=255)
                                 gcaGTATAGTGATGATCCACCGACGCGGGGATGAGCAGCCTGTTATCGCGTCTGGATCATAAACCGAAGGTCTTCGGTTCAAATAACCAAGCCGGAAACAATTAAAATTAGATGTATTTACGCATTACTCTGACGCGGGGTGGAGCAGCCt                                                                              <  8:1572103/149‑1 (MQ=9)
                                                               tGGAGTAGCCTTGTAGCTCGGCTGTATCATAACCCGAATGTCGTCTGTTGAAAACCTTCCGCCGCAAGGAATTACACATTGTTGAATTCATGCATTCCTGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAgg                                                 <  1:1489941/150‑1 (MQ=2)
                                                                                                           cgGTTCAAATACGGCCGGCGGAACCAATTAAACAGAGATGAAGTATCCCACTCCCGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGaa    <  5:67516/151‑1 (MQ=14)
                                                                                                             gtgcgaaacagcccccgGCCACCAATTACAAAGGGATGAAGAGACGCAGTACGGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGGCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAgc  <  5:550552/139‑2 (MQ=14)
                                                                                                             gTTCCAAACCGCCCCGCGCAACCAATAACAATGCGCTGAAGTAACCGCATACCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGt  <  5:550547/151‑1 (MQ=14)
                                                                                                                             cGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATc                                  >  7:1303961/1‑103 (MQ=255)
                                                                                                                                                       |                                                                                                              
GGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGTAAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGT  >  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: