Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913689,5140AC44.7% 141.7 / 18.0 85G456G (GGT→GGGlntapolipoprotein N‑acyltransferase
Reads supporting (aligned to +/- strand):  ref base A (34/13);  new base C (0/38);  total (34/51)
Fisher's exact test for biased strand distribution p-value = 2.34e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGC  >  NC_000913/689364‑689661
                                                                                                                                                      |                                                                                                                                                   
aCTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCa                                                                                                                                                     <  1:806952/151‑1 (MQ=255)
       aaCACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCg                                                                                                                                              <  2:354532/151‑1 (MQ=255)
       aaCACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCg                                                                                                                                              >  1:1325549/1‑151 (MQ=255)
       aaCACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAATCCGGGGGGCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCg                                                                                                                                              >  2:900529/1‑151 (MQ=255)
         cacaGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCggg                                                                                                                                            >  1:159696/1‑151 (MQ=255)
            agcagcAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGcc                                                                                                                                          >  2:1100051/1‑150 (MQ=255)
            agcagcAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCa                                                                                                                                         <  2:983903/151‑1 (MQ=255)
              cagcaAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAAt                                                                                                                                       <  2:907631/151‑1 (MQ=255)
                 caAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCAc                                                                                                                                    >  1:455001/1‑151 (MQ=255)
                    aaCCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCAccgc                                                                                                                                 <  2:1721338/151‑1 (MQ=255)
                          aCAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAAt                                                                                                                           >  1:62649/1‑151 (MQ=255)
                            aaggCCGTCAGCACCCACAGCGGCCAGTGGCGGGTACGTGCGTATGGTGTGAGTCGGGTGGTCGGCGTCACGTTAGTGTTTAACACCTCGCGGGTGAACGGGGGGATCTTCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATg                                                                                                                          <  2:492636/147‑1 (MQ=255)
                            aaTGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGc                                                                                                                         >  1:1628470/1‑151 (MQ=255)
                               gCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGGATCCGGGGGGCGGCGGCATGATTGAGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCgt                                                                                                                      >  2:759581/1‑151 (MQ=255)
                                ccGTCAGCCCCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCgt                                                                                                                      <  1:1139923/150‑1 (MQ=255)
                                 cGTCAGCCCCCACAGCGGCCAGTGGCGGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGGGGTTAACCCCTCGCGGGTGACCTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAACCACCGCCGTAATGCCgttg                                                                                                                    <  2:1150545/151‑1 (MQ=255)
                                  gTCAGCCCCCACAGCGGCCAGTTGCGGGTACGGGCGTATGGTGTGGGCCCGGTGGTGGGCGCCACGTTAGTGGTTAACCCCTCGGGGGTGAACTGGGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCg                                                                                                                       <  1:528921/147‑1 (MQ=255)
                                  gTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGGGTCCGGTGGTCGGCCTCATGTTTGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCgttgt                                                                                                                   >  2:1700195/1‑151 (MQ=255)
                                  gTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCAAATCACCGCCGTAATGCCgttg                                                                                                                    <  1:1191474/150‑1 (MQ=255)
                                       cccccACAGCGGCCAGTGGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGttgttg                                                                                                                 <  1:1165604/146‑1 (MQ=255)
                                       cACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTACCCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGt                                                                                                               <  1:1441061/150‑1 (MQ=255)
                                       cACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTAGGTCTCGGGGGCGGCCGCACGATAGAGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGtg                                                                                                              >  2:1322885/1‑151 (MQ=255)
                                        aCCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTAAATCCCGTGGGCGGCCGCACCATAGAGAATAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGtgc                                                                                                             >  2:755996/1‑151 (MQ=255)
                                               gCGGCCAGTTGCGGGTAGGTGGGTATGGTGGGGGTCGGGTGGTCGGCTTCACGTTAGTGGTTACCCCCCCGGGGGTGAACTGGGGGATCACCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTgcgc                                                                                                        <  2:145836/149‑1 (MQ=255)
                                                 ggcagtttGCGGGTACGTGGTTATGGTGGGAGTCGGGGGGTGGGCGTCACGTTGGTGTTTAACCCCCCGCGGGTAAACTGCGGGATCATCGCTTGAATCCCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCaa                                                                                                      <  2:1129026/143‑1 (MQ=255)
                                                   cctgTTGCCGGTGCGGGCGTGTGGGGTGAGCCCGGTGGTCGGCGCCCCGTTATTGGTTAACCCCGCGCGGGTAACCGGGGGGATCCTCGCCTGACTCCCCCCCGCGGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAg                                                                                                   <  2:381897/147‑1 (MQ=255)
                                                      gtggCCGGTACGGGCGTATGGTGTGAGTCGGGTGGGCGGCGTCCCGTTAGTGTTTACCCCCCCGGGGGTAAACGGGGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTg                                                                                                 <  2:171538/146‑1 (MQ=255)
                                                                cgggCGTATGTTGTGAGTCCGTTGGTCGGCGTCACTTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCt                                                                                     <  1:1218914/148‑1 (MQ=255)
                                                                 gTGCTTATGGTGTGAGTCCGGTGGTCGGCGTCGCTTTAGTGGTTACCCCCCCGCGGGTGACCGGCGGGTTCACCGCCTGAATCTCCCCCGGCGGCCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCt                                                                                     <  2:639514/150‑1 (MQ=255)
                                                                   gCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGATCACCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCa                                                                                  <  2:360596/151‑1 (MQ=255)
                                                                   gCGTAGGGTGTGAGTCCGGTGGTGGGCGTCACGTTAGTGGTTAACCCCTCGCGGGTGAACTGCGGGACCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCa                                                                                  <  1:1089378/151‑1 (MQ=255)
                                                                     gtttGGTGTGAGTCCGGTGGTCGGCGTCACGTTGGTGTTTAACCCCTCGCGGGGGATCTGCGGGATCACCCCCTAAACCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAgc                                                                                <  1:887682/148‑1 (MQ=255)
                                                                     gTATGGGGTGGGTCGGTTGGTCGGCTGCATGTTAGTGTTTAACCCCTGGCGGGTGAACTGCGGGTTCATCGCCTCAATCTCACCCTGGGGCCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTcc                                                                                   <  1:1127071/148‑1 (MQ=255)
                                                                      tATGGTGTGAGTCGGGTGGTGGGCGTCCCGTTAGTGGTTACCGCCTCGCGGGTGAACGGGGGGACCACCGCCTGAACCCCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAgc                                                                                <  2:26355/150‑1 (MQ=255)
                                                                         ggTGTGAGTCGGGTGGTGGGCGTCACGTTAGTGTTAAACCCCTCGGGGGTGAACGGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgcac                                                                            <  1:494570/151‑1 (MQ=255)
                                                                         ggTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTACCCCCGCGGGGGTAAACGGCGGGATCATCGCCTAAATCTCCCCCTGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGcgca                                                                             <  1:1748472/150‑1 (MQ=255)
                                                                         ggTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAgcg                                                                               <  1:1391032/148‑1 (MQ=255)
                                                                                gTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCTGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTcg                                                                     <  2:180359/151‑1 (MQ=255)
                                                                                 tCCGGTGGTCGGCTTCACGTTAGTGGTTAACCCCTCGCGGGTGAACTGGGGGATCACCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTcg                                                                     <  1:569426/150‑1 (MQ=255)
                                                                                 tCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTcgc                                                                    >  1:42354/1‑151 (MQ=255)
                                                                                    ggtggtGGGCGTCATTTTAGGGTTTAACCCCGCGCGGGTGAAGTGCGGGATCATCGCCTGAATCTCCCCCTGCGGCCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTcgcg                                                                   <  1:1401938/149‑1 (MQ=255)
                                                                                    ggtggtCGGCGTTACGTTAGTGGTTAACCCCTCGGGGGTGAACTGCGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTcgc                                                                    <  1:98622/148‑1 (MQ=255)
                                                                                    ggtggtCGGCGTCACGTTAGTGTTTAACCCCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCCCCCCCTGCGGGCCAACCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTc                                                                      <  2:661069/146‑1 (MQ=255)
                                                                                           ggCGTCCCGTTAGTGGTTACCCCCTCGGGGGTGAACTGCGGGCTCACCGCCTGAATCTCCCCCTGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTgg                                                           <  2:354401/150‑1 (MQ=255)
                                                                                           ggCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTgg                                                           >  2:12926/1‑150 (MQ=255)
                                                                                            gCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGa                                                          >  1:1640384/1‑150 (MQ=255)
                                                                                             cGTCACGTTAGTGGTTACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgt                                                       >  1:1183934/1‑151 (MQ=255)
                                                                                              gTCACGTTAGTGGTTACCCCCCCGCGGGTGAACTGCGGGATCATCGCCTAAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAg                                                        <  1:580730/150‑1 (MQ=255)
                                                                                               tCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgtg                                                      >  2:1724035/1‑151 (MQ=255)
                                                                                               tCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgtg                                                      >  1:334780/1‑151 (MQ=255)
                                                                                                catgTTGGTGGTTACCCCCTGGCGGGTGAACTGCGGGCACTTCGCTTAAATCTCCCCCTGGGGCCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgtg                                                      <  1:453410/147‑1 (MQ=255)
                                                                                                cACGTTAGTGTTTAACCCCTCGCGGGTAAACTGGGGGATCACCGCCTGAATCCCCCCCTGCGGCCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgt                                                       <  2:1349702/149‑1 (MQ=255)
                                                                                                cACGTTAGTGGTTAACCCCTCGGGGGGGAACTGGGGGATCATCGCCTGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgtg                                                      <  2:590984/150‑1 (MQ=255)
                                                                                                cACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgtgt                                                     >  1:1355996/1‑151 (MQ=255)
                                                                                                cACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAgtgt                                                     >  1:814635/1‑151 (MQ=255)
                                                                                                 aCGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGtt                                                    >  2:1186986/1‑151 (MQ=255)
                                                                                                  cGTTAGTGGTTACCCCCGCGCGGGTGAACTGCGGAATCACCCCCGAAACCACACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTg                                                   <  1:1761392/151‑1 (MQ=255)
                                                                                                   tttATTGTTTACCCCCCCGCGGGTAAACTGGGGGATCATCCCCTGAACCCCCCCCTGCGGGCCAACCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGtt                                                    <  2:620419/148‑1 (MQ=255)
                                                                                                   gTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGc                                                  >  1:1270083/1‑151 (MQ=255)
                                                                                                       gtgtTTAACCCCTCGCGGGTGAACTGCGGGATCACCGCCGGAATCTCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCAt                                               <  1:952091/150‑1 (MQ=255)
                                                                                                         ggTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATgg                                             >  2:317268/1‑150 (MQ=255)
                                                                                                          gTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGAc                                           >  1:1163900/1‑151 (MQ=255)
                                                                                                          gTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGAc                                           >  1:776843/1‑151 (MQ=255)
                                                                                                           ttAACCCCCCGGGGGTGACCTGGGGGATCACCGCCTGAATCTCCCCCTGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGAc                                           <  2:744070/150‑1 (MQ=255)
                                                                                                                  cTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAAt                                       >  1:1348992/1‑147 (MQ=255)
                                                                                                                     gcgGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGAttt                                 >  1:672300/1‑150 (MQ=255)
                                                                                                                        ggTGAACTGGGGGATCACCGCCTGAATCCCCCCCTGGGGGCCAACCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCa                             <  2:579988/151‑1 (MQ=255)
                                                                                                                         gTGAACGGGGGGATCATCGCCGGAACCTCCCCCGGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCaa                            <  2:595221/151‑1 (MQ=255)
                                                                                                                             aCTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAAcc                         >  1:1529941/1‑150 (MQ=255)
                                                                                                                              cTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCAc                       >  1:1347799/1‑151 (MQ=255)
                                                                                                                              cTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCAc                       >  2:1256488/1‑151 (MQ=255)
                                                                                                                              cTGCGGGATCATCGCCTGAATCCCCCCCTGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCa                        <  2:575372/150‑1 (MQ=255)
                                                                                                                                gggggATCACCGCCTGAATCCCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCAc                       <  1:549829/147‑1 (MQ=255)
                                                                                                                                  gggTTCACCGCCTGAATCCCCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGc                     <  2:226290/149‑1 (MQ=255)
                                                                                                                                  gggATCTTCGCCTGAATCTCCCCCTGCGGGCCAATCTCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCa                        <  1:1383322/146‑1 (MQ=255)
                                                                                                                                  gggATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCAt                   <  1:123027/151‑1 (MQ=255)
                                                                                                                                      tcatcGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGtt               >  2:17221/1‑151 (MQ=255)
                                                                                                                                       catcCCCGGAACCCCCCCCGGCGGGCCAACCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGtt               <  2:763861/150‑1 (MQ=255)
                                                                                                                                        atcGCCTGATTCTCCCCCGGCGGGCAAATCACCGCCGTAATGCCGTTGTTGGTGCTGTGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg             <  1:1257522/151‑1 (MQ=255)
                                                                                                                                          ccccTAAATCCCCCCCGGCGGGCCACCCGCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGgag           <  2:1723274/149‑1 (MQ=255)
                                                                                                                                           gCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTg              >  1:912500/1‑147 (MQ=255)
                                                                                                                                            ccTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg             >  1:1353009/1‑147 (MQ=255)
                                                                                                                                            ccTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTgg             >  1:1588326/1‑147 (MQ=255)
                                                                                                                                                   ctcCCCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTca    <  2:1499859/149‑1 (MQ=255)
                                                                                                                                                   ctcACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTcagc  >  2:936186/1‑151 (MQ=255)
                                                                                                                                                     cACCCTGCGGGCAAATCCCCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTcag   <  1:1377886/148‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
ACTCATCAACACAGCAGCAAAACCAAACAATGCCGTCAGCACCCACAGCGGCCAGTTGCCGGTACGTGCGTATGGTGTGAGTCCGGTGGTCGGCGTCACGTTAGTGGTTAACACCTCGCGGGTGAACTGCGGGATCATCGCCTGAATCTCACCCTGCGGGCCAATCACCGCCGTAATGCCGTTGTTGGTGCTGCGCAACAGTGGGCGCGCCAGCTCCAGCGCACGCATTCGCGCCATCTGGAAGTGTTGCCATGGACCAATAGATTTACCAAACCACGCATCGTTGGAGATAGTCAGC  >  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 < 33 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: