Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913689,5140AC47.0% 141.9 / 23.8 100G456G (GGT→GGG) lntapolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (37/16);  new base C (0/47);  total (37/63)
Fisher's exact test for biased strand distribution p-value = 4.61e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 103 total aligned reads displayed.

ACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGC  >  NC_000913/689364‑689661
                                                                                                                                                         |                                                                                                                                                   
aCTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCa                                                                                                                                                     <  1:1122657/151‑1 (MQ=255)
aCTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCa                                                                                                                                                     <  1:179375/151‑1 (MQ=255)
 cTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCa                                                                                                                                                     <  2:1156052/150‑1 (MQ=255)
  tcatcaACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCAcc                                                                                                                                                   >  2:347487/1‑151 (MQ=255)
      caACACTGCAGCAAAACCAAACACTGCCGTCAGCACCCACAGCGGCCAGTTGCGGGTACGTGAGTATGATGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTg                                                                                                                                                <  1:797314/150‑1 (MQ=255)
      caACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGc                                                                                                                                               >  2:719803/1‑151 (MQ=255)
       aaCACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTgag                                                                                                                                              <  2:780489/151‑3 (MQ=255)
       aaCACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCg                                                                                                                                              <  2:242230/151‑1 (MQ=255)
          acaGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGc                                                                                                                                           >  2:1411514/1‑151 (MQ=255)
           cagcagCAAAACCAAACATTGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGcc                                                                                                                                          <  1:1790186/151‑1 (MQ=255)
            agcagcAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCa                                                                                                                                         <  1:440597/151‑1 (MQ=255)
             gcagcaAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCaa                                                                                                                                        <  2:382022/151‑1 (MQ=255)
               agcaAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTAGATCTCGGGGGCCGCCG‑CATGATTGAAGATAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAAt                                                                                                                                       >  1:1248895/1‑150 (MQ=255)
               agcaAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATc                                                                                                                                      <  1:1441532/151‑1 (MQ=255)
                  aaaaCCAAACAATGCCGTCAGCCCCCACAGCGGCCAGTTGCGGGTACGGGCGTAGGGTGT‑GAGTCCGGTGGTCGGCGCCACGTTAGTGTTTAACACCTCGC‑GGGTGAACTGCGGGATC‑TTCGCCTGAATCTCCCCCTGCGGGCCAATCAcc                                                                                                                                   <  1:1850263/151‑1 (MQ=255)
                   aaaCCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCAccg                                                                                                                                  <  1:1133793/151‑1 (MQ=255)
                       cAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGt                                                                                                                              >  1:423908/1‑151 (MQ=255)
                        aaaCAATACCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCCGGGGGCCGCGTCATGATTGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGt                                                                                                                              >  1:372761/1‑150 (MQ=255)
                             aTGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGGGGGCGGCGTCACGATAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGcc                                                                                                                        >  1:437702/1‑151 (MQ=255)
                                  gTCAGCCCCCACGGCGGCCAGTTGCCGGTACGTGCGTAGGGTGT‑GAGTCGGGTGGTCGGCGTCCCGTTAGTGGTTAACCCCTCGC‑GGGTGAACGGGGGGATC‑ATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGcc                                                                                                                        <  1:1326069/146‑1 (MQ=255)
                                  gTCAGCCCCCACAGCGGCCAGTTGCGGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTAAACACCTCGC‑GGGTGAACGGCGGGATC‑ACCGCCTAAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCgttgt                                                                                                                   <  1:674208/151‑1 (MQ=255)
                                  gTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCgttg                                                                                                                    >  2:512824/1‑150 (MQ=255)
                                   tCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCgttgtt                                                                                                                  >  2:596833/1‑151 (MQ=255)
                                     aGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCCCCCTGCGGGCCAACCACCGCCGTAATGCCGttgttg                                                                                                                 <  2:956001/150‑1 (MQ=255)
                                       cccccACAGCGGCCAGTTGCGGGTACGTGCGTAGGGTGT‑GATTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCgttgt                                                                                                                   <  1:512094/144‑1 (MQ=255)
                                       cccccACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGt                                                                                                               <  2:161566/148‑1 (MQ=255)
                                       cACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑AAATCCGGGGGTCCGCGGCACGATAGAGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGtg                                                                                                              >  1:953034/1‑151 (MQ=255)
                                         cccTCAGCGGCCAGTGGCGGGTACGTTCGTATGGTGT‑GAGTCCGTGGGTGGGGCTCTGTTTAGTGGTTACCCCCCCGC‑GGGAAAACTGC‑GGATCAACCCCCTTATTCCCCCCCGGCGGGCCAATCGCCGCCGTAATGCCGTTGTTGGtgc                                                                                                             <  1:1577788/150‑1 (MQ=255)
                                         cccACAGCGGCCGTTTGCCGGTCGGTGCGTAGGGTGT‑GGGTCCGGTGGCCGGCGCCCCGTTAGGGTTAAACCCCTCGC‑GGGTGACCGGCGGGATC‑ACCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGtgc                                                                                                             <  1:1770940/150‑1 (MQ=255)
                                         cccACAGCGGCCATTTGCCGTTACGTGCGTAGGGTGT‑GAGCCCGGTGGTCGGCGTCACGTTAGTGGTTACCCCCTCGC‑GGGTGACCGGCGGGATC‑TTCGCCTGAATCTCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGtgct                                                                                                            <  1:1107428/151‑1 (MQ=255)
                                          ccACAGCGGCGAGTGGCGGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGC‑GGGTAAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGtgctg                                                                                                           <  2:326268/151‑1 (MQ=255)
                                           cacaGCGGCCATTGGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGC‑GGGTGAACGGCGGGATC‑ATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGtgctgc                                                                                                          <  2:374199/151‑1 (MQ=255)
                                             caGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTgcgc                                                                                                        >  2:144757/1‑151 (MQ=255)
                                                 ggCCAGTTGCCGGTAGGTGCGTATGGTGT‑GAGCCCGGTGGTCGGCTTCACGTTAGTGGTTAACCCCTCGG‑GGGTAAACTGGGGGATC‑ACCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAAc                                                                                                     <  1:354453/150‑1 (MQ=255)
                                                   ccAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGC‑GGGTGAACTGGGGGATC‑ATCCCCTAAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGt                                                                                                  <  1:1957600/151‑1 (MQ=255)
                                                   ccAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGGCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGt                                                                                                  >  1:1283387/1‑151 (MQ=255)
                                                   ccAGTGGCCGGTACGTGCGTAGGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGGGGTTAACACCTCGG‑GGGTGAACGGCGGGATC‑ATCGCCTGAATCCCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGt                                                                                                  <  1:40412/151‑1 (MQ=255)
                                                    cAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTg                                                                                                 <  2:521084/151‑1 (MQ=255)
                                                          ccGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcgcg                                                                                           >  2:1512247/1‑151 (MQ=255)
                                                          ccGGTACGTGCGTATGGTGT‑GAGTCCAGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGC‑GGGTGAACTGGGGGATC‑ATCGCCTGAATCTCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcgcg                                                                                           <  1:442433/151‑1 (MQ=255)
                                                            ggTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcgcgc                                                                                          <  2:922014/150‑1 (MQ=255)
                                                            ggTACGTGCGTATGGTGT‑AAGTCTCGGGGGCGGCGGCACGATAGTAGATAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGcc                                                                                         >  1:647064/1‑151 (MQ=255)
                                                             gTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGgcgcgc                                                                                          <  1:1464984/149‑1 (MQ=255)
                                                                cGTGCTTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCCCCCTGGGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCa                                                                                        <  2:863303/148‑1 (MQ=255)
                                                                cGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCt                                                                                     >  2:1258139/1‑151 (MQ=255)
                                                                   gCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTcc                                                                                   <  1:1089792/150‑1 (MQ=255)
                                                                      tATGGTGT‑GAGTCGGGTGGTCGGCGTCACGTTAGTGGTTACCCCCCCGC‑GGGTGACCTGGGGGATC‑ATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAgc                                                                                <  1:635271/150‑1 (MQ=255)
                                                                      tATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCCCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCa                                                                                  <  1:1549764/148‑1 (MQ=255)
                                                                         ggTGT‑GAGCCGGGTGTTCGGCGTCCCTTTAGTGGTTAACACCTCGG‑GGGGAACCTGCAGGATC‑ATCGCAAGAATCTCCCCCGGCGGGCAAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAgcg                                                                               <  2:1248214/148‑1 (MQ=255)
                                                                          gtgt‑GAGTCGGGTGGTGGGCGTCCCTTTAGTGGTTAACACCGCGC‑GGGTGAACTGGGGGATC‑ATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgcacg                                                                           <  1:1021705/151‑1 (MQ=255)
                                                                            gt‑gAGTCCGGTGGTGGGCGCCACGTTAGTGTTTAACACCCCGC‑GGGTGAACTGCGGGTTC‑ATCGCTTAAATCTACCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgcacgc                                                                          <  2:203210/150‑1 (MQ=255)
                                                                               gAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCAtt                                                                       >  1:671942/1‑151 (MQ=255)
                                                                               gAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCAtt                                                                       >  2:1118747/1‑151 (MQ=255)
                                                                                 gTCCGGTGGTCGGGGTCACGTTAGGGTTTAACACCCCGG‑GGGTAAACGGCGGGATC‑ACCGCCTGAACCCCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTc                                                                      <  1:1093814/150‑1 (MQ=255)
                                                                                   ccGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTcgcg                                                                   >  1:358535/1‑151 (MQ=255)
                                                                                     gggggccggggccACGTTAGGGGTAACCCCCTCGC‑GGGGGAACTGGGGGATC‑ACCGCCTGAATCCCCCCCTGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCAtt                                                                       <  1:1222640/133‑1 (MQ=255)
                                                                                        ggtCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATc                                                              >  2:1664156/1‑151 (MQ=255)
                                                                                         gtCGGGGTCACGTTAGTGGTTAACCCCTCGC‑GGGTAAACTGGGGGATC‑ACCCCCTGAATCCCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCt                                                             <  1:916394/151‑1 (MQ=255)
                                                                                         gtCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTAACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGGGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCCCGCATTCGCGCCATCt                                                             >  1:1665356/1‑151 (MQ=255)
                                                                                           cGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTgg                                                           >  1:243533/1‑151 (MQ=255)
                                                                                               gTCACGTTAGTGGTTAACCCCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAg                                                        <  2:813252/150‑1 (MQ=255)
                                                                                                 cccGTTAGTGGTTACCCCCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAg                                                        <  1:752215/146‑1 (MQ=255)
                                                                                                 cACGTTAGTGTTTAACCCCTCGC‑GGGTGAACTGCGGGACC‑ATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgtg                                                      <  1:616758/150‑1 (MQ=255)
                                                                                                  aCGTTAGTGGTTAACCCCTCGG‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCCCCCTGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgt                                                       <  2:1691451/148‑1 (MQ=255)
                                                                                                   cGTTAGTGTTTACCCCCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCCCCCTGCGGCCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGtt                                                    <  1:66305/150‑1 (MQ=255)
                                                                                                     ttAGTGGTTAACGCC‑CGCGGGGTGAACTGCGGGATC‑TTCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGc                                                  <  2:1094398/150‑1 (MQ=255)
                                                                                                     ttAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCCGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTg                                                   <  2:940509/149‑1 (MQ=255)
                                                                                                      tAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGcc                                                 >  2:166165/1‑150 (MQ=255)
                                                                                                       aGTGGTTAACCCCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCa                                                <  2:1396089/150‑1 (MQ=255)
                                                                                                        gTGGTTAACCCCTCGC‑GGGTGAATTGCGGGATC‑ATCGCCTGAATCTCCCCCGGGGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATg                                              <  2:74216/151‑1 (MQ=255)
                                                                                                        gTGGTTAACCCCTCGC‑GGGGGAACTGCGGGATC‑ATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCAt                                               <  2:1155535/150‑1 (MQ=255)
                                                                                                         ttgtaaaCCCCTCGC‑GGGTGATCTGGGGTATC‑ATCGCCTAAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCa                                                <  2:1277539/143‑1 (MQ=255)
                                                                                                           gTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATggtc                                           >  1:1708155/1‑149 (MQ=255)
                                                                                                               accccTCGC‑GGGTGAACTGCGGGTTC‑TTCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCa                                         <  2:794198/146‑1 (MQ=255)
                                                                                                                 acCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAAt                                       >  1:1323045/1‑149 (MQ=255)
                                                                                                                  ccTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAAt                                       >  1:1794400/1‑148 (MQ=255)
                                                                                                                  ccTCGC‑GGGGGAACTGCGGGATC‑ACCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGa                                    <  2:1874358/151‑1 (MQ=255)
                                                                                                                   cccgc‑gGGTGAACTGCGGGATC‑ATCGCCTGAACCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAt                                   <  2:83793/149‑1 (MQ=255)
                                                                                                                      gc‑gGGTGAACTGCGGGATC‑ATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAttt                                 <  2:1395817/150‑1 (MQ=255)
                                                                                                                      gc‑gGGTGAACTGCGGGATC‑ATCGCCTGAATCTCCCCCTGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAttt                                 <  1:1273276/150‑1 (MQ=255)
                                                                                                                      gc‑gGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAtt                                  >  1:1350369/1‑149 (MQ=255)
                                                                                                                         ggggGAACTGCGGGTTC‑ATCGCCTGAACCCCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTAc                               <  1:1497485/150‑1 (MQ=255)
                                                                                                                          ggTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTAcc                              >  1:366956/1‑150 (MQ=255)
                                                                                                                                  gCGGGATC‑ATCGCCTGAACCTCCCCCTGGGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACg                      <  1:801985/150‑1 (MQ=255)
                                                                                                                                    gggTTC‑ATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCAt                   <  1:1079490/151‑1 (MQ=255)
                                                                                                                                    gggATC‑ATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGc                     <  2:21427/149‑1 (MQ=255)
                                                                                                                                    gggATC‑ACCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCa                    <  1:1218649/150‑1 (MQ=255)
                                                                                                                                         c‑atcGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTg              >  1:1253751/1‑151 (MQ=255)
                                                                                                                                         c‑atcGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTg              >  2:1295346/1‑151 (MQ=255)
                                                                                                                                           atcGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg             >  1:1146501/1‑151 (MQ=255)
                                                                                                                                             cGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg             >  2:1422658/1‑149 (MQ=255)
                                                                                                                                              gCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGga            <  1:1040230/149‑1 (MQ=255)
                                                                                                                                               ccTGAATCTCCCCCTGGGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGAt         <  2:617207/151‑1 (MQ=255)
                                                                                                                                               ccTGAATCCCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGAt         <  1:1443510/151‑1 (MQ=255)
                                                                                                                                                  gAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGt      >  1:473427/1‑151 (MQ=255)
                                                                                                                                                  gAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGt      >  2:1580231/1‑151 (MQ=255)
                                                                                                                                                  gAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGt      >  1:1072861/1‑151 (MQ=255)
                                                                                                                                                    aTCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTc     >  1:583219/1‑150 (MQ=255)
                                                                                                                                                      ctcCCCCTGCGGGCCAATCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTcag   <  1:1125142/150‑1 (MQ=255)
                                                                                                                                                       tcACCCTGCGGCCCCATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCCAACCACGCATCGTTGGAGATAGTcagc  >  1:1855587/1‑150 (MQ=255)
                                                                                                                                                         |                                                                                                                                                   
ACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGT‑GAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGC‑GGGTGAACTGCGGGATC‑ATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGC  >  NC_000913/689364‑689661

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