Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,947,4670AC31.3% 207.0 / 122.6 116intergenic (+4/‑30)metZ/metWtRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base A (5/74);  new base C (3/33);  total (9/107)
Fisher's exact test for biased strand distribution p-value = 7.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 150 total aligned reads displayed.

CAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGG  >  NC_000913/2947318‑2947603
                                                                                                                                                     |                                                                                                                                         
cAAATCCGTGATCGGGTTAAAAAAAAGGTTGCATGAAAACTCGAGCGTAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTat                                                                                                                                          <  8:327965/151‑2 (MQ=255)
cAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTat                                                                                                                                          >  1:707886/1‑150 (MQ=255)
cAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaa                                                                                                                                          <  3:487007/151‑1 (MQ=255)
cAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaa                                                                                                                                          <  1:1435914/151‑1 (MQ=255)
cAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaa                                                                                                                                          <  5:703660/151‑1 (MQ=255)
cAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaa                                                                                                                                          <  5:1084223/151‑1 (MQ=255)
cAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaa                                                                                                                                          <  4:947238/151‑1 (MQ=255)
cAAATCCGTGATCGGGGTAAAAAAAAGGTTGAATGAAAAATCGAGCGGTGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTat                                                                                                                                          <  4:1256085/151‑2 (MQ=255)
 aaaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTaaa                                                                                                                                         <  4:900945/151‑1 (MQ=255)
 aaaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTAAAAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATaaa                                                                                                                                          >  4:251893/1‑150 (MQ=255)
  aaTGCGTGATCGGGGTAAAAAAAAGGGTCCATGAAAACGCGAGCGGAGTATACTGCGCCTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaat                                                                                                                                        <  2:1245966/151‑2 (MQ=18)
  aaTCGGTGATCGGGGTAAAACAAAGGGGGCAGGAAAACGCTATCGTAGTATAGTGCGAATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCGCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattg                                                                                                                                        <  2:801748/151‑4 (MQ=18)
  aaTCCGTGATGGGGTTAAAATTACGGTTGCATGAATACTCTATGGTTGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCGTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCGGCAACCAATTattg                                                                                                                                        <  6:1033728/151‑4 (MQ=17)
  aaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCATGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaat                                                                                                                                        >  4:1184258/1‑150 (MQ=255)
  aaTCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattg                                                                                                                                        <  8:119183/151‑4 (MQ=255)
   atgccggATCTGGATAAAACACACGGTGCAGGAAATCTAGTGCGGAGTATACAGCGCTTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCaaatt                                                                                                                                       <  2:1324552/145‑3 (MQ=11)
   atgcGTGATGGTGGTAAAATAAAGGTTGCATGAATTCGAGAGCGGAGTTTAGTGCGCATCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattga                                                                                                                                       <  8:1348228/148‑5 (MQ=17)
     gggTGATCGGGTTATAATAAAGGTTGAATGAAAAATAGATCGTAGTATATTGCGATTACAGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAtttg                                                                                                                                     <  7:97791/149‑2 (MQ=11)
     gcctGAGCGGGGTAAAAAAAAGCTTGCAGCAAAAAGCGATCGCAGTATAGTCCCCCTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAtttg                                                                                                                                     <  5:393627/148‑2 (MQ=14)
     gcGTGATCGGGTTAAAATAAAGGTTGAAGGAAAACCCGAGCGGAGTATACTCCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAtttg                                                                                                                                     <  4:589775/150‑2 (MQ=17)
     ccgggctCGGGGTAAAAAAAAGGGTGCATCAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAAtttg                                                                                                                                     <  2:678894/145‑2 (MQ=18)
     ccGTGATCGTGTTAAAATTCAGGTTGCATGAATAATCGTGCGGAGTATAGTGCGCTTACTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaa                                                                                                                                      <  6:241227/150‑6 (MQ=14)
     ccGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTatagaac                                                                                                                                     >  8:563403/1‑145 (MQ=255)
     ccGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAttt                                                                                                                                     <  2:586978/151‑1 (MQ=255)
     ccGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAttt                                                                                                                                     <  1:1553717/151‑1 (MQ=255)
     ccGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAAttt                                                                                                                                     <  1:1298787/151‑1 (MQ=255)
      cGTGATCGGTGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGCGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtg                                                                                                                                    <  6:1415814/151‑1 (MQ=255)
      cGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtg                                                                                                                                    <  5:846912/151‑1 (MQ=255)
       gTGATCGGGGTAAAAACCAGGGTGCATGAAATCGCGAGCGTAGTATATTGCGCATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTtgat                                                                                                                                  <  4:1606734/151‑1 (MQ=17)
        tGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacacc                                                                                                                                  >  5:1457219/1‑142 (MQ=255)
        tGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGGTCAAATCCGGCCCCCGCAACCAATAAAAAtttatt                                                                                                                                  >  1:1470949/1‑148 (MQ=255)
         gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccc                                                                                                                                 >  3:1235896/1‑141 (MQ=255)
         gATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTtgatc                                                                                                                                 >  5:731727/1‑150 (MQ=255)
          aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaacaccct                                                                                                                                >  3:107089/1‑140 (MQ=255)
          aTCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaaaaccc                                                                                                                                 >  4:1298835/1‑140 (MQ=255)
            ggggtTAAAAATAAGGTTGCATGAATTCGCTAGCGGAGTATAGTGCGCATCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGt                                                                                                                             <  2:132826/150‑1 (MQ=17)
            cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGGCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGaa                                                                                                                               <  4:291398/150‑1 (MQ=255)
            cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGaa                                                                                                                               <  1:327221/150‑1 (MQ=255)
            cGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGaa                                                                                                                               <  2:807956/150‑1 (MQ=255)
            cGGCGTAAAAAAAAGCGTGAAGGAAAAAGAGAGCGGAGTAGAGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGaa                                                                                                                               <  5:1318837/149‑1 (MQ=14)
             ggggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTatagaacaccctaa                                                                                                                              >  7:823479/1‑137 (MQ=255)
              gggTAAAAACAAGGGTGCATGAAAACCCGAGCGGAGTATAGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aa                                                                                                                          <  6:339675/151‑1 (MQ=17)
              gggTAAAAAAAAGGTTGCATGAAAACTCTTGCGGAGTATTGTGCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aa                                                                                                                          <  6:731621/151‑1 (MQ=17)
              gggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TGAAATTTGGTGAAgga                                                                                                                            >  1:340283/1‑148 (MQ=255)
               gttAAAAATCAGGTGGCATCAAATCTAGAGTGGTGAATATTGCCAATACAGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aa                                                                                                                          <  3:995029/148‑1 (MQ=11)
               ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aa                                                                                                                          <  7:1664659/151‑1 (MQ=255)
               ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGGCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aa                                                                                                                          <  8:1019551/151‑1 (MQ=255)
                  aaaagaaaGGGTGAATGAATAATAGTGCGTAGTATAGTGCGACTCCCGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aaaagc                                                                                                                      <  2:1589222/151‑4 (MQ=11)
                  aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAg                                                                                                                        <  3:1623487/150‑1 (MQ=255)
                  aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAg                                                                                                                        <  3:763068/150‑1 (MQ=255)
                  aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAg                                                                                                                        <  8:921274/150‑1 (MQ=255)
                  aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCAGCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAg                                                                                                                        <  4:731291/150‑1 (MQ=255)
                   aaataaaGGGTGAATGAAATCTATAGCGGAGTATAGTGCGATTACAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aaaagca                                                                                                                     <  4:1610503/151‑5 (MQ=11)
                   aaaaacaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCa                                                                                                                      <  2:213100/151‑1 (MQ=255)
                   aaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCa                                                                                                                      <  5:385160/151‑1 (MQ=255)
                    aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGaa                                                                                                                               >  7:593954/1‑142 (MQ=255)
                    aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGaa                                                                                                                               <  2:592622/142‑1 (MQ=255)
                    aaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAg                                                                                                                     <  8:1281911/151‑1 (MQ=255)
                     aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTTATGAAGT‑aaaagctcg                                                                                                                   >  1:917750/1‑145 (MQ=18)
                     aaaaaGGGTGCATGAAAACGCGTGCGGAGTATATTGCGCATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aaa                                                                                                                         <  6:917374/145‑1 (MQ=17)
                     aaaaaGGCTGCATGAAAACACGAGCGGAGTATACTCCGCATCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGt                                                                                                                    <  5:1537649/151‑1 (MQ=14)
                      aaaaGGTTGCATGAAAACGCGATCGGAGTATACTCCGCATCCATGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTa                                                                                                                   <  5:1183619/151‑1 (MQ=14)
                      aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTa                                                                                                                   <  4:863288/151‑1 (MQ=255)
                       ccagcgtgcaccaaaaCCCGAGCGGAGAATAGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTAc                                                                                                                  <  7:265285/139‑1 (MQ=11)
                       aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTAc                                                                                                                  <  1:1176940/151‑1 (MQ=255)
                       aaaGGTTGAATGAATTAGCGATAGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTAc                                                                                                                  <  3:854223/151‑1 (MQ=14)
                         aGGTTCCATGAAAAAGCGTTCGGAGTATACTGCGATTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAGTCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACgg                                                                                                                <  6:360356/151‑1 (MQ=11)
                          gggtGCATGAAAAATCGAGCGGAGTATAGTCCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aa                                                                                                                          <  2:871110/136‑1 (MQ=14)
                          ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGt                                                                                                               <  4:230916/151‑1 (MQ=255)
                           gTTGAATGAAAACTATAGCGTAGTATATTGCGCATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTg                                                                                                              <  6:1630803/151‑1 (MQ=11)
                           gTTGAATGAAAACGATTGTGTAGTATATTGCGAATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTg                                                                                                              <  6:1159979/151‑1 (MQ=11)
                             tgaaTGAAAACTCGTGCGGTGTATAGACCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGa                                                                                                             <  1:1053676/147‑1 (MQ=11)
                             tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGa                                                                                                             <  3:838444/150‑1 (MQ=255)
                             ggCATGAAAACCCGAGCGGAGTATACTGCGCCTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGAc                                                                                                            <  6:109527/150‑1 (MQ=11)
                              gaaTGAAAAAGCGAGGGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGAc                                                                                                            <  1:30666/148‑1 (MQ=14)
                              gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGAc                                                                                                            <  1:212441/150‑1 (MQ=255)
                              gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGAc                                                                                                            <  6:201517/150‑1 (MQ=255)
                                aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGAcgcg                                                                                                         <  3:533556/151‑1 (MQ=255)
                                  taaatggcgTGCGGAGTATAGTGCGCATGGACGGACTCGGGGTGGAGCAGCCTGGTAGCTCGTGGGGCTCATAACGCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCgg                                                                                                        <  6:1428451/144‑1 (MQ=17)
                                  gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCggg                                                                                                       <  7:538414/151‑1 (MQ=255)
                                  gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCggg                                                                                                       <  6:1271039/151‑1 (MQ=255)
                                  gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCgg                                                                                                        <  7:1608326/150‑1 (MQ=255)
                                  gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCgg                                                                                                        <  4:191133/150‑1 (MQ=255)
                                  caaatCGCTAGCGGAGAATTGTTCGCATCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGCCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCggg                                                                                                       <  7:1042783/150‑1 (MQ=11)
                                      aCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTgg                                                                                                   <  2:1232929/151‑1 (MQ=255)
                                          gAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTtgat                                                                                                                                  >  8:781187/1‑116 (MQ=255)
                                           aGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGGCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGagcag                                                                                              <  4:207716/151‑1 (MQ=255)
                                               gAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTg                                                                                          <  5:222360/151‑1 (MQ=255)
                                               gAGTATAGTGCGCATCCACGGACGCGGGGTAGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTg                                                                                          <  1:622222/151‑1 (MQ=255)
                                               ccgaactacccccctCCCCGGACGCGGGGTGGAGCAGCCTTGTAGCTCGTCTTCCTCATAACCCGAAGGTCTTCCGTTCAAATCCGTCCCCCGCAACAAATGAATATTTTTGTATTT‑TAAGGACTACTGTGACGCGGGGTGGAGCAGCCTg                                                                                          <  8:828810/137‑1 (MQ=2)
                                                 gTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGt                                                                                        <  5:947337/151‑1 (MQ=255)
                                                   ataggcagcATCCACAGACCCGTGTGGGAGAAGGCTGTTCTCTCGGCTGGATCATAAGCCGAAGGGCGTAGGGTCTTAACCGGCCCCCGCAACGAATTAAACTTCGCTGAAGT‑ATCCCACTACAGGGACGCGGGGTGGAGCAGCCTGGTa                                                                                       <  5:859627/143‑1 (MQ=9)
                                                     aggcaccATCCACCGACCAGGGGGGGAGCAGCCTGGTCGCTCGTCGCTCTCACAACCCCAATGCCGGCGGTTCAAAACCCCACCCCGCAACCAATTACCCAGCGCAGAAGT‑CACGCCCTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCt                                                                                    <  2:1063613/145‑1 (MQ=9)
                                                      gTGGGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattgaac                                                                                                                                     <  7:789673/102‑7 (MQ=11)
                                                          gcATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTattga                                                                                                                                       >  3:81521/1‑92 (MQ=25)
                                                              ccACGGACGAGGGGGGGAGCAGCCTGGTAGCTCGGCGGGCTCATAACCCGAAGGGCGGCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AACGCAGTACCGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCt                                                                           <  4:359784/151‑1 (MQ=14)
                                                                  ccaCGCTGTGTGGAGGAGCCTGTTCGGTCGTCGGGAAAATATCCCGAAGGCCGGCTGTTCAAATAAGGCCCCCTCAACCAATAAACTATTGATGAAGT‑AAAGGATTCCGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATa                                                                       <  5:1624018/149‑1 (MQ=2)
                                                                           gTGGAGCAGGCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGGCGGTTCAAATCCCGCCCCCGCAACCAATTAAAATTTGCTGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAg                                                               <  5:1204160/150‑1 (MQ=255)
                                                                               cgcagcCGCGTAGCTCGTCGGGCTAATAAGCCGAAGGGCGTCGGCGCCCATACCCCCCCCGCAACGAATTACACATTGCTGAATA‑AACGCACTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGgtcgt                                                          <  2:332413/150‑1 (MQ=2)
                                                                                                                                       ccccGCAACCAATTAACTAGTGTTGAAGT‑ATAGCATTCCGGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACgg  <  5:186524/151‑1 (MQ=17)
                                                                                                                                                     |                                                                                                                                         
CAAATCCGTGATCGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGG  >  NC_000913/2947318‑2947603

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: ‑