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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,443,0930CA28.8% 570.8 / 218.4 277L440I (CTC→ATC) feaBphenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (115/82);  new base A (34/46);  total (149/128)
Fisher's exact test for biased strand distribution p-value = 1.73e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 288 total aligned reads displayed.

GGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGG  >  CP009273/1442945‑1443242
                                                                                                                                                    |                                                                                                                                                     
ggtAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAAtct                                                                                                                                                      >  1:3256973/1‑150 (MQ=255)
   aaaTCCCGATGCTAAATTGCGCTTATCTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAAtatctc                                                                                                                                                   <  2:4411801/150‑1 (MQ=255)
     aTCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTccc                                                                                                                                                 <  1:4430070/150‑1 (MQ=255)
     aTCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTccc                                                                                                                                                 <  2:650043/150‑1 (MQ=255)
     aTCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTccc                                                                                                                                                 <  1:1193254/150‑1 (MQ=255)
      tCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGACTATGGCTTAACTGCCAGTGTCTGGACGCACAATCTCTCCCa                                                                                                                                                >  2:5071483/1‑150 (MQ=255)
      tCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCa                                                                                                                                                <  1:4919222/150‑1 (MQ=255)
       cccGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAg                                                                                                                                               >  1:3361770/1‑150 (MQ=255)
       cccGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAg                                                                                                                                               <  1:2576841/150‑1 (MQ=255)
            tGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGctct                                                                                                                                          >  1:1041362/1‑150 (MQ=255)
              cTAACTTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTgg                                                                                                                                        <  2:2830980/150‑1 (MQ=255)
                  aTTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAata                                                                                                                                    >  2:926944/1‑150 (MQ=255)
                   ttGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAatat                                                                                                                                   >  1:770527/1‑150 (MQ=255)
                   ttGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAatat                                                                                                                                   <  2:859699/150‑1 (MQ=255)
                    tGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAAtata                                                                                                                                  <  2:4803458/150‑1 (MQ=255)
                      cgcTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGc                                                                                                                                <  1:864039/150‑1 (MQ=255)
                      cgcTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGc                                                                                                                                <  1:2185689/150‑1 (MQ=255)
                       gcTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGACTATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCg                                                                                                                               >  2:3330742/1‑150 (MQ=255)
                        cTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGa                                                                                                                              <  1:2404579/150‑1 (MQ=255)
                        cTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGa                                                                                                                              >  2:2595724/1‑150 (MQ=255)
                        cTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGa                                                                                                                              >  1:599755/1‑150 (MQ=255)
                        cTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGa                                                                                                                              >  1:3578065/1‑150 (MQ=255)
                         ttAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCTTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGAt                                                                                                                             >  2:777322/1‑150 (MQ=255)
                           aaCTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCg                                                                                                                           <  1:4732497/150‑1 (MQ=255)
                              tCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCtt                                                                                                                        >  1:3177452/1‑150 (MQ=255)
                                  gatgagGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATATCGATCGCTTAcag                                                                                                                    >  2:3762771/4‑150 (MQ=255)
                                  gAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGACTATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATCGCGATCGCTTAcag                                                                                                                    >  2:4717723/1‑150 (MQ=255)
                                  gAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTAcag                                                                                                                    >  2:4587011/1‑150 (MQ=255)
                                  gAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTAcag                                                                                                                    >  2:3965845/1‑150 (MQ=255)
                                   aaGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTAcagg                                                                                                                   >  1:2474094/1‑150 (MQ=255)
                                    agagGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTAcaggc                                                                                                                  <  1:4780914/150‑1 (MQ=255)
                                     gagGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCTGGCTCTGGAATATAGCGATCGCTTAcaggca                                                                                                                 <  1:310526/150‑1 (MQ=255)
                                     gagGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTAcaggca                                                                                                                 >  1:4096144/1‑150 (MQ=255)
                                      agGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTAcaggcag                                                                                                                >  2:3909047/1‑150 (MQ=255)
                                      agGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTAcaggcag                                                                                                                >  2:4154725/1‑150 (MQ=255)
                                       ggTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTAcaggcagg                                                                                                               <  1:1564051/150‑1 (MQ=255)
                                         tgtTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGa                                                                                                             >  2:2390210/1‑150 (MQ=255)
                                         tgtTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGa                                                                                                             <  2:4570788/150‑1 (MQ=255)
                                          gtTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGAc                                                                                                            >  2:1813972/1‑150 (MQ=255)
                                             tGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGgt                                                                                                         >  2:3822070/1‑150 (MQ=255)
                                             tGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGgt                                                                                                         >  2:3116602/1‑150 (MQ=255)
                                               gTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGgtgt                                                                                                       <  1:2683866/150‑1 (MQ=255)
                                                tCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGtgtg                                                                                                      <  1:4091824/150‑1 (MQ=255)
                                                tCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGtgtg                                                                                                      <  2:5037104/150‑1 (MQ=255)
                                                tCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGtgtg                                                                                                      >  2:3384219/1‑150 (MQ=255)
                                                 ccGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTgg                                                                                                     >  1:4346592/1‑150 (MQ=255)
                                                   ggtggtAAACCTGGTGCGATAAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGAGTGGGt                                                                                                   >  2:4321401/1‑150 (MQ=255)
                                                   ggtggtAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCACAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGt                                                                                                   >  1:1068370/1‑150 (MQ=255)
                                                   ggtggtAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGt                                                                                                   >  1:4123864/1‑150 (MQ=255)
                                                   ggtggtAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGt                                                                                                   >  1:4695009/1‑150 (MQ=255)
                                                   ggtggtAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGt                                                                                                   >  1:3548110/1‑150 (MQ=255)
                                                   ggtggtAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGt                                                                                                   >  1:1923815/1‑150 (MQ=255)
                                                      ggtAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTaaa                                                                                                >  2:1915761/1‑150 (MQ=255)
                                                            ccTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGTTCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCa                                                                                          >  2:4337263/1‑150 (MQ=255)
                                                             cTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCAt                                                                                         >  1:3041247/1‑150 (MQ=255)
                                                              tGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATa                                                                                        <  2:2705750/150‑1 (MQ=255)
                                                                 tGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCt                                                                                     <  1:1781374/150‑1 (MQ=255)
                                                                  gCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCtt                                                                                    >  1:4471321/1‑150 (MQ=255)
                                                                     aGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAAt                                                                                 <  1:2371581/150‑1 (MQ=255)
                                                                             aTGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTa                                                                         >  1:3906146/1‑150 (MQ=255)
                                                                               ggAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAAc                                                                       <  2:1102130/150‑1 (MQ=255)
                                                                                    agagGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTAcc                                                                  >  1:651747/1‑150 (MQ=255)
                                                                                    agagGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTAcc                                                                  >  2:2241927/1‑150 (MQ=255)
                                                                                        gCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGttt                                                              <  2:1931191/150‑1 (MQ=255)
                                                                                        gCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGttt                                                              <  2:4457355/150‑1 (MQ=255)
                                                                                         cGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTtg                                                             >  2:1683376/1‑150 (MQ=255)
                                                                                            tACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTtggtg                                                          <  2:2430649/150‑1 (MQ=255)
                                                                                            tACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTtggtg                                                          >  2:2912690/1‑150 (MQ=255)
                                                                                              cAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTggg                                                        >  1:1395071/1‑150 (MQ=255)
                                                                                                   ggCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGaa                                                   >  1:2108830/1‑150 (MQ=255)
                                                                                                     cAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGc                                                 >  2:212867/1‑150 (MQ=255)
                                                                                                       aaCGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAg                                               <  1:3842318/150‑1 (MQ=255)
                                                                                                            cacGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAgg                                          >  2:567134/1‑150 (MQ=255)
                                                                                                               gtaaTATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAAc                                       <  1:1041463/148‑1 (MQ=255)
                                                                                                                   tatGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGc                                   <  2:4020805/150‑1 (MQ=255)
                                                                                                                    atGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGcc                                  <  1:2316237/150‑1 (MQ=255)
                                                                                                                      ggCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGt                                >  1:924806/1‑150 (MQ=255)
                                                                                                                       gCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTg                               >  2:4250344/1‑150 (MQ=255)
                                                                                                                        cTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGa                              >  2:1209088/1‑150 (MQ=255)
                                                                                                                        cTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGa                              >  2:1712301/1‑150 (MQ=255)
                                                                                                                         ttAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTCACTTACCGGTTGGTGGGATGAAGCTGTCAGGACCGGGCCGGGag                             >  2:3057053/1‑149 (MQ=255)
                                                                                                                            aCTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGAtttt                          >  1:3510784/1‑150 (MQ=255)
                                                                                                                             cTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCGTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTg                         <  1:3070224/150‑1 (MQ=255)
                                                                                                                               gCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGc                       <  1:4638680/150‑1 (MQ=255)
                                                                                                                                ccAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGcc                      <  1:5097821/150‑1 (MQ=255)
                                                                                                                                 cAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGccc                     <  1:3116602/150‑1 (MQ=255)
                                                                                                                                 cAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGccc                     >  2:1586014/1‑150 (MQ=255)
                                                                                                                                  aGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGcccc                    <  1:1978960/150‑1 (MQ=255)
                                                                                                                                   gtgtCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCg                   <  1:4163280/150‑1 (MQ=255)
                                                                                                                                      tCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGAct                <  2:1886929/150‑1 (MQ=255)
                                                                                                                                       cTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActg               >  2:3438751/1‑150 (MQ=255)
                                                                                                                                           aCGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActggctg           <  1:3996367/150‑1 (MQ=255)
                                                                                                                                           aCGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActggctg           <  2:3193292/150‑1 (MQ=255)
                                                                                                                                             gCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGa         <  1:2434646/150‑1 (MQ=255)
                                                                                                                                                aaaTCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACgg      >  1:548136/1‑150 (MQ=255)
                                                                                                                                                aaaTCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGGTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACgg      >  1:289749/1‑150 (MQ=255)
                                                                                                                                                  aTCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGtt    <  1:1683376/150‑1 (MQ=255)
                                                                                                                                                   tctctcCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTg   <  2:3543925/150‑1 (MQ=255)
                                                                                                                                                   tctctcCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTg   >  2:34895/1‑150 (MQ=255)
                                                                                                                                                    ctctcCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTATCTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTgg  <  1:2659945/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGG  >  CP009273/1442945‑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: