Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,443,093 C→A 15.9% L440I (CTC→ATC)  feaB → phenylacetaldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,443,0930CA15.9% 765.2 / 96.1 250L440I (CTC→ATC) feaBphenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (105/105);  new base A (19/21);  total (124/126)
Fisher's exact test for biased strand distribution p-value = 8.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 265 total aligned reads displayed.

TGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGG  >  CP009273/1442944‑1443242
                                                                                                                                                     |                                                                                                                                                     
 ggtAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAAtct                                                                                                                                                      >  2:1455089/1‑150 (MQ=255)
 ggtAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAAtct                                                                                                                                                      >  1:3441779/1‑150 (MQ=255)
      aTCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTccc                                                                                                                                                 <  1:532050/150‑1 (MQ=255)
       tCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCa                                                                                                                                                <  2:2798019/150‑1 (MQ=255)
       tCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCa                                                                                                                                                <  1:420397/150‑1 (MQ=255)
        cccGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAg                                                                                                                                               <  1:4313711/150‑1 (MQ=255)
        cccGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAg                                                                                                                                               <  2:3065563/150‑1 (MQ=255)
        cccGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAg                                                                                                                                               <  1:1972485/150‑1 (MQ=255)
         ccGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAgg                                                                                                                                              <  2:1444067/150‑1 (MQ=255)
         ccGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAgg                                                                                                                                              <  2:59838/150‑1 (MQ=255)
           gATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGct                                                                                                                                            <  2:3800446/150‑1 (MQ=255)
             tGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGctct                                                                                                                                          >  2:1554461/1‑150 (MQ=255)
             tGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGctct                                                                                                                                          >  2:2071177/1‑150 (MQ=255)
              gCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTg                                                                                                                                         <  1:2351199/150‑1 (MQ=255)
               cTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTgg                                                                                                                                        >  1:3335502/1‑150 (MQ=255)
                tAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGa                                                                                                                                       >  1:822681/1‑150 (MQ=255)
                 aaaTTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGaa                                                                                                                                      <  1:4285913/150‑1 (MQ=255)
                  aaTTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAat                                                                                                                                     >  1:977501/1‑150 (MQ=255)
                      gcgcTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAg                                                                                                                                 <  2:311624/150‑1 (MQ=255)
                        gcTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCg                                                                                                                               <  1:1294827/150‑1 (MQ=255)
                          ttAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGAt                                                                                                                             <  1:1317264/150‑1 (MQ=255)
                          ttAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGAt                                                                                                                             >  2:574765/1‑150 (MQ=255)
                           tAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATc                                                                                                                            >  2:2578658/1‑150 (MQ=255)
                           tAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATc                                                                                                                            <  1:1920347/150‑1 (MQ=255)
                            aaCTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCCGGCTCTGGAATATAGCGATCg                                                                                                                           >  2:3054137/1‑150 (MQ=255)
                              cTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCt                                                                                                                         >  1:2169294/1‑150 (MQ=255)
                                 gTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTAc                                                                                                                      >  2:2456432/1‑150 (MQ=255)
                                 gTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTAc                                                                                                                      >  1:783077/1‑150 (MQ=255)
                                   gAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTAcag                                                                                                                    <  2:710924/150‑1 (MQ=255)
                                   gAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTAcag                                                                                                                    <  1:1792528/150‑1 (MQ=255)
                                   gAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTAcag                                                                                                                    <  2:4195463/150‑1 (MQ=255)
                                           gtTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGAc                                                                                                            <  1:3461537/150‑1 (MQ=255)
                                            gttGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACg                                                                                                           >  1:1264893/2‑150 (MQ=255)
                                             ttGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACgg                                                                                                          >  1:286838/1‑150 (MQ=255)
                                             ttGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACgg                                                                                                          <  2:3465831/150‑1 (MQ=255)
                                              tGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGgt                                                                                                         <  2:4063035/150‑1 (MQ=255)
                                                gTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGgtgt                                                                                                       >  2:1982678/1‑150 (MQ=255)
                                                 tCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGtgtg                                                                                                      <  1:590087/150‑1 (MQ=255)
                                                   cGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTggg                                                                                                    >  1:3513376/1‑150 (MQ=255)
                                                    ggtggtAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGt                                                                                                   >  2:3604700/1‑150 (MQ=255)
                                                      tggtAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTaa                                                                                                 <  1:1114479/150‑1 (MQ=255)
                                                      tggtAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTaa                                                                                                 <  2:855829/150‑1 (MQ=255)
                                                      tggtAAACCTGGTGCGAGTAGAGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTATCTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACTGGCAGGGTCGGTGTGGGtga                                                                                                 >  2:2645958/1‑148 (MQ=255)
                                                          aaaCCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAg                                                                                             >  2:3921712/1‑150 (MQ=255)
                                                          aaaCCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAg                                                                                             >  2:1042319/1‑150 (MQ=255)
                                                          aaaCCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAg                                                                                             <  2:1393779/150‑1 (MQ=255)
                                                          aaaCCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAg                                                                                             <  2:1393345/150‑1 (MQ=255)
                                                               tGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATa                                                                                        <  2:3563477/150‑1 (MQ=255)
                                                                ggTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATAc                                                                                       >  1:3858542/1‑150 (MQ=255)
                                                                 gTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATAcc                                                                                      <  2:3286789/150‑1 (MQ=255)
                                                                  tGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCt                                                                                     <  2:3988048/150‑1 (MQ=255)
                                                                  tGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCt                                                                                     >  1:2421037/1‑150 (MQ=255)
                                                                      aGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAAt                                                                                 <  2:800479/150‑1 (MQ=255)
                                                                       gTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAAtt                                                                                >  2:3956282/1‑150 (MQ=255)
                                                                            ggATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGc                                                                           >  1:3769026/1‑150 (MQ=255)
                                                                                 gagaAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACt                                                                      >  2:811827/1‑150 (MQ=255)
                                                                                  agaagaGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACtt                                                                     <  2:2985979/150‑1 (MQ=255)
                                                                                    aagaGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTAc                                                                   >  1:1856697/1‑150 (MQ=255)
                                                                                    aagaGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTAc                                                                   >  1:3362977/1‑150 (MQ=255)
                                                                                      gagGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCg                                                                 <  1:2616223/150‑1 (MQ=255)
                                                                                       agGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGt                                                                >  2:372769/1‑150 (MQ=255)
                                                                                        ggCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGtt                                                               <  1:322617/150‑1 (MQ=255)
                                                                                        ggCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCCGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGtt                                                               >  2:2875005/1‑150 (MQ=255)
                                                                                         gCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGttt                                                              <  1:2237677/150‑1 (MQ=255)
                                                                                          cGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTtg                                                             <  1:563159/150‑1 (MQ=255)
                                                                                           gTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTGAATTGACGCTAACTTACCGTTtgg                                                            >  2:4274636/1‑150 (MQ=255)
                                                                                                 aCTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGAt                                                      <  1:3165422/150‑1 (MQ=255)
                                                                                                      cAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGc                                                 <  1:862325/150‑1 (MQ=255)
                                                                                                      cAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGc                                                 <  1:1572347/150‑1 (MQ=255)
                                                                                                       aaaCGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCa                                                >  1:1243310/1‑150 (MQ=255)
                                                                                                       aaaCGACACGGAATATGGCTTAACTGCCAGTGTCGGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGct                                                >  1:4109591/1‑149 (MQ=255)
                                                                                                        aaCGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTCATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAg                                               >  2:4133537/1‑150 (MQ=255)
                                                                                                        aaCGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGctg                                               <  2:3420288/150‑3 (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)
                                                                                                             cacGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAgg                                          >  2:639828/1‑150 (MQ=255)
                                                                                                                ggAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAAc                                       >  2:1191404/1‑150 (MQ=255)
                                                                                                                 gAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACg                                      >  1:4146172/1‑150 (MQ=255)
                                                                                                                  aaTATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACgg                                     <  1:3554708/150‑1 (MQ=255)
                                                                                                                  aaTATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACgg                                     <  1:1361433/150‑1 (MQ=255)
                                                                                                                   atatGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACggg                                    >  1:2768507/1‑150 (MQ=255)
                                                                                                                     atGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGcc                                  <  2:674894/150‑1 (MQ=255)
                                                                                                                       ggCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGt                                >  2:1806520/1‑150 (MQ=255)
                                                                                                                         cTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGa                              <  1:3983978/150‑1 (MQ=255)
                                                                                                                            aaCTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGAttt                           >  2:1769155/1‑150 (MQ=255)
                                                                                                                             aCTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGAtttt                          >  1:2991631/1‑150 (MQ=255)
                                                                                                                                gCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGc                       <  1:1453643/150‑1 (MQ=255)
                                                                                                                                  cAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGACGCAGTCAGGAACGGGCCGTGATTTTGGccc                     >  1:1035923/1‑150 (MQ=255)
                                                                                                                                   aGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGcccc                    <  1:1345812/150‑1 (MQ=255)
                                                                                                                                    gtgtCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCg                   <  1:2123251/150‑1 (MQ=255)
                                                                                                                                     tgtCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGa                  >  2:4012334/1‑150 (MQ=255)
                                                                                                                                     tgtCTGGACGCAAAATCACAACCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAAAAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAATCAGTCAGGAACGGTACTTGATTTTTGCCCCGa                  >  2:3262418/1‑150 (MQ=255)
                                                                                                                                               cAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGAc        >  1:3691964/1‑150 (MQ=255)
                                                                                                                                               cAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGAc        >  2:167341/1‑150 (MQ=255)
                                                                                                                                                aaaaTATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACg       >  2:111505/1‑150 (MQ=255)
                                                                                                                                                   aTCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGtt    >  1:4057296/1‑150 (MQ=255)
                                                                                                                                                     ctctcCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTgg  >  1:3683166/1‑150 (MQ=255)
                                                                                                                                                     ctctcCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTgg  <  2:225506/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGG  >  CP009273/1442944‑1443242

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: