Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,294,769 C→A 5.6% P488P (CCG→CCT napA ← nitrate reductase, periplasmic, large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,294,7690CA 5.6% 1240.6 / 33.3 319P488P (CCG→CCTnapAnitrate reductase, periplasmic, large subunit
Reads supporting (aligned to +/- strand):  ref base C (138/163);  new base A (9/9);  total (147/172)
Fisher's exact test for biased strand distribution p-value = 8.10e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 338 total aligned reads displayed.

GCGTTACCGTAAGCGCCCTCTTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCACTTCTTCTCGCAGATATCAC  >  CP009273/2294620‑2294917
                                                                                                                                                     |                                                                                                                                                    
gcgttACCGTAAGCGCCCTCTTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTc                                                                                                                                                      <  2:3081403/150‑1 (MQ=255)
 cgttACCGTAAGCGCCCTCTTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTcg                                                                                                                                                     >  2:1125139/1‑150 (MQ=255)
  gttACCGTAAGCGCCCTCTTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCCTTGATGTTcgg                                                                                                                                                    >  2:2565792/1‑150 (MQ=255)
     aCCGTAAGCGCCCTCTTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTcggccc                                                                                                                                                 >  2:3941520/1‑150 (MQ=255)
         tAAGCGCCCTCTTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTcggcccggcc                                                                                                                                             >  2:3362171/1‑150 (MQ=255)
         tAAGCGCCCTCTTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTcggcccggcc                                                                                                                                             >  1:2439950/1‑150 (MQ=255)
               cccTCTTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAtg                                                                                                                                       >  1:295186/1‑150 (MQ=255)
                ccTCGTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAtgt                                                                                                                                      <  2:59400/150‑1 (MQ=255)
                 ctctTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAtgtt                                                                                                                                     <  1:3929884/150‑1 (MQ=255)
                  tctTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCAtgttg                                                                                                                                    >  1:3255223/1‑150 (MQ=255)
                    tttCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATgttgtt                                                                                                                                  >  2:386368/1‑150 (MQ=255)
                       ctctACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGt                                                                                                                               <  1:2824664/150‑1 (MQ=255)
                         ctACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTac                                                                                                                             >  2:1507385/1‑150 (MQ=255)
                          tACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTaca                                                                                                                            <  1:4537115/150‑1 (MQ=255)
                            cccACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTAcaca                                                                                                                          >  1:3071376/1‑150 (MQ=255)
                            cccACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTAcaca                                                                                                                          >  1:2189593/1‑150 (MQ=255)
                             ccACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACAt                                                                                                                         >  1:550139/1‑150 (MQ=255)
                             ccACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACAt                                                                                                                         <  1:4320795/150‑1 (MQ=255)
                                 aTTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTc                                                                                                                     <  1:4147710/150‑1 (MQ=255)
                                   tGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCa                                                                                                                   >  2:92055/1‑150 (MQ=255)
                                     cGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGt                                                                                                                 <  1:1891631/150‑1 (MQ=255)
                                         cGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAAc                                                                                                             >  1:905389/1‑150 (MQ=255)
                                         cGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAAc                                                                                                             <  1:1041233/150‑1 (MQ=255)
                                         cGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTAGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAAc                                                                                                             >  2:2220594/1‑150 (MQ=255)
                                          ggCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTAGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACa                                                                                                            >  1:2085440/1‑150 (MQ=255)
                                             aGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTc                                                                                                         >  2:1889206/1‑150 (MQ=255)
                                             aGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTc                                                                                                         >  1:3631373/1‑150 (MQ=255)
                                                 tCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGtt                                                                                                     >  2:1773344/1‑150 (MQ=255)
                                                 tCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGtt                                                                                                     >  2:849131/1‑150 (MQ=255)
                                                  cAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGttt                                                                                                    <  1:4456142/150‑1 (MQ=255)
                                                   aaGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTg                                                                                                   >  1:1829685/1‑150 (MQ=255)
                                                     gTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACAGTCAGTTTGcc                                                                                                 >  1:145101/1‑150 (MQ=255)
                                                      tCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCAGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCg                                                                                                >  1:2038268/1‑150 (MQ=255)
                                                       cggcggCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGt                                                                                               >  2:3829386/1‑150 (MQ=255)
                                                        ggcggcCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTc                                                                                              >  1:3181867/1‑150 (MQ=255)
                                                         gcggcCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCt                                                                                             <  2:3652743/150‑1 (MQ=255)
                                                             ccAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCa                                                                                         <  2:3523837/150‑1 (MQ=255)
                                                              cAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAg                                                                                        <  1:3357809/150‑1 (MQ=255)
                                                               aGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAgc                                                                                       <  2:2636705/150‑1 (MQ=255)
                                                               aGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAgc                                                                                       <  2:1624132/150‑1 (MQ=255)
                                                                  gcACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCa                                                                                    <  2:248793/150‑1 (MQ=255)
                                                                      tgactgTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGt                                                                                <  2:1757693/150‑1 (MQ=255)
                                                                       gactgTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTc                                                                               >  1:42515/1‑150 (MQ=255)
                                                                         ctgTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCtt                                                                             <  1:2569082/150‑1 (MQ=255)
                                                                         ctgTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCtt                                                                             <  1:2725182/150‑1 (MQ=255)
                                                                         ctgTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCtt                                                                             <  2:1089460/150‑1 (MQ=255)
                                                                           gTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTtgt                                                                           <  1:1436473/150‑1 (MQ=255)
                                                                            tCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTtgtg                                                                          <  2:2072049/150‑1 (MQ=255)
                                                                              ggATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGcc                                                                        <  1:4443411/150‑1 (MQ=255)
                                                                              ggATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGcc                                                                        >  2:1351624/1‑150 (MQ=255)
                                                                                aTACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCAc                                                                      >  2:2796736/1‑150 (MQ=255)
                                                                                  aCGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACcg                                                                    >  1:687122/1‑150 (MQ=255)
                                                                                       tCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGc                                                               >  2:4556112/1‑150 (MQ=255)
                                                                                        cGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCa                                                              >  1:1316605/1‑150 (MQ=255)
                                                                                        cGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCa                                                              >  2:373549/1‑150 (MQ=255)
                                                                                        cGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCa                                                              <  1:4603657/150‑1 (MQ=255)
                                                                                         ggAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATCCGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAg                                                             <  1:3552818/150‑1 (MQ=255)
                                                                                         ggAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAg                                                             <  2:4383878/150‑1 (MQ=255)
                                                                                         ggAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAg                                                             <  2:1984104/150‑1 (MQ=255)
                                                                                         ggAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTAGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAg                                                             <  1:3844993/150‑1 (MQ=255)
                                                                                          gagaCGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGa                                                            <  2:4763052/150‑1 (MQ=255)
                                                                                            gaCGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGAcc                                                          <  2:4692282/150‑1 (MQ=255)
                                                                                              cGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGa                                                        <  2:3277129/150‑1 (MQ=255)
                                                                                                atgatgAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTAGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGAtt                                                      <  2:1233450/150‑1 (MQ=255)
                                                                                                  gatgaAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGAtttt                                                    <  2:3199647/150‑1 (MQ=255)
                                                                                                  gatgaAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGAtttt                                                    <  2:4063511/150‑1 (MQ=255)
                                                                                                  gatgaAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGAtttt                                                    <  2:2802642/150‑1 (MQ=255)
                                                                                                    tgaAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCg                                                  <  2:1792700/150‑1 (MQ=255)
                                                                                                       aGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCg                                               >  1:504454/1‑150 (MQ=255)
                                                                                                       aGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCg                                               >  1:1341616/1‑150 (MQ=255)
                                                                                                       aGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCg                                               >  2:171390/1‑150 (MQ=255)
                                                                                                           gcgcGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAAt                                           <  2:4292716/150‑1 (MQ=255)
                                                                                                           gcgcGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAAt                                           >  2:4324932/1‑150 (MQ=255)
                                                                                                             gcgGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATgg                                         <  2:2091411/150‑1 (MQ=255)
                                                                                                                gATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGc                                      >  1:578663/1‑150 (MQ=255)
                                                                                                                 aTCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGcc                                     >  2:2996045/1‑150 (MQ=255)
                                                                                                                  tCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCg                                    <  1:235571/150‑1 (MQ=255)
                                                                                                                     cgcCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTc                                 >  1:829773/1‑150 (MQ=255)
                                                                                                                         aGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGa                             >  1:1388389/1‑150 (MQ=255)
                                                                                                                           cccGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGata                           <  2:3223085/150‑1 (MQ=255)
                                                                                                                             cGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATAtt                         >  1:2882027/1‑150 (MQ=255)
                                                                                                                               gCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTcc                       <  1:1016716/150‑1 (MQ=255)
                                                                                                                               gCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTcc                       <  1:1764952/150‑1 (MQ=255)
                                                                                                                                cATACGCTCTTCATTGATGTTAGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCa                      >  1:3654601/1‑150 (MQ=255)
                                                                                                                                  tACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCAct                    >  2:1497692/1‑150 (MQ=255)
                                                                                                                                   aCGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCActt                   >  1:1529702/1‑150 (MQ=255)
                                                                                                                                     gCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCActtct                 <  2:2345487/150‑1 (MQ=255)
                                                                                                                                     gCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCActtct                 <  2:3848167/150‑1 (MQ=255)
                                                                                                                                      ctctTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCActtctt                <  1:114924/150‑1 (MQ=255)
                                                                                                                                      ctctTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCActtctt                <  2:777921/150‑1 (MQ=255)
                                                                                                                                        ctTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCACTtcttct              >  1:756234/1‑150 (MQ=255)
                                                                                                                                        ctTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCACTtcttct              >  2:4614805/1‑150 (MQ=255)
                                                                                                                                        ctTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCACTtcttct              >  1:2796330/1‑150 (MQ=255)
                                                                                                                                         ttCATTGATGTTAGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCACTTCTtctc             <  2:52066/150‑1 (MQ=255)
                                                                                                                                          tCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCACTTCTTCTCg            >  2:1013855/1‑150 (MQ=255)
                                                                                                                                             ttGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCACTTCTTCTCGCAg         <  2:1292366/150‑1 (MQ=255)
                                                                                                                                               gATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCACTTCTTCTCGCAGat       <  2:3631589/150‑1 (MQ=255)
                                                                                                                                                aTGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCACTTCTTCTCGCAGata      <  2:4151318/150‑1 (MQ=255)
                                                                                                                                                 tGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCACTTCTTCTCGCAGatat     >  2:2900886/1‑150 (MQ=255)
                                                                                                                                                   ttCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCACTTCTTCTCGCAGATATCa   <  1:4352755/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCGTTACCGTAAGCGCCCTCTTTCTCTACCCACATTGCGGTCGGCAGGATCAAGTCGGCGGCCAGCGCACTGACTGTCGGATACGGATCGGAGACGATGATGAAGTTGCGCGGATCGCGCCAGCCCGGCATACGCTCTTCATTGATGTTCGGCCCGGCCTGCATGTTGTTGGTACACATGGTCCAGTAAACATTCAGTTTGCCGTCTTTCAGCGCACGGTCTTGTGCCACCGCGTGCAGACCGATTTTCGCCGGAATGGTGCCGCTCGGGATATTCCACTTCTTCTCGCAGATATCAC  >  CP009273/2294620‑2294917

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: