Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,809,895 A→T 11.9% Y591F (TAC→TTC)  katE → catalase HPII, heme d‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,809,8950AT11.9% 735.9 / 77.6 243Y591F (TAC→TTC) katEcatalase HPII, heme d‑containing
Reads supporting (aligned to +/- strand):  ref base A (109/105);  new base T (11/18);  total (120/123)
Fisher's exact test for biased strand distribution p-value = 2.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.69e-01

Only 100 of 249 total aligned reads displayed.

GCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTG  >  CP009273/1809746‑1810043
                                                                                                                                                     |                                                                                                                                                    
gcgcGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTa                                                                                                                                                      <  2:2606591/150‑1 (MQ=255)
gcgcGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTa                                                                                                                                                      <  2:2606393/150‑1 (MQ=255)
  gcgTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACg                                                                                                                                                    <  1:2238115/150‑1 (MQ=255)
      tgttgACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCAt                                                                                                                                                >  2:766206/1‑150 (MQ=255)
         tgACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCTCCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTcc                                                                                                                                             >  2:125469/1‑150 (MQ=255)
          gACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCt                                                                                                                                            <  1:2036484/150‑1 (MQ=255)
           aCCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTg                                                                                                                                           <  2:852488/150‑1 (MQ=255)
           aCCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTg                                                                                                                                           >  2:3706376/1‑150 (MQ=255)
           aCCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTg                                                                                                                                           >  2:1382355/1‑150 (MQ=255)
             cAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGAc                                                                                                                                         >  1:602445/1‑150 (MQ=255)
              aGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTTCGCCATTCCTGACg                                                                                                                                        <  2:3291872/150‑1 (MQ=255)
              aGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACg                                                                                                                                        <  2:187994/150‑1 (MQ=255)
                cTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGt                                                                                                                                      >  1:3740169/1‑150 (MQ=255)
                   gcgcATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGAt                                                                                                                                   <  2:2982509/150‑1 (MQ=255)
                   gcgcATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGAt                                                                                                                                   <  2:4380449/150‑1 (MQ=255)
                    cgcATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTTCGCCATTCCTGACGGTGAtg                                                                                                                                  <  2:711879/150‑1 (MQ=255)
                    cgcATATTGATCTCACTCTGGCCCAGGCGGAGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCATAAGTTTGTACGCCATTCCTGACGGTGAtg                                                                                                                                  >  1:2256771/1‑150 (MQ=255)
                        tatTGATCTAACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATAACCCCACATCCGGACGTCAACGGTCTGAAAAAGTATCCATCCTTAAGTTTGTACTCCATTCCTGACGGTGATGTGaa                                                                                                                              >  1:2144762/1‑150 (MQ=255)
                             aTCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTc                                                                                                                         >  2:1097205/1‑150 (MQ=255)
                              tctcACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTTCGCCATTCCTGACGGTGATGTGAAAGGTcg                                                                                                                        <  2:389098/150‑1 (MQ=255)
                                   ctctGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTgg                                                                                                                   <  2:1323713/150‑1 (MQ=255)
                                   ctctGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTgg                                                                                                                   >  1:4248010/1‑150 (MQ=255)
                                    tctGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGt                                                                                                                  <  2:198170/150‑1 (MQ=255)
                                     ctGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTTCGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTa                                                                                                                 <  1:1597922/150‑1 (MQ=255)
                                          ccAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAt                                                                                                            >  1:1158387/1‑150 (MQ=255)
                                            aGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGGTGTGAAAGGTCGCGTGGTAGCGAttt                                                                                                          >  2:3908220/1‑150 (MQ=255)
                                            aGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGAttt                                                                                                          >  1:2297795/1‑150 (MQ=255)
                                                   ggCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGTACGACAACGGTCTGAAAAAGGATCCAACCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTaa                                                                                                   <  1:4529793/150‑1 (MQ=255)
                                                   ggCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTaa                                                                                                   <  1:2197193/150‑1 (MQ=255)
                                                    gCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTTCGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAat                                                                                                  <  1:3394433/150‑1 (MQ=255)
                                                    gCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAat                                                                                                  <  2:2788006/150‑1 (MQ=255)
                                                     cGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAatg                                                                                                 >  2:1653116/1‑150 (MQ=255)
                                                       aaaaaTCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTTCGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAatgat                                                                                               <  2:4260244/150‑1 (MQ=255)
                                                       aaaaaTCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAatgat                                                                                               >  1:2032927/1‑150 (MQ=255)
                                                              tcGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTga                                                                                        >  2:3473417/1‑150 (MQ=255)
                                                                ggTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTgaga                                                                                      >  1:3845888/1‑150 (MQ=255)
                                                                ggTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTgaga                                                                                      >  1:2419929/1‑150 (MQ=255)
                                                                 gTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTTCGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGAt                                                                                     <  1:1965170/150‑1 (MQ=255)
                                                                 gTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGAt                                                                                     >  1:1389683/1‑150 (MQ=255)
                                                                    tCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCgg                                                                                  <  1:4224272/150‑1 (MQ=255)
                                                                      gAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAGAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCa                                                                                <  1:1533366/150‑1 (MQ=255)
                                                                      gAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCa                                                                                >  2:2054803/1‑150 (MQ=255)
                                                                            actgacGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCtt                                                                          <  2:3412365/150‑1 (MQ=255)
                                                                              tgacGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCt                                                                        <  1:2347912/150‑1 (MQ=255)
                                                                              tgacGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCt                                                                        <  2:588644/150‑1 (MQ=255)
                                                                                acgacCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTgg                                                                      >  1:1426539/1‑150 (MQ=255)
                                                                                 cgacCAGCTGAATATCACACCACCTACGGACGTCAACGGTCTGAAAAGGGATCCATCCTTAAGTTTGTACGCCATTACTGACGGTGATGTGAAATGTCGCGTGTTAGCGATTTTACTAAATGATGAAGTGAGAACGGCAGACCTTCTGGc                                                                     >  2:1722876/1‑150 (MQ=255)
                                                                                  gacCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGcc                                                                    <  2:4178155/150‑1 (MQ=255)
                                                                                    ccAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCAt                                                                  <  2:2845714/150‑1 (MQ=255)
                                                                                      aGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTTCGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATtc                                                                >  1:2511509/1‑150 (MQ=255)
                                                                                      aGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATtc                                                                >  1:2729913/1‑150 (MQ=255)
                                                                                       gCTGAATATAACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATtct                                                               <  2:4210473/150‑1 (MQ=255)
                                                                                            atatCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAgg                                                          >  1:3697264/1‑150 (MQ=255)
                                                                                              atCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGgcg                                                        >  2:1437132/1‑150 (MQ=255)
                                                                                                cACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTTCGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCt                                                      >  2:1180886/1‑150 (MQ=255)
                                                                                                 aCCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTg                                                     <  1:3870308/150‑1 (MQ=255)
                                                                                                    ccACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAg                                                  <  2:3742982/150‑1 (MQ=255)
                                                                                                     cACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAgg                                                 <  1:4334147/150‑1 (MQ=255)
                                                                                                         tCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCaa                                             >  2:2591323/1‑150 (MQ=255)
                                                                                                         tCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCaa                                             <  1:252860/150‑1 (MQ=255)
                                                                                                           cGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAg                                           <  1:2925583/150‑1 (MQ=255)
                                                                                                           cGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAg                                           <  1:821695/150‑1 (MQ=255)
                                                                                                               cGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTTCGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGt                                       <  2:2118813/150‑1 (MQ=255)
                                                                                                                gTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGtt                                      >  2:4478756/1‑150 (MQ=255)
                                                                                                                gTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGtt                                      >  1:4054984/1‑150 (MQ=255)
                                                                                                                 tCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTc                                     >  1:1632853/1‑150 (MQ=255)
                                                                                                                 tCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTc                                     >  2:4547496/1‑150 (MQ=255)
                                                                                                                 tCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCATGCCGGTGATGTGAAAGGTCGCGTGGTCGCGCTTTTACTTAATGATGACGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGtgc                                     >  2:3974485/1‑148 (MQ=255)
                                                                                                                  cAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCa                                    >  1:3227007/1‑150 (MQ=255)
                                                                                                                   aaCGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCAt                                   <  2:1449505/150‑1 (MQ=255)
                                                                                                                     cGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGc                                 <  2:4341117/150‑1 (MQ=255)
                                                                                                                        tCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCaa                              >  1:4022941/1‑150 (MQ=255)
                                                                                                                          tGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAAc                            <  1:620194/150‑1 (MQ=255)
                                                                                                                          tGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAAc                            >  1:3551085/1‑150 (MQ=255)
                                                                                                                           gAAAAAGGATCCATCCTTAAGTTTGTTCGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACt                           >  2:3425629/1‑150 (MQ=255)
                                                                                                                           gAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACt                           <  1:182360/150‑1 (MQ=255)
                                                                                                                           gAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACt                           >  1:2871285/1‑150 (MQ=255)
                                                                                                                            aaaaaGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTg                          <  1:2846962/150‑1 (MQ=255)
                                                                                                                             aaaaGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGc                         <  1:1747653/150‑1 (MQ=255)
                                                                                                                                aGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGctct                      >  1:1852785/1‑150 (MQ=255)
                                                                                                                                aGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGctct                      >  1:1503832/1‑150 (MQ=255)
                                                                                                                                aGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACATAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGctct                      <  2:545223/150‑1 (MQ=255)
                                                                                                                                   atccatccTTAAGTTTGTTCGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACt                   >  2:2182643/1‑150 (MQ=255)
                                                                                                                                   atccatccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACt                   <  2:3560398/150‑1 (MQ=255)
                                                                                                                                    tccatccTTAAGTTTGTTCGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTc                  <  2:909867/150‑1 (MQ=255)
                                                                                                                                    tccatccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTc                  <  1:772792/150‑1 (MQ=255)
                                                                                                                                     ccatccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTcc                 >  2:2831772/1‑150 (MQ=255)
                                                                                                                                     ccatccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTcc                 <  1:1077279/150‑1 (MQ=255)
                                                                                                                                      catccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTccc                >  1:2816614/1‑150 (MQ=255)
                                                                                                                                      catccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTccc                <  1:2732334/150‑1 (MQ=255)
                                                                                                                                        tccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGa              >  1:1782858/1‑150 (MQ=255)
                                                                                                                                         ccTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATTATGAAGTGAGATCTGCAGACCTTCTGGCCAATCTCAAGGCGCTGAAGGCCAAAGGCTTTCATGCCAAACTGCTCTACTCCCGac             >  2:1534478/1‑149 (MQ=255)
                                                                                                                                          cTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAAt            <  1:1176676/150‑1 (MQ=255)
                                                                                                                                           ttAAGTTTGTTCGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATg           >  1:3417092/1‑150 (MQ=255)
                                                                                                                                           ttAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATg           >  1:1685093/1‑150 (MQ=255)
                                                                                                                                              aGTTTGTTCGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGt        >  2:1791295/1‑150 (MQ=255)
                                                                                                                                               gTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACGCCCGAATGGGTg       <  2:2256771/150‑1 (MQ=255)
                                                                                                                                                 ttGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGaa     <  2:2132464/150‑1 (MQ=255)
                                                                                                                                                 ttGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGaa     <  1:68856/150‑1 (MQ=255)
                                                                                                                                                    tACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTg  <  1:2528162/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGCGTTGTTGACCAGCTGGCGCATATTGATCTCACTCTGGCCCAGGCGGTGGCGAAAAATCTCGGTATCGAACTGACTGACGACCAGCTGAATATCACCCCACCTCCGGACGTCAACGGTCTGAAAAAGGATCCATCCTTAAGTTTGTACGCCATTCCTGACGGTGATGTGAAAGGTCGCGTGGTAGCGATTTTACTTAATGATGAAGTGAGATCGGCAGACCTTCTGGCCATTCTCAAGGCGCTGAAGGCCAAAGGCGTTCATGCCAAACTGCTCTACTCCCGAATGGGTGAAGTG  >  CP009273/1809746‑1810043

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: