Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,443,196 A→T 10.2% Q474L (CAG→CTG)  feaB → phenylacetaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,443,1960AT10.2% 808.1 / 67.8 254Q474L (CAG→CTG) feaBphenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (115/113);  new base T (15/10);  total (130/124)
Fisher's exact test for biased strand distribution p-value = 4.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.32e-01

Only 100 of 257 total aligned reads displayed.

AAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTC  >  CP009273/1443047‑1443345
                                                                                                                                                     |                                                                                                                                                     
 aaCGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTCATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAg                                                                                                                                                      >  2:4133537/1‑150 (MQ=255)
   cGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTc                                                                                                                                                    <  1:785340/150‑1 (MQ=255)
   cGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTc                                                                                                                                                    >  2:1628813/1‑150 (MQ=255)
    aacacGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCa                                                                                                                                                   <  1:4355695/149‑1 (MQ=255)
      cacGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTCACTTACCGTTTGGTGGGATGAAGCAGTCAgg                                                                                                                                                 >  2:683573/1‑150 (MQ=255)
      cacGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAgg                                                                                                                                                 >  2:4157487/1‑150 (MQ=255)
      cacGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTGTGGTGGGATGAAGCAGTCAgg                                                                                                                                                 >  2:640613/1‑150 (MQ=255)
         ggAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAAc                                                                                                                                              >  2:336560/1‑150 (MQ=255)
          gAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACg                                                                                                                                             >  2:948447/1‑150 (MQ=255)
          gAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACg                                                                                                                                             >  1:4146172/1‑150 (MQ=255)
           aaTATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACgg                                                                                                                                            <  1:1361433/150‑1 (MQ=255)
            atatGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACggg                                                                                                                                           >  1:833820/1‑150 (MQ=255)
            atatGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACggg                                                                                                                                           <  1:1125229/150‑1 (MQ=255)
                ggCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGt                                                                                                                                       >  2:3993843/1‑150 (MQ=255)
                ggCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGt                                                                                                                                       >  2:1806520/1‑150 (MQ=255)
                   ttAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGAt                                                                                                                                    >  2:4043678/1‑150 (MQ=255)
                     aaCTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGAttt                                                                                                                                  >  2:3592089/1‑150 (MQ=255)
                         gCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGGTGAAGCAGTCAGGAACGGGCCGTGATTTTGGc                                                                                                                              <  2:3270256/150‑1 (MQ=255)
                         gCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGc                                                                                                                              <  1:1453643/150‑1 (MQ=255)
                              tgtCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGa                                                                                                                         <  1:1815945/150‑1 (MQ=255)
                                   ggACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActggc                                                                                                                    <  2:798538/150‑1 (MQ=255)
                                    gACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActggct                                                                                                                   <  2:1893066/150‑1 (MQ=255)
                                        cAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGAc                                                                                                               >  1:3691963/1‑150 (MQ=255)
                                         aaaaTATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACg                                                                                                              >  2:111505/1‑150 (MQ=255)
                                           actctctcCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGt                                                                                                            >  2:1543100/3‑150 (MQ=255)
                                              ctctcCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTgg                                                                                                         <  2:225506/150‑1 (MQ=255)
                                              ctctcCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTgg                                                                                                         <  2:608495/150‑1 (MQ=255)
                                              ctctcCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTgg                                                                                                         >  1:3682753/1‑150 (MQ=255)
                                              ctctcCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTgg                                                                                                         >  1:3512457/1‑150 (MQ=255)
                                               tctcCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGAGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgt                                                                                                        >  1:1249309/1‑150 (MQ=255)
                                                 tcCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGgtgt                                                                                                      >  1:4015266/1‑150 (MQ=255)
                                                    cAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGaa                                                                                                   >  1:3452439/1‑150 (MQ=255)
                                                      ggCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAAc                                                                                                 <  1:1021207/150‑1 (MQ=255)
                                                      ggCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAAc                                                                                                 <  1:3635974/150‑1 (MQ=255)
                                                      ggCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAAc                                                                                                 <  1:2898435/150‑1 (MQ=255)
                                                      ggCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACGTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAAc                                                                                                 <  1:3635926/150‑1 (MQ=255)
                                                       gCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACt                                                                                                <  1:3827724/150‑1 (MQ=255)
                                                         tctGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTaa                                                                                              <  1:3054137/150‑1 (MQ=255)
                                                             gAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCg                                                                                          <  2:2448285/150‑1 (MQ=255)
                                                              aaTATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCgg                                                                                         >  1:497148/1‑150 (MQ=255)
                                                               atatAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgt                                                                                        <  1:575735/150‑1 (MQ=255)
                                                                tataGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtg                                                                                       >  1:1544499/1‑150 (MQ=255)
                                                                  taGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtg                                                                                     >  1:94564/1‑150 (MQ=255)
                                                                     cGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtgtgt                                                                                  >  2:183642/1‑150 (MQ=255)
                                                                     cGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGgtgtgtgt                                                                                  >  1:375122/1‑150 (MQ=255)
                                                                      gATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTa                                                                                 >  1:2818323/1‑150 (MQ=255)
                                                                        tCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACg                                                                               >  1:2945468/1‑150 (MQ=255)
                                                                         cGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACgg                                                                              >  2:1872111/1‑150 (MQ=255)
                                                                          gCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGt                                                                             <  1:3390141/150‑1 (MQ=255)
                                                                          gCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGt                                                                             >  1:96214/1‑150 (MQ=255)
                                                                             tACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTAtt                                                                          >  1:4304310/1‑150 (MQ=255)
                                                                             tACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTAtt                                                                          <  2:369877/150‑1 (MQ=255)
                                                                             gaCAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTAtt                                                                          <  1:681963/149‑1 (MQ=255)
                                                                              aCAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTa                                                                         <  1:1435390/150‑1 (MQ=255)
                                                                              aCAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTa                                                                         <  2:4236205/150‑1 (MQ=255)
                                                                                      ggACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGtt                                                                 <  2:1201702/150‑1 (MQ=255)
                                                                                        aCGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCg                                                               <  1:4159706/150‑1 (MQ=255)
                                                                                        aCGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCg                                                               >  2:3117226/1‑150 (MQ=255)
                                                                                            tgtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCa                                                           >  2:3049872/1‑150 (MQ=255)
                                                                                            tgtgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCa                                                           <  2:833820/150‑1 (MQ=255)
                                                                                              tgGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATa                                                         <  1:2479541/150‑1 (MQ=255)
                                                                                                 gTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGt                                                      <  2:3230924/150‑1 (MQ=255)
                                                                                                 gTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGt                                                      >  2:1273421/1‑150 (MQ=255)
                                                                                                      cAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaaaa                                                 >  1:1354603/1‑150 (MQ=255)
                                                                                                       aGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaaaaa                                                >  2:4083091/1‑150 (MQ=255)
                                                                                                       aGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTaaaaaa                                                >  1:4332395/1‑150 (MQ=255)
                                                                                                        gCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAAc                                               <  1:2291450/150‑1 (MQ=255)
                                                                                                         ccATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTCTTAATCTGGTTCGCTCATAAGTAAAAAACg                                              >  2:3509354/1‑150 (MQ=255)
                                                                                                         ccATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACg                                              >  1:1592941/1‑150 (MQ=255)
                                                                                                             aCCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCAc                                          >  2:2863901/1‑150 (MQ=255)
                                                                                                               cTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCt                                        >  1:1693746/1‑150 (MQ=255)
                                                                                                                 tAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTgg                                      >  1:100604/1‑150 (MQ=255)
                                                                                                                 tAATGGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTgg                                      <  2:543581/150‑1 (MQ=255)
                                                                                                                  aaTTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGt                                     <  2:4057296/150‑1 (MQ=255)
                                                                                                                    ttGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGGTCGCTCATAAGTAAAAAACGGCACCTGGTGc                                   >  2:1660430/1‑150 (MQ=255)
                                                                                                                     tGACGCTCACTTCCCGTTGGGTGGGATGACGCAGTCAGGAACGGGCGGTTTTTTGGGCCCCGAGTGGCTGGAGGGTTGGTGTGAAACTAAGTGGGTGTGTGTAGGGTATTAATCTGGTGCGCTCATAATTAAAAAACGGCACCTGGTGcc                                  <  1:3007265/150‑1 (MQ=255)
                                                                                                                      gACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCg                                 >  1:879844/1‑150 (MQ=255)
                                                                                                                      gACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCg                                 <  2:3168123/150‑1 (MQ=255)
                                                                                                                        cGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTgaagtt                               >  2:2749888/1‑145 (MQ=255)
                                                                                                                         gCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGATGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGttt                              <  2:1351534/150‑1 (MQ=255)
                                                                                                                          cTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGtttt                             <  1:2139992/150‑1 (MQ=255)
                                                                                                                            aaCTTACGGTTGGGTGGTATGAAGCAGTCATGAACGGGCCTTTATTTTGGCCCCGACTGGGTGGACGGTTGGTGTGAAACTACGTCGGTGTGTGTACGTTTTTAATCTGTTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGtttttt                           <  1:2713614/150‑1 (MQ=255)
                                                                                                                              cTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTg                         <  1:3592089/150‑1 (MQ=255)
                                                                                                                              cTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTg                         <  2:4214710/150‑1 (MQ=255)
                                                                                                                              cTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTg                         <  1:1769155/150‑1 (MQ=255)
                                                                                                                                tACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTCTTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTc                       >  2:2791413/1‑150 (MQ=255)
                                                                                                                                tACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTTTc                       >  1:2092853/1‑150 (MQ=255)
                                                                                                                                  ccGTTTGGGGGGATGACGCAGTCAGGACCGGGCCGTTATTTGGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTACGTCGGTGTGTGTCCGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTg                     <  1:420498/150‑1 (MQ=255)
                                                                                                                                   cGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGTTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGa                    <  1:3073032/150‑1 (MQ=255)
                                                                                                                                       tggtggGATGACGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACGGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGAATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACa                <  1:2775450/150‑1 (MQ=255)
                                                                                                                                        ggtggGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACaa               >  2:1743859/1‑150 (MQ=255)
                                                                                                                                         gtggGATGACCCGGTCAGGCCCGGCCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTTAACATATGGCTGTGTGTGTCCGGTATTAATCTGGTTCGCAAATAACTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAAt              <  1:2645958/150‑1 (MQ=255)
                                                                                                                                           gggATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCa            <  2:3768939/150‑1 (MQ=255)
                                                                                                                                           gggATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGGCGTAATCTGGTTCGCTCAGAAGTAAAAAACGGCACCTGGTGCCGTTTTTGTGGCTGAAACCATcc            >  1:2805240/1‑149 (MQ=255)
                                                                                                                                               tGAAGCAGTCAGGACCGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTa        <  1:2456432/150‑1 (MQ=255)
                                                                                                                                                 aaGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtc      >  1:2033982/1‑150 (MQ=255)
                                                                                                                                                  aGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtct     >  1:1738745/1‑150 (MQ=255)
                                                                                                                                                   gCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtctt    <  2:3858542/150‑1 (MQ=255)
                                                                                                                                                   gCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtctt    <  1:336560/150‑1 (MQ=255)
                                                                                                                                                     aGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTAtcttta  >  1:3975862/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGGTGTGAAACTAAGTCGGTGTGTGTACGGTATTAATCTGGTTCGCTCATAAGTAAAAAACGGCACCTGGTGCCGTTTTTTTGTCTGAAACAATCACTTATCTTTC  >  CP009273/1443047‑1443345

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