Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 29,139 G→T 100% G256C (GGT→TGT)  dapB → dihydrodipicolinate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP00927329,1390GT95.2% 954.2 / 32.8 304G256C (GGT→TGT) dapBdihydrodipicolinate reductase
Reads supporting (aligned to +/- strand):  ref base G (7/8);  new base T (105/184);  total (112/192)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 304 total aligned reads displayed.

GGTTTTGCCACCGTGCGTGCAGGTGACATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTGGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAATTTTTGGCCCTTTTATTTTTGGTGT  >  CP009273/28992‑29284
                                                                                                                                                   |                                                                                                                                                 
     tGCCACCGTGCGTGCAGGTGACATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGa                                                                                                                                            >  2:4533747/1‑150 (MQ=255)
     tGCCACCGTGCGTGCAGGTGACATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGa                                                                                                                                            >  2:2183309/1‑150 (MQ=255)
       ccACCGTGCGTGCAGGTGACATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGaaa                                                                                                                                          <  1:555526/150‑1 (MQ=255)
         cccGTGCGTGCAGGAGACATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGc                                                                                                                                        <  1:4413544/149‑1 (MQ=255)
           cgtgcgtgCAGGTGACATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCgg                                                                                                                                      >  1:3026680/1‑150 (MQ=255)
            gtgcgtgcAGGTGACATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGt                                                                                                                                     >  1:1034685/1‑150 (MQ=255)
            gtgcgtgcAGGTGACATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGt                                                                                                                                     >  1:2866450/1‑150 (MQ=255)
            gtgcgtgcAGGTGACATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGt                                                                                                                                     >  2:2651190/1‑150 (MQ=255)
             tgcgtgcAGGTGACATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCGATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTc                                                                                                                                    <  2:4396094/150‑1 (MQ=255)
             tgcgtgcAGGTGACATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTc                                                                                                                                    <  1:103805/150‑1 (MQ=255)
              gcgtgcAGGTGACATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCt                                                                                                                                   <  1:4532429/150‑1 (MQ=255)
               cgtgcAGGTGACATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCtt                                                                                                                                  <  1:3061467/150‑1 (MQ=255)
                 tgcAGGTGACATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCtttt                                                                                                                                >  1:1221220/1‑150 (MQ=255)
                    aGGTGACATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGa                                                                                                                             >  2:2617666/1‑150 (MQ=255)
                        gACATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATg                                                                                                                         <  1:503005/150‑1 (MQ=255)
                        gACATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATg                                                                                                                         <  1:2328024/150‑1 (MQ=255)
                         aCATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGTCATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGc                                                                                                                        <  1:2082286/150‑1 (MQ=255)
                           aTCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCga                                                                                                                      >  2:350354/1‑150 (MQ=255)
                           aTCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCga                                                                                                                      >  2:4490133/1‑150 (MQ=255)
                             cGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCgaga                                                                                                                    <  2:615229/150‑1 (MQ=255)
                             cGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCgaga                                                                                                                    <  2:614742/150‑1 (MQ=255)
                               ttGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATg                                                                                                                  <  1:4272303/150‑1 (MQ=255)
                               ttGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATg                                                                                                                  <  2:1444847/150‑1 (MQ=255)
                               ttGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATg                                                                                                                  <  2:1416120/150‑1 (MQ=255)
                                 ggTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCGCCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTa                                                                                                                <  2:2408044/150‑1 (MQ=255)
                                 ggTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTa                                                                                                                >  1:4299983/1‑150 (MQ=255)
                                 ggTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTa                                                                                                                >  2:3591497/1‑150 (MQ=255)
                                   tGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACt                                                                                                              >  1:3622937/1‑150 (MQ=255)
                                    gAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACtt                                                                                                             <  2:2486933/150‑1 (MQ=255)
                                    gAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACtt                                                                                                             <  2:4205575/150‑1 (MQ=255)
                                       cATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGAt                                                                                                          <  1:3672944/150‑1 (MQ=255)
                                       cATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGAt                                                                                                          <  2:2368077/150‑1 (MQ=255)
                                              cgATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCaata                                                                                                   <  2:3319783/150‑1 (MQ=255)
                                              cgATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCaata                                                                                                   >  1:4043544/1‑150 (MQ=255)
                                              cgATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCaata                                                                                                   >  1:3052051/1‑150 (MQ=255)
                                              cgATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCaata                                                                                                   >  1:1977501/1‑150 (MQ=255)
                                               gATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCaataa                                                                                                  >  1:1353061/1‑150 (MQ=255)
                                                aTGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCaataat                                                                                                 >  2:3927442/1‑150 (MQ=255)
                                                 tGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAAtt                                                                                                >  2:1416582/1‑150 (MQ=255)
                                                      gCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTGGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTaa                                                                                           <  1:1006720/150‑1 (MQ=255)
                                                        cGATTTTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAAca                                                                                         <  2:1597737/150‑2 (MQ=255)
                                                               ggCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAat                                                                                  >  2:415861/1‑150 (MQ=255)
                                                                  gAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATAtt                                                                               >  1:706995/1‑150 (MQ=255)
                                                                  gAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATAtt                                                                               <  1:4180910/150‑1 (MQ=255)
                                                                   aGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATAttt                                                                              <  2:1284385/150‑1 (MQ=255)
                                                                      gTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGtt                                                                           >  2:3806522/1‑150 (MQ=255)
                                                                       tCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACGTTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTa                                                                          >  2:1795295/1‑150 (MQ=255)
                                                                        cTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTGGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTAt                                                                         <  2:2998904/150‑1 (MQ=255)
                                                                         tGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATg                                                                        >  1:625797/1‑150 (MQ=255)
                                                                           gagaTCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGt                                                                      <  2:535953/150‑1 (MQ=255)
                                                                           gagaTCACCCAGATGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGt                                                                      <  2:3675958/150‑1 (MQ=255)
                                                                             gaTCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGc                                                                    >  2:2766716/1‑150 (MQ=255)
                                                                              aTCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCa                                                                   >  1:4208145/1‑150 (MQ=255)
                                                                                     aTAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATaa                                                            >  1:1178891/1‑150 (MQ=255)
                                                                                      tAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTATGATCGGCTTTGTGGTTGTGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCTATACTTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAac                                                           <  2:715035/150‑1 (MQ=255)
                                                                                      tAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAac                                                           >  1:1988347/1‑150 (MQ=255)
                                                                                       aaGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAaca                                                          >  1:1754861/1‑150 (MQ=255)
                                                                                             tCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACattta                                                    >  2:2085971/1‑150 (MQ=255)
                                                                                             tCCAGCAGTATGACATTTGCTAACGGCGCGGTAAGATCTGCTTTGTGTTTGAGTTGTAAGGAAAGCTGTCTTTTTGATATGCGAGATGTACTTGATCTAAATAATTTGTAACCACAAAATATTTGTTATGGTGCTAAAATAACACattta                                                    <  2:4280854/150‑1 (MQ=255)
                                                                                              ccAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACatttaa                                                   >  1:4497045/1‑150 (MQ=255)
                                                                                                aGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACatttaatt                                                 <  2:3131000/150‑1 (MQ=255)
                                                                                                aGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACatttaatt                                                 >  2:2672593/1‑150 (MQ=255)
                                                                                                 gCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACatttaattt                                                >  1:3310/1‑150 (MQ=255)
                                                                                                 gCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACatttaattt                                                >  1:3847/1‑150 (MQ=255)
                                                                                                 gCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACatttaattt                                                <  1:163686/150‑1 (MQ=255)
                                                                                                  ccGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACatttaattta                                               <  2:1564860/150‑1 (MQ=255)
                                                                                                  ccGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTGGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACatttaattta                                               <  1:4330351/150‑1 (MQ=255)
                                                                                                    gTATGACATTTGCTACCGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTAtt                                             <  2:3488384/150‑1 (MQ=255)
                                                                                                    gTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTAtt                                             <  2:3488161/150‑1 (MQ=255)
                                                                                                      aTGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGa                                           >  2:4048522/1‑150 (MQ=255)
                                                                                                      aTGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGa                                           >  2:1086357/1‑150 (MQ=255)
                                                                                                       tGACATTTGCTTACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGATAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGAt                                          <  2:85928/150‑1 (MQ=255)
                                                                                                       tGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGAt                                          <  1:3292439/150‑1 (MQ=255)
                                                                                                        gACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGAtt                                         >  2:1876528/1‑150 (MQ=255)
                                                                                                        gACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGAtt                                         >  1:2129941/1‑150 (MQ=255)
                                                                                                          cATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATtat                                       >  2:1517874/1‑150 (MQ=255)
                                                                                                           aTTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTGGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATtata                                      >  1:4184673/1‑150 (MQ=255)
                                                                                                             ttGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATaaa                                    <  1:2282209/150‑1 (MQ=255)
                                                                                                             ttGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATaaa                                    <  1:2101321/150‑1 (MQ=255)
                                                                                                                cTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAggg                                 >  2:1370021/1‑150 (MQ=255)
                                                                                                                 tAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGc                                >  1:3562549/1‑150 (MQ=255)
                                                                                                                   aCGGCGCGGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCtt                              <  1:1607817/150‑1 (MQ=255)
                                                                                                                       cgcgGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAAt                          <  2:1213881/150‑1 (MQ=255)
                                                                                                                        gcgGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAAtt                         <  2:367426/150‑1 (MQ=255)
                                                                                                                         cgGTAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAAttt                        <  1:90212/150‑1 (MQ=255)
                                                                                                                            tAAGATCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAATTTTTg                     <  1:2609163/150‑1 (MQ=255)
                                                                                                                             acgttCGGCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCACAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAATTTTTgg                    >  2:4300403/5‑150 (MQ=255)
                                                                                                                                    gCTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAATTTTTGGCCCtttt             <  2:2036102/150‑1 (MQ=255)
                                                                                                                                     cTTTGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAATTTTTGGCCCTTTTa            >  1:4072000/1‑150 (MQ=255)
                                                                                                                                      tttGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAATTTTTGGCCCTTTTAt           >  1:1060679/1‑150 (MQ=255)
                                                                                                                                       ttGTGGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAATTTTTGGCCCTTTTAtt          <  1:3277180/150‑1 (MQ=255)
                                                                                                                                        tgtttttGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAATTTTTGGCCCTTTTAttt         <  2:3695998/145‑1 (MQ=255)
                                                                                                                                         gtgGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAATTTTTGGCCCTTTTAtttt        <  1:2622267/150‑1 (MQ=255)
                                                                                                                                         gtgGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAATTTTTGGCCCTTTTAtttt        <  1:1643003/150‑1 (MQ=255)
                                                                                                                                          tgGTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAATTTTTGGCCCTTTTAttttt       >  1:3748784/1‑150 (MQ=255)
                                                                                                                                           ggTTGAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAATTTTTGGCCCTTTTATTTTTg      >  2:694379/1‑150 (MQ=255)
                                                                                                                                              tGAGTTGTAAGGAAAGCGGTCATTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAATTTTTGGCCCTTTTATTTTTGgtg   <  2:3134630/150‑1 (MQ=255)
                                                                                                                                               gAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACCTTTAATTTATTGATTATAAAGGGCTTTAATTTTTGGCCCTTTTATTTTTGgtgt  <  2:74168/150‑1 (MQ=255)
                                                                                                                                               gAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAATTTTTGGCCCTTTTATTTTTGgtgt  <  1:3017284/150‑1 (MQ=255)
                                                                                                                                               gAGTTGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAATTTTTGGCCCTTTTATTTTTGgtgt  <  2:1034685/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GGTTTTGCCACCGTGCGTGCAGGTGACATCGTTGGTGAACATACCGCGATGTTTGCCGATATTGGCGAGCGTCTGGAGATCACCCATAAGGCGTCCAGCCGTATGACATTTGCTAACGGCGCGGTAAGATCGGCTTTGTGGTTGAGTGGTAAGGAAAGCGGTCTTTTTGATATGCGAGATGTACTTGATCTCAATAATTTGTAACCACAAAATATTTGTTATGGTGCAAAAATAACACATTTAATTTATTGATTATAAAGGGCTTTAATTTTTGGCCCTTTTATTTTTGGTGT  >  CP009273/28992‑29284

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: