Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,291,194 C→A 20.0% C182F (TGC→TTC)  napC ← quinol dehydrogenase, electron source for NapAB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,291,1940CA20.0% 914.6 / 151.5 314C182F (TGC→TTC) napCquinol dehydrogenase, electron source for NapAB
Reads supporting (aligned to +/- strand):  ref base C (152/99);  new base A (35/28);  total (187/127)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 328 total aligned reads displayed.

TCAGCGTAAATGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGCGCCGCCACTCATTCTG  >  CP009273/2291045‑2291342
                                                                                                                                                     |                                                                                                                                                    
tCAGCGTAAATGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGc                                                                                                                                                      <  1:3401176/150‑1 (MQ=255)
  aGCGTAAATGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGaaa                                                                                                                                                    >  1:39124/1‑150 (MQ=255)
      tAAATGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGcaatcaa                                                                                                                                                >  2:2982432/1‑150 (MQ=255)
        aaTGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATa                                                                                                                                              >  1:4303062/1‑150 (MQ=255)
         aTGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATAc                                                                                                                                             <  1:4003871/150‑1 (MQ=255)
         aTGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATAc                                                                                                                                             >  1:3676888/1‑150 (MQ=255)
          tGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACa                                                                                                                                            >  1:4469742/1‑150 (MQ=255)
               aGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGttt                                                                                                                                       >  2:2518388/1‑150 (MQ=255)
               aGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGttt                                                                                                                                       >  2:2851333/1‑150 (MQ=255)
                gCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAAACTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTATATGGCAATCAATACAGGTTTg                                                                                                                                      >  2:4237723/1‑150 (MQ=255)
                  cACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGcc                                                                                                                                    >  2:1994951/1‑150 (MQ=255)
                    cTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCa                                                                                                                                  >  2:4009702/1‑150 (MQ=255)
                    cTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCa                                                                                                                                  >  2:3595908/1‑150 (MQ=255)
                     tAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCAt                                                                                                                                 >  2:2144488/1‑150 (MQ=255)
                     tAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCAt                                                                                                                                 <  2:4137744/150‑1 (MQ=255)
                      aaaTAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATc                                                                                                                                >  2:3518265/1‑150 (MQ=255)
                       aaTAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCt                                                                                                                               <  1:415183/150‑1 (MQ=255)
                        aTAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCtt                                                                                                                              >  2:750189/1‑150 (MQ=255)
                         tAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCttt                                                                                                                             >  2:2300768/1‑150 (MQ=255)
                          aaGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTc                                                                                                                            >  2:3081249/1‑150 (MQ=255)
                            ggTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCAc                                                                                                                          >  1:2729942/1‑150 (MQ=255)
                            ggTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCAc                                                                                                                          >  2:1019517/1‑150 (MQ=255)
                             gttcgttcATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCAcc                                                                                                                         >  1:4596399/1‑150 (MQ=255)
                                 gttcATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCt                                                                                                                     >  1:770969/1‑150 (MQ=255)
                                 gttcATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCt                                                                                                                     <  2:1486670/150‑1 (MQ=255)
                                 gttcATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCt                                                                                                                     <  1:1261482/150‑1 (MQ=255)
                                  ttcATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATCTCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTg                                                                                                                    <  2:2778118/150‑1 (MQ=255)
                                  ttcATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTg                                                                                                                    >  1:2248309/1‑150 (MQ=255)
                                        ccGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATg                                                                                                              >  2:4070145/1‑150 (MQ=255)
                                              acaAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCtt                                                                                                        >  1:3294989/1‑150 (MQ=255)
                                               caAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTc                                                                                                       <  1:445572/150‑1 (MQ=255)
                                               caAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTc                                                                                                       >  1:4226100/1‑150 (MQ=255)
                                               caAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTc                                                                                                       >  1:1078697/1‑150 (MQ=255)
                                                aaaGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTcg                                                                                                      <  1:1485939/150‑1 (MQ=255)
                                                 aaGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTcgc                                                                                                     >  2:1222066/1‑150 (MQ=255)
                                                 aaGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTcgc                                                                                                     >  2:4190371/1‑150 (MQ=255)
                                                       ctctctGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAAc                                                                                               <  2:25095/150‑1 (MQ=255)
                                                        tctctGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAAcc                                                                                              >  1:3806114/1‑150 (MQ=255)
                                                         ctctGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCg                                                                                             >  2:3931215/1‑150 (MQ=255)
                                                          tctGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGa                                                                                            <  2:773131/150‑1 (MQ=255)
                                                          tctGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGa                                                                                            <  1:216397/150‑1 (MQ=255)
                                                           ctGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGAt                                                                                           >  2:3762149/1‑150 (MQ=255)
                                                                   aaGCATACCCACGCAAGAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGc                                                                                   <  1:2414640/150‑1 (MQ=255)
                                                                    aGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCg                                                                                  >  2:4522637/1‑150 (MQ=255)
                                                                       aTACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGtt                                                                               >  2:2132110/1‑150 (MQ=255)
                                                                          cccACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTa                                                                            <  2:3689160/150‑1 (MQ=255)
                                                                          cccACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTa                                                                            >  2:3952155/1‑150 (MQ=255)
                                                                                cAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCata                                                                      >  1:1425738/1‑150 (MQ=255)
                                                                                 aaTAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCatat                                                                     <  1:2326893/150‑1 (MQ=255)
                                                                                  aTAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCAtata                                                                    <  2:3356698/150‑1 (MQ=255)
                                                                                  aTAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCAtata                                                                    >  1:1843692/1‑150 (MQ=255)
                                                                                   tAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATAc                                                                   >  1:4468506/1‑150 (MQ=255)
                                                                                   tAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATAc                                                                   <  2:1245951/150‑1 (MQ=255)
                                                                                    aaCCCTGTTAAAAACCTGGCTCGACTTCTCGCATATCGGGCAGCTTGTGAGCTATCCCTTTATGGAAATCAATAAATGTTTGCACATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACt                                                                  <  1:818971/150‑1 (MQ=255)
                                                                                    aaCCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACt                                                                  <  2:3637345/150‑1 (MQ=255)
                                                                                    aaCCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACt                                                                  <  2:2597374/150‑1 (MQ=255)
                                                                                     aCCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTc                                                                 >  2:635564/1‑150 (MQ=255)
                                                                                       ccTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGa                                                               >  1:3542515/1‑150 (MQ=255)
                                                                                       ccTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGa                                                               <  2:4620992/150‑1 (MQ=255)
                                                                                       ccTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGa                                                               >  1:563878/1‑150 (MQ=255)
                                                                                       ccTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGa                                                               >  2:1527274/1‑150 (MQ=255)
                                                                                         tGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAg                                                             >  2:4299746/1‑150 (MQ=255)
                                                                                         tGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAg                                                             >  1:4415324/1‑150 (MQ=255)
                                                                                           ttAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGtt                                                           <  2:2045903/150‑1 (MQ=255)
                                                                                           ttAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGtt                                                           <  1:834648/150‑1 (MQ=255)
                                                                                           ttAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGtt                                                           <  1:2518224/150‑1 (MQ=255)
                                                                                             aaaaaCCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTtgt                                                         >  1:1530489/1‑150 (MQ=255)
                                                                                                aaCCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGAc                                                      >  2:2711606/1‑150 (MQ=255)
                                                                                                aaCCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGAc                                                      >  2:3098215/1‑150 (MQ=255)
                                                                                                 aCCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACa                                                     >  2:3244696/1‑150 (MQ=255)
                                                                                                    tGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGtt                                                  >  1:3620415/1‑150 (MQ=255)
                                                                                                     ggCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTa                                                 <  2:4577055/150‑1 (MQ=255)
                                                                                                     ggCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTa                                                 >  2:1333184/1‑150 (MQ=255)
                                                                                                      gCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTAc                                                <  2:2440688/150‑1 (MQ=255)
                                                                                                        tCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACgg                                              >  2:1659867/1‑150 (MQ=255)
                                                                                                        tCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACgg                                              >  1:2820108/1‑150 (MQ=255)
                                                                                                           aCTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCAc                                           >  2:1149745/1‑150 (MQ=255)
                                                                                                            cTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACt                                          >  1:3522852/1‑150 (MQ=255)
                                                                                                              tCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTcc                                        <  2:4350120/150‑1 (MQ=255)
                                                                                                              tCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTcc                                        <  2:1414076/150‑1 (MQ=255)
                                                                                                                aCGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTg                                      >  1:1991606/1‑150 (MQ=255)
                                                                                                                  gCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCg                                    <  1:2577542/150‑1 (MQ=255)
                                                                                                                   cATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGa                                   >  1:4636305/1‑150 (MQ=255)
                                                                                                                   cATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGa                                   <  1:1810393/150‑1 (MQ=255)
                                                                                                                    atatCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAg                                  >  2:4562075/1‑150 (MQ=255)
                                                                                                                             aGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTcc                         >  2:1592451/1‑150 (MQ=255)
                                                                                                                              gCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCt                        <  1:1331484/150‑1 (MQ=255)
                                                                                                                                 tgtgCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCa                     <  1:2002313/150‑1 (MQ=255)
                                                                                                                                  gtgCGCTATCCCTTTATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCAt                    <  1:1607209/150‑1 (MQ=255)
                                                                                                                                    gcgcTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATgc                  <  2:182076/150‑1 (MQ=255)
                                                                                                                                        tATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGcgccg              <  1:1040093/150‑1 (MQ=255)
                                                                                                                                              tttATGGAAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGCGCCGCCACTc        <  1:3528082/150‑1 (MQ=255)
                                                                                                                                               ttATGGCCATCAATACAGGTTTTCCCATCTTTCACCGCCTGGTCATTCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGCGCCGCCACTCa       <  2:3830210/150‑1 (MQ=255)
                                                                                                                                                tATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGCGCCGCCACTCAt      >  1:4517458/1‑150 (MQ=255)
                                                                                                                                                 aTGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGCGCCGCCACTCAtt     >  2:3505067/1‑150 (MQ=255)
                                                                                                                                                  tGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCTTGCGCCGCCACTCATTc    >  2:3069134/1‑150 (MQ=255)
                                                                                                                                                  tGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGCGCCGCCACTCATTc    >  1:313341/1‑150 (MQ=255)
                                                                                                                                                   ggCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGCGCCGCCACTCATTCt   <  1:211398/150‑1 (MQ=255)
                                                                                                                                                   ggCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGCGCCGCCACTCATTCt   <  1:1680192/150‑1 (MQ=255)
                                                                                                                                                    gCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGCGCCGCCACTCATTCtg  <  1:981564/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCAGCGTAAATGACAAGCCACTAAATAAGGTTCGTTCATCCCGCTCACAAAGTAACTCTCTGGCTTCAAGCATACCCACGCAATAACCCTGTTAAAAACCTGGCTCGACTTCACGCATATCGGGCAGCTTGTGCGCTATCCCTTTATGGCAATCAATACAGGTTTGCCCATCTTTCACCGCCTGGTCATGCATCTTCGCGGCAACCGATTTCTGGGCGGTTGTATCCATATACTCGAAGTTGTGACAGTTACGGCACTCCTGCGAGTTATTGTCCTTCATGCGCCGCCACTCATTCTG  >  CP009273/2291045‑2291342

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: