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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,443,0930CA12.6% 777.4 / 63.0 235L440I (CTC→ATC) feaBphenylacetaldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (85/120);  new base A (15/15);  total (100/135)
Fisher's exact test for biased strand distribution p-value = 4.31e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

Only 100 of 245 total aligned reads displayed.

GGTAAATCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGTTGG  >  CP009273/1442945‑1443242
                                                                                                                                                    |                                                                                                                                                     
   aaaTCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAAtctctc                                                                                                                                                   <  1:1723216/150‑1 (MQ=255)
     aTCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTccc                                                                                                                                                 >  1:1338521/1‑150 (MQ=255)
      tCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCa                                                                                                                                                >  1:601091/1‑150 (MQ=255)
      tCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCa                                                                                                                                                >  2:81417/1‑150 (MQ=255)
      tCCCGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCa                                                                                                                                                <  1:4060893/150‑1 (MQ=255)
         cGATGCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGc                                                                                                                                             >  1:1431050/1‑150 (MQ=255)
         cGATGCTAAATTGCGCTTAACTCGAGAAGAGGTTTTAGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACTACACGGAATATGGCTTAACTGCCATTGTCTGGACGCAAAATCTCTCCCAGGc                                                                                                                                             <  2:1819145/150‑1 (MQ=255)
             gCTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTg                                                                                                                                         >  1:4505115/1‑150 (MQ=255)
              cTAAATTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTgg                                                                                                                                        <  2:675127/150‑1 (MQ=255)
                aaaTTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGaa                                                                                                                                      <  1:1567795/150‑1 (MQ=255)
                 aaTTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAat                                                                                                                                     <  2:4114526/150‑1 (MQ=255)
                 aaTTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAat                                                                                                                                     <  1:2418476/150‑1 (MQ=255)
                  aTTGCGCTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAata                                                                                                                                    <  2:1012850/150‑1 (MQ=255)
                       gcTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCg                                                                                                                               >  2:583540/1‑150 (MQ=255)
                        cTTAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGa                                                                                                                              <  2:2678792/150‑1 (MQ=255)
                          tAACTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATc                                                                                                                            <  2:3694481/150‑1 (MQ=255)
                            aCTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGc                                                                                                                          <  2:698152/150‑1 (MQ=255)
                             cTCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCt                                                                                                                         >  2:1350365/1‑150 (MQ=255)
                              tCGTGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCtt                                                                                                                        <  2:3309987/150‑1 (MQ=255)
                                 tGAAGAGTTGTTTGGTCCGGTGGTTAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCATACGACACGGAATATTGCTTAACTGCCAGTGTCTGGAAGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTAca                                                                                                                     <  2:1186619/150‑1 (MQ=255)
                                 tGAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTAca                                                                                                                     <  2:4290574/150‑1 (MQ=255)
                                  gAAGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTAcag                                                                                                                    >  1:2066780/1‑150 (MQ=255)
                                   aaGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTAcagg                                                                                                                   <  1:1777690/150‑1 (MQ=255)
                                   aaGAGGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTAcagg                                                                                                                   <  1:1490815/150‑1 (MQ=255)
                                    agagGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTAcaggc                                                                                                                  <  1:169265/150‑1 (MQ=255)
                                    agagGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTAcaggc                                                                                                                  >  1:1397568/1‑150 (MQ=255)
                                      agGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTAcaggcag                                                                                                                >  2:165276/1‑150 (MQ=255)
                                      agGTGTTTGGTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTAcaggcag                                                                                                                <  2:2639928/150‑1 (MQ=255)
                                              ggTCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGgtg                                                                                                        <  2:489987/150‑1 (MQ=255)
                                                tCCGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGtgtg                                                                                                      <  1:1283305/150‑1 (MQ=255)
                                                 ccGGTGGTTAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTgg                                                                                                     <  2:1694646/150‑1 (MQ=255)
                                                 ccGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTgg                                                                                                     >  1:3610955/1‑150 (MQ=255)
                                                  cGGTGGTAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTggg                                                                                                    <  1:1744799/150‑1 (MQ=255)
                                                    gtggtAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCTCTTACAGGCAGGGACGGTGTGGGTa                                                                                                  >  1:884001/1‑150 (MQ=255)
                                                     tggtAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTaa                                                                                                 <  2:3719514/150‑1 (MQ=255)
                                                        tAAACCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACa                                                                                              <  1:3467112/150‑1 (MQ=255)
                                                          aaCCTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGc                                                                                            >  2:3588668/1‑150 (MQ=255)
                                                            ccTGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCa                                                                                          >  2:886902/1‑150 (MQ=255)
                                                              tGGTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATa                                                                                        <  1:3624975/150‑1 (MQ=255)
                                                                gTGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATAcc                                                                                      <  2:3379739/150‑1 (MQ=255)
                                                                 tGCGAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCt                                                                                     <  1:2209839/150‑1 (MQ=255)
                                                                    taGTAGCGGATGGAGAATAGGCGTTACAAGTGGAAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCATGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTaa                                                                                  >  2:4403846/2‑150 (MQ=255)
                                                                    gAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTaa                                                                                  <  2:4211538/150‑1 (MQ=255)
                                                                    gAGTAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTaa                                                                                  <  2:4211615/150‑1 (MQ=255)
                                                                     aGTAGCGGATGGAGAAGATGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAAt                                                                                 >  1:2054647/1‑150 (MQ=255)
                                                                       tAGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTg                                                                               >  2:1610114/1‑150 (MQ=255)
                                                                        aGCGGATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGa                                                                              <  1:1084580/150‑1 (MQ=255)
                                                                           ggATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGc                                                                           >  2:2695546/1‑150 (MQ=255)
                                                                            gATGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCt                                                                          <  2:4017958/150‑1 (MQ=255)
                                                                              tGGAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTaa                                                                        <  2:4366393/150‑1 (MQ=255)
                                                                               ggAGAAGAGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAAc                                                                       >  2:2387544/1‑150 (MQ=255)
                                                                                 agaagaGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCCGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTCAACAGCCATCCCTTAATTGACGCTAACtt                                                                     >  2:3838467/1‑150 (MQ=255)
                                                                                 agaagaGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACtt                                                                     >  2:3836824/1‑150 (MQ=255)
                                                                                   aagaGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTAc                                                                   <  2:523931/150‑1 (MQ=255)
                                                                                   aagaGGCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTAc                                                                   <  2:4407573/150‑1 (MQ=255)
                                                                                        gCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGttt                                                              <  1:1569368/150‑1 (MQ=255)
                                                                                        gCGTTACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGttt                                                              <  1:28296/150‑1 (MQ=255)
                                                                                           ttACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTtggt                                                           <  1:173336/150‑1 (MQ=255)
                                                                                            tACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTtggtg                                                          <  2:2022747/150‑1 (MQ=255)
                                                                                            tACAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTtggtg                                                          <  2:2420413/150‑1 (MQ=255)
                                                                                             aCAACTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTtggtgg                                                         <  1:1521361/150‑1 (MQ=255)
                                                                                               aaCTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGa                                                       >  2:379129/1‑150 (MQ=255)
                                                                                                aCTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGAt                                                      >  1:3732605/1‑150 (MQ=255)
                                                                                                 cTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATg                                                     <  2:86349/150‑1 (MQ=255)
                                                                                                 cTGGCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATg                                                     >  2:3342439/1‑150 (MQ=255)
                                                                                                   ggCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGaa                                                   >  1:2256966/1‑150 (MQ=255)
                                                                                                    gCAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAg                                                  >  1:330144/1‑150 (MQ=255)
                                                                                                     cAAACGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGc                                                 >  1:2951407/1‑150 (MQ=255)
                                                                                                      aaaCGACCCGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCGTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCa                                                <  1:3482308/150‑1 (MQ=255)
                                                                                                      aaaCGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCa                                                <  1:3482275/150‑1 (MQ=255)
                                                                                                      aaaCGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCa                                                <  1:1235833/150‑1 (MQ=255)
                                                                                                       aaCGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAg                                               <  1:1868502/150‑1 (MQ=255)
                                                                                                       aaCGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGctg                                               <  1:4114531/150‑3 (MQ=255)
                                                                                                        aCGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTGTGGTGGGATGAAGCAGt                                              >  2:1497836/1‑150 (MQ=255)
                                                                                                         cGACACGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTc                                             <  1:426495/150‑1 (MQ=255)
                                                                                                          gACACGGAATATGGCTTAACTGCCTGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCa                                            <  2:3422641/150‑1 (MQ=255)
                                                                                                             acGGAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGa                                         >  2:3153966/1‑150 (MQ=255)
                                                                                                                gAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACg                                      <  2:642013/150‑1 (MQ=255)
                                                                                                                gAATATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACg                                      <  2:3310842/150‑1 (MQ=255)
                                                                                                                 aaTATGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACgg                                     <  1:519026/150‑1 (MQ=255)
                                                                                                                   tatGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTGTGGTGGGATGAAGCAGTCAGGAACGGGc                                   >  1:1084211/1‑150 (MQ=255)
                                                                                                                   tatGGCTTAACTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACATTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGc                                   >  2:3582263/1‑150 (MQ=255)
                                                                                                                    atGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGcc                                  >  2:2201811/1‑150 (MQ=255)
                                                                                                                    atGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGcc                                  >  2:2201179/1‑150 (MQ=255)
                                                                                                                     tGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCg                                 >  2:280437/1‑150 (MQ=255)
                                                                                                                     tGGCTTAACTGCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCg                                 >  1:4146508/1‑150 (MQ=255)
                                                                                                                           aaCTGCCAGTGTCTGGACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGAttt                           >  1:991389/1‑150 (MQ=255)
                                                                                                                             ccgCCCGTGTCTGGACCCAACACCTCACCCACGCTCTCGAAGAGAGCGATCGCTTACCGGCAGGGACGGTGTGCGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTg                         <  2:2587215/148‑1 (MQ=255)
                                                                                                                               gCCAGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGc                       >  1:4334540/1‑150 (MQ=255)
                                                                                                                                  aGTGTCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGTCAGGTACGGTGTGGGTAAACAGCCATACCTAAATTGACGCTAACTTACCGTTTGGTGGGATGAAGAAGTCAGTAACGGGCCGTGATTTTGGcacc                    >  2:2126984/1‑150 (MQ=255)
                                                                                                                                      tCTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGAct                >  1:1716158/1‑150 (MQ=255)
                                                                                                                                       cTGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActg               >  2:3215564/1‑150 (MQ=255)
                                                                                                                                        tGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActgg              >  1:3207664/1‑150 (MQ=255)
                                                                                                                                        tGGACGCAAAATCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGActgg              <  1:81956/150‑1 (MQ=255)
                                                                                                                                        tGGACGCAAAATATCTCCCAGGCTCTGGAATATCGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGActgg              >  2:1460507/1‑150 (MQ=255)
                                                                                                                                         ggACGCAAAATATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGActggc             <  2:1522930/150‑1 (MQ=255)
                                                                                                                                                aaaTCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCCGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACgg      >  2:3246052/1‑150 (MQ=255)
                                                                                                                                                aaaTATCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCTGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACgg      >  1:2207353/1‑150 (MQ=255)
                                                                                                                                                 aaTCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGt     <  1:2665046/150‑1 (MQ=255)
                                                                                                                                                 aaTCTCTCCCAGGCTCTGGAATATAGCGATCGCTTACAGGCAGGGACGGTGTGGGTAAACAGCCATACCTTAATTGACGCTAACTTACCGTTTGGTGGGATGAAGCAGTCAGGAACGGGCCGTGATTTTGGCCCCGACTGGCTGGACGGt     <  2:3490424/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: