Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,165,758 C→G 5.4% intergenic (+73/‑327) ycfP → / → ndh putative UPF0227 family esterase/respiratory NADH dehydrogenase 2/cupric reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,165,7580CG 5.4% 297.0 / 6.9 93intergenic (+73/‑327)ycfP/ndhputative UPF0227 family esterase/respiratory NADH dehydrogenase 2/cupric reductase
Reads supporting (aligned to +/- strand):  ref base C (45/43);  new base G (2/3);  total (47/46)
Fisher's exact test for biased strand distribution p-value = 6.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

GAAGAGCAGACGCACAAATTCAAGAATATCTCCCCGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAACTCTGTTTTTTGATCTCACCCGGTAA  >  NC_000913/1165611‑1165907
                                                                                                                                                   |                                                                                                                                                     
gAAGAGCAGACGCACAAATTCAAGAATATCTCCCCGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCAc                                                                                                                                                     >  2:309799/1‑150 (MQ=255)
 aaGAGCAGACGCACAAATTCAAGAATATCTCCCCGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCCTCAGGTAACCATGCCGGGCCCAGATGCGGTCCCAATTATGTGGGTCGGATAAGGCGTCCAcg                                                                                                                                                    >  1:634698/1‑150 (MQ=255)
 aaGAGCAGACGCACAAATTCAAGAATATCTCCCCGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCAcg                                                                                                                                                    <  2:1308771/150‑1 (MQ=255)
  agagCAGACGCACAAATTCAAGAATATCTCCCCGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCAcgc                                                                                                                                                   >  1:931491/1‑150 (MQ=255)
  agagCAGACGCACAAATTCAAGAATATCTCCCCGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATAAGGTAACCTTGCCGGTACCAGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGGACAcgc                                                                                                                                                   >  1:828598/1‑150 (MQ=255)
     gCAGACGCACAAATTCAAGAATATCTCCCCGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTTACCTTTCCGGTACCTGATTCGCGCCGATTTCTGTGGGTCGGATAAGGCGGCCAcgccgc                                                                                                                                                >  1:497720/1‑150 (MQ=255)
        gACGCACAAATTCAAGAATATCTCCCCGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATGCGGCGTGGAcgccttatc                                                                                                                                             <  2:1108459/150‑6 (MQ=255)
          cGCACAAATTCAAGAATATCTCCCCGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATTCGCTCCGAATTCTGAGGGTCGGATAAGGCGTCCACGCCGCATCCg                                                                                                                                           >  2:1055132/1‑150 (MQ=255)
            cacaAATTCAAGAATATCTCCCCGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGAc                                                                                                                                         >  1:538969/1‑150 (MQ=255)
            cacaAATTCAAGAATATCTCCCCGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGAc                                                                                                                                         <  1:1138105/150‑1 (MQ=255)
            cacaAATTCAAGAATATCTCCCCGCATATACAGCGCAATAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCAGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATGCGGCGTGGACGCCTTATCCGAc                                                                                                                                         <  2:260226/150‑1 (MQ=255)
                  ttCAAGAATATCTCCCCGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGa                                                                                                                                   <  1:1304708/150‑1 (MQ=255)
                   tCAAGAATATCTCCCCGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAg                                                                                                                                  <  2:341348/150‑1 (MQ=255)
                   tCAAGAATATCTCCCCGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCAGCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAg                                                                                                                                  <  2:277721/150‑1 (MQ=255)
                         atatCTCCCCGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCaa                                                                                                                            >  2:1127572/1‑150 (MQ=255)
                         atatCTCCCCGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTAGGATAAGGCGTACACGCCGCATCCGACAGTCGAGCATCaa                                                                                                                            >  2:1051748/1‑150 (MQ=255)
                           aactcCCCGCATTTACAGAGCATTAAAGCGTTCAAAACCCTCGGCTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCCCACCGAATTCTGTGGGTCGGATAAGGCCACCACGCCGCATCCGACAGTCGAGCAACACTg                                                                                                                          >  2:8627/3‑150 (MQ=255)
                             ctcCCCGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATAAGGTAACCATGCCGGTAACGGATTCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGCGCATCAATGcc                                                                                                                        >  2:526995/1‑150 (MQ=255)
                              tcCCCGCATTTACAGCGCATTAAAGCGTTCAATACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTAACTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCt                                                                                                                       <  1:1539107/150‑1 (MQ=255)
                              accccGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCt                                                                                                                       <  1:1370530/149‑1 (MQ=255)
                               ccccGCATTTACAGCTCATTTAAGAGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCAGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCatgcattg                                                                                                                      <  2:1519946/150‑8 (MQ=255)
                               ccccGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTg                                                                                                                      >  1:1201855/1‑150 (MQ=255)
                                 ccGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATGCGGCGTGGACGCCTTATCCGAc                                                                                                                                         <  1:382887/129‑1 (MQ=255)
                                 ccGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATGCGGCGTGGACGCCTTATCCGAc                                                                                                                                         >  2:382887/1‑129 (MQ=255)
                                     atatACAGCCCACTACAGCGTTCAAAACCCACGGGTAAATCCCCTCGTCGCATCAGGCAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATgcgc                                                                                                                <  2:715022/147‑1 (MQ=255)
                                      tttACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGct                                                                                                               <  2:1115752/150‑1 (MQ=255)
                                              gcATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCa                                                                                                       <  2:785374/150‑1 (MQ=255)
                                                 ttAAAGCGTTCAAAACCCTCGGGTAAATGGACTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGAAAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGc                                                                                                    <  2:154556/150‑1 (MQ=255)
                                                      gCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTAcc                                                                                               <  2:931491/150‑1 (MQ=255)
                                                      gCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATAAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTAGAGCATCACTGCCTGATGCGCTTCTTATCAGGCCTAcc                                                                                               >  2:739272/1‑150 (MQ=255)
                                                      gCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCAGCAGGTAACCTTGCCGGTACCTGACGCGCTCCGAATTCTGTGGGGCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTAcc                                                                                               <  1:444863/150‑1 (MQ=255)
                                                      gCGTTCAAAACCCACGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTAcc                                                                                               <  1:178891/150‑1 (MQ=255)
                                                           cAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGCTGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGCTGCGCTTCTTATCAGGCCGACCGAACg                                                                                          >  2:1155147/1‑150 (MQ=255)
                                                           cAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATTCGCTTCTTATCAGGCCTACCGAACg                                                                                          >  2:1114088/1‑150 (MQ=255)
                                                               aCCCTCGTTTAAATGCCCTCGTCTCTTCTGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTTGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCt                                                                                      <  2:1703301/150‑1 (MQ=255)
                                                                     gggTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATac                                                                                >  1:276075/1‑150 (MQ=255)
                                                                      ggTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGGTCCGAATTCTGTGGGTCGGATTAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATacc                                                                               >  1:649561/1‑149 (MQ=255)
                                                                      ggTAAATGCCCTCGGCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATaca                                                                               <  1:154772/150‑1 (MQ=255)
                                                                        tAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACAcc                                                                             <  1:505762/150‑1 (MQ=255)
                                                                         aaaTGCCCTCGTCGCATCAGGTAACCTTGCCGGTTCCAGGTGCGCGCCGAATTCTGTGGGTCGGATAAGGCGTCCCCGCCCCATCCGACAGTCGAGCATCAATTCCTGATGCGCTTATTATCAGGCCTACCGCACGCCCAGCATAAAccc                                                                            >  2:594301/1‑150 (MQ=255)
                                                                         aaaTGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATTCGCATCTTATCAGGCCTACCGAACGCCCTGCATagcccc                                                                            >  2:48013/1‑145 (MQ=255)
                                                                          aaTGCCCTCGTCGCTTCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCTTCCGACAGTCGAGCATCAATGCCTGATGCGATTCTTATCAGGCCTACCGAACGCCCTGCATACAcccc                                                                           <  2:528934/150‑1 (MQ=255)
                                                                           aTGCCCTCGTCGCATCAGGTAACGTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCt                                                                          <  2:1192168/150‑1 (MQ=255)
                                                                           aTGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCt                                                                          >  1:915815/1‑150 (MQ=255)
                                                                            tGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCAGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAGCGCCCTGCATACCCCCCTc                                                                         >  1:1556272/1‑150 (MQ=255)
                                                                             gCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACATTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCa                                                                        >  2:220542/1‑150 (MQ=255)
                                                                              cccTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACCCCCCTccc                                                                       >  2:670378/1‑148 (MQ=255)
                                                                              cccTCGTCGCATCAGGTAACCATGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCCTACGGCAGTCGAGCATCAATGCCTGATGCGCTTCGTATCAGGCATACCGCACGCCCTGCCTACACCGCTCAc                                                                       >  1:769822/1‑150 (MQ=255)
                                                                                  cgtcgCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCta                                                                   <  2:538969/150‑1 (MQ=255)
                                                                                  cgtcgCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCATGCATACACCCCTCACTCta                                                                   >  1:1398415/1‑150 (MQ=255)
                                                                                    tcgCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCtata                                                                 <  1:207486/150‑1 (MQ=255)
                                                                                     cgCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCCACAGTCGAGCATCAATGCCGGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACCCCCCGCCCGCTatat                                                                >  1:717090/1‑150 (MQ=255)
                                                                                        aTCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGTGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCAc                                                             <  2:1595891/150‑1 (MQ=255)
                                                                                        aTCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCAc                                                             <  2:430651/150‑1 (MQ=255)
                                                                                        aTCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCAc                                                             >  2:1668886/1‑150 (MQ=255)
                                                                                         tCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATACGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCCCTCGATATcgcg                                                            >  1:674937/1‑147 (MQ=255)
                                                                                             gTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGGCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTca                                                        <  2:41772/150‑1 (MQ=255)
                                                                                                 ccTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACaaa                                                    >  2:383708/1‑150 (MQ=255)
                                                                                                  cTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCCCAAAt                                                   >  1:837302/1‑150 (MQ=255)
                                                                                                   ttGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAAtt                                                  <  1:1676112/150‑1 (MQ=255)
                                                                                                     gCCGGTTCCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCg                                                <  2:766562/150‑1 (MQ=255)
                                                                                                     gCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCg                                                <  2:828598/150‑1 (MQ=255)
                                                                                                     gCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCg                                                <  1:1342950/150‑1 (MQ=255)
                                                                                                         gTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCa                                            <  1:309799/150‑1 (MQ=255)
                                                                                                            ccTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAaat                                         <  1:732608/150‑1 (MQ=255)
                                                                                                              tGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAAACCCTGCAGCGAGTCTTATCAGGCCTACCGAACGCCCAGCATACACCCCGCACTCTAAATCACTCACACAAATTCGCTCAaataa                                       >  2:301655/1‑150 (MQ=255)
                                                                                                               gATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCAGACACCCCTCACTCTATATGACTCTCACAAATTCGCTCAaatact                                      >  2:705492/1‑148 (MQ=255)
                                                                                                                aTGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAataata                                     <  2:1201855/150‑1 (MQ=255)
                                                                                                                   cgcTCCGAATTCTGTGGGTCGTATAAGGCGGCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATaataaac                                  <  2:1708995/150‑1 (MQ=255)
                                                                                                                    gcTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATaataaaca                                 >  1:718997/1‑150 (MQ=255)
                                                                                                                          aaTTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATaataaacaataaac                           <  2:1364746/150‑1 (MQ=255)
                                                                                                                            ttCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTTTATCACTCTCACAAATTCGCTCAAATAATAAACAATACActg                         >  1:78425/1‑149 (MQ=255)
                                                                                                                            ttCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAATCAATAAActg                         >  1:418445/1‑149 (MQ=255)
                                                                                                                             tCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAActct                        <  1:1114088/150‑1 (MQ=255)
                                                                                                                              cTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACCCCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAActctt                       >  1:1120872/1‑149 (MQ=255)
                                                                                                                               tgtgGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAACTCTGt                      <  2:276075/150‑1 (MQ=255)
                                                                                                                                  gggTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAACTCTGtttt                   >  1:1061264/1‑150 (MQ=255)
                                                                                                                                   ggTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAACTCTGatttt                  >  2:688362/1‑150 (MQ=255)
                                                                                                                                   ggTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAACTCGGttttt                  >  2:1731363/1‑150 (MQ=255)
                                                                                                                                     tCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAACTCTGTTTTTTg                <  2:950249/150‑1 (MQ=255)
                                                                                                                                       ggatgcggCGTGGACGCCTTATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTTCCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAACTCTGTTTTTTGAt              >  2:862117/7‑150 (MQ=255)
                                                                                                                                       ggATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAACTCTGTTTTTTGAt              >  2:619231/1‑150 (MQ=255)
                                                                                                                                       ggATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAAAAATAAACTCTGTTTTTTGAt              >  2:1493922/1‑150 (MQ=255)
                                                                                                                                       gccacccgCGTCCACGCCGCATCCGACAGTCGAGCATCAATCCCGGATGCGCTTCTTATCAGGCCTACCGAACCCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAAAAACCTCTGTTTTTAGAt              <  2:629078/143‑1 (MQ=255)
                                                                                                                                         aTAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAACTCTGTTTTTTGAtct            >  1:554733/1‑150 (MQ=255)
                                                                                                                                          tAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCTCAAATTCGCTCAAATAATAAACAATAAACTCTGTTTTTTGAtctc           <  2:634698/150‑1 (MQ=255)
                                                                                                                                           aaGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAACTCTGTTTTTTGATCTCa          <  1:1551960/150‑1 (MQ=255)
                                                                                                                                            aGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCTCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAACTCTGTTTTTTGATCTCAc         <  1:1446264/150‑1 (MQ=255)
                                                                                                                                            aGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAGCGCTGTTTTTTGATCTCAc         >  2:927496/1‑150 (MQ=255)
                                                                                                                                            aGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAACTCTGTTTTTTGATCTCAc         >  1:132271/1‑150 (MQ=255)
                                                                                                                                             ggCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAACTCTGTTTTTTGATCTCAcc        >  2:1022308/1‑150 (MQ=255)
                                                                                                                                              gCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCATGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAACTCTGTTTTTTGATCTCAccc       <  2:1398415/150‑1 (MQ=255)
                                                                                                                                                gTCCACGCCGCATCCGACAGGCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAACTCTGTTTTTTGATCTCACCCgg     <  1:1667153/150‑1 (MQ=255)
                                                                                                                                                   cACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGACCGCCCTGCATACACCCCTCACTCTATATGACTCTCACAAATTCGCTCAAATAATACACAATAAACTCTGTTTTTTGATCTCACCCGGTaa  >  2:1667040/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GAAGAGCAGACGCACAAATTCAAGAATATCTCCCCGCATTTACAGCGCATTAAAGCGTTCAAAACCCTCGGGTAAATGCCCTCGTCGCATCAGGTAACCTTGCCGGTACCTGATGCGCTCCGAATTCTGTGGGTCGGATAAGGCGTCCACGCCGCATCCGACAGTCGAGCATCAATGCCTGATGCGCTTCTTATCAGGCCTACCGAACGCCCTGCATACACCCCTCACTCTATATCACTCTCACAAATTCGCTCAAATAATAAACAATAAACTCTGTTTTTTGATCTCACCCGGTAA  >  NC_000913/1165611‑1165907

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: