Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,234,574 A→T 5.2% E80D (GAA→GAT dadX → alanine racemase, catabolic, PLP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,234,5740AT 5.2% 1027.0 / 34.3 288E80D (GAA→GATdadXalanine racemase, catabolic, PLP‑binding
Reads supporting (aligned to +/- strand):  ref base A (124/149);  new base T (8/7);  total (132/156)
Fisher's exact test for biased strand distribution p-value = 6.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.45e-01

Only 100 of 295 total aligned reads displayed.

TGGTCGGTGGTAAAAGCGAACGCTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCG  >  CP009273/1234425‑1234723
                                                                                                                                                     |                                                                                                                                                     
   tCGGTGGTAAAAGCGAACGCTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGa                                                                                                                                                    <  2:1583482/150‑1 (MQ=255)
           aaaaGCGAACGCTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCa                                                                                                                                            <  1:4634862/150‑1 (MQ=255)
           aaaaGCGAACGCTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCa                                                                                                                                            >  1:311365/1‑150 (MQ=255)
            aaaGCGAACGCTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCAt                                                                                                                                           <  2:928533/150‑1 (MQ=255)
            aaaGCGAACGCTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCAt                                                                                                                                           >  1:4158804/1‑150 (MQ=255)
            aaaGCGAACGCTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCAt                                                                                                                                           <  2:2137480/150‑1 (MQ=255)
             aaGCGAACGCTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATg                                                                                                                                          <  1:383760/150‑1 (MQ=255)
              aGCGAACGCTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGc                                                                                                                                         >  2:3417170/1‑150 (MQ=255)
               gCGAACGCTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCt                                                                                                                                        >  1:2503396/1‑150 (MQ=255)
                cGAACGCTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTc                                                                                                                                       >  1:2492719/1‑150 (MQ=255)
                cGAACGCTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTc                                                                                                                                       <  2:3155913/150‑1 (MQ=255)
                    cGCTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGa                                                                                                                                   <  2:2951041/150‑1 (MQ=255)
                      cTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATc                                                                                                                                 <  1:3922642/150‑1 (MQ=255)
                      cTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATc                                                                                                                                 <  1:746611/150‑1 (MQ=255)
                      cTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATc                                                                                                                                 >  1:3646485/1‑150 (MQ=255)
                      cTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATc                                                                                                                                 >  2:5036030/1‑150 (MQ=255)
                      cTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATc                                                                                                                                 <  2:3624767/150‑1 (MQ=255)
                       ttACGGGCATGGTATAGAGCGTCTGTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTTCGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCt                                                                                                                                <  2:3225446/150‑1 (MQ=255)
                          cGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGga                                                                                                                             >  1:53592/1‑150 (MQ=255)
                          cGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGga                                                                                                                             >  2:5218541/1‑150 (MQ=255)
                             gCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGAt                                                                                                                          <  2:53087/150‑1 (MQ=255)
                             gCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGAt                                                                                                                          <  1:2981087/150‑1 (MQ=255)
                              cATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGAtt                                                                                                                         <  1:3744207/150‑1 (MQ=255)
                               aTGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGATGGATTTTTCCATGCTCAGGATCTGGAGAttt                                                                                                                        <  1:708289/150‑1 (MQ=255)
                               aTGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGAttt                                                                                                                        <  2:4307983/150‑1 (MQ=255)
                                tGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTa                                                                                                                       >  1:3452016/1‑150 (MQ=255)
                                tGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTa                                                                                                                       <  1:4510746/150‑1 (MQ=255)
                                  gTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATg                                                                                                                     >  1:26702/1‑150 (MQ=255)
                                      tGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCa                                                                                                                 <  2:236827/150‑1 (MQ=255)
                                       gAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAg                                                                                                                >  1:614355/1‑150 (MQ=255)
                                        aGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGc                                                                                                               >  2:1861127/1‑150 (MQ=255)
                                        aGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGc                                                                                                               <  1:2645861/150‑1 (MQ=255)
                                            tctCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCg                                                                                                           <  2:4885475/148‑1 (MQ=255)
                                                tGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTg                                                                                                       <  2:4141299/150‑1 (MQ=255)
                                                tGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTg                                                                                                       >  2:4789332/1‑150 (MQ=255)
                                                tGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTg                                                                                                       <  1:730276/150‑1 (MQ=255)
                                                tGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTg                                                                                                       <  1:3355582/150‑1 (MQ=255)
                                                    gcgcGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGacca                                                                                                   <  2:1670349/150‑1 (MQ=255)
                                                    gcgcGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGacca                                                                                                   <  2:2065332/150‑1 (MQ=255)
                                                     cgcgATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGaccac                                                                                                  <  1:3692030/150‑1 (MQ=255)
                                                                ccACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTAcaca                                                                                       >  1:4975865/1‑150 (MQ=255)
                                                                 cACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAg                                                                                      <  1:2721808/150‑1 (MQ=255)
                                                                      aTGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACt                                                                                 >  2:1817133/1‑150 (MQ=255)
                                                                        ggCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTgg                                                                               >  2:3974011/1‑150 (MQ=255)
                                                                          cTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCa                                                                             >  2:4448360/1‑150 (MQ=255)
                                                                          cTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCa                                                                             >  1:2176334/1‑150 (MQ=255)
                                                                          cTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCa                                                                             >  2:3858432/1‑150 (MQ=255)
                                                                          cTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCa                                                                             <  1:3789016/150‑1 (MQ=255)
                                                                          cTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCa                                                                             <  2:5140634/150‑1 (MQ=255)
                                                                          cTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCa                                                                             >  2:5081026/1‑150 (MQ=255)
                                                                           tttGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAg                                                                            <  1:1481439/150‑1 (MQ=255)
                                                                           tttGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAg                                                                            >  1:4735112/1‑150 (MQ=255)
                                                                             tGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCt                                                                          <  2:3299917/150‑1 (MQ=255)
                                                                             tGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCt                                                                          <  2:120377/150‑1 (MQ=255)
                                                                              gCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTc                                                                         <  1:1291025/150‑1 (MQ=255)
                                                                                aTTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCaa                                                                       <  1:2526480/150‑1 (MQ=255)
                                                                                aTTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCaa                                                                       >  1:2840624/1‑150 (MQ=255)
                                                                                   gCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGAGGTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCCAAGc                                                                    >  2:1808130/1‑150 (MQ=255)
                                                                                     gtAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCACCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCAc                                                                  <  2:4882193/149‑1 (MQ=255)
                                                                                      tAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGATGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACt                                                                 >  1:3324463/1‑150 (MQ=255)
                                                                                      tAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGATGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACt                                                                 <  2:1036113/150‑1 (MQ=255)
                                                                                      tAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACt                                                                 >  2:1895800/1‑150 (MQ=255)
                                                                                       aaCCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTg                                                                <  1:4349963/150‑1 (MQ=255)
                                                                                         ccTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCa                                                              <  1:3714377/150‑1 (MQ=255)
                                                                                           tGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCaaa                                                            >  1:3092249/1‑150 (MQ=255)
                                                                                             gAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGATGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAAt                                                          <  2:2514767/150‑1 (MQ=255)
                                                                                               agagGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATgc                                                        <  2:2792102/150‑1 (MQ=255)
                                                                                                gagGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATgcg                                                       <  2:3080615/150‑1 (MQ=255)
                                                                                                  ggCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGcgcg                                                     >  2:1181867/1‑150 (MQ=255)
                                                                                                   gCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGATGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCgg                                                    <  1:4960428/150‑1 (MQ=255)
                                                                                                     aaTAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCACAATGCGCGGCt                                                  >  2:3523379/1‑150 (MQ=255)
                                                                                                     aaTAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCt                                                  <  1:1217563/150‑1 (MQ=255)
                                                                                                      aTAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTa                                                 <  1:2317536/150‑1 (MQ=255)
                                                                                                      aTAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTa                                                 <  2:1478981/150‑1 (MQ=255)
                                                                                                          cGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAg                                             >  1:4021387/1‑150 (MQ=255)
                                                                                                             tACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCAc                                          >  1:1099456/1‑150 (MQ=255)
                                                                                                              aCGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCAcc                                         <  1:5084095/150‑1 (MQ=255)
                                                                                                              aCGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCAcc                                         >  2:1345633/1‑150 (MQ=255)
                                                                                                               cGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCg                                        >  1:1013067/1‑150 (MQ=255)
                                                                                                                 tGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGtt                                      >  1:3523716/1‑150 (MQ=255)
                                                                                                                   aGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCAAAAAGCACCGTTgg                                    >  1:2278620/1‑150 (MQ=255)
                                                                                                                    gcgcGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGa                                   >  1:800765/1‑150 (MQ=255)
                                                                                                                     cgcgGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGat                                  >  1:2586028/1‑150 (MQ=255)
                                                                                                                         gCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATAttt                              <  2:1322408/150‑1 (MQ=255)
                                                                                                                           tGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTAt                            <  2:5186937/150‑1 (MQ=255)
                                                                                                                             gAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCt                          >  1:2430222/1‑150 (MQ=255)
                                                                                                                              aaaGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCtt                         <  2:4666327/150‑1 (MQ=255)
                                                                                                                                 ggACCGATCCTGATGCTGGATGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTaaa                      >  1:1060908/1‑150 (MQ=255)
                                                                                                                                   aCCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGt                    >  1:608163/1‑150 (MQ=255)
                                                                                                                                   aCCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGt                    <  2:2588878/150‑1 (MQ=255)
                                                                                                                                      gATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGaa                 <  2:4684018/150‑1 (MQ=255)
                                                                                                                                         ccTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAg              <  1:4183544/150‑1 (MQ=255)
                                                                                                                                            gATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTgg           <  1:1892659/150‑1 (MQ=255)
                                                                                                                                             aTGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTggg          >  2:2258045/1‑150 (MQ=255)
                                                                                                                                             aTGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTggg          <  1:4876850/150‑1 (MQ=255)
                                                                                                                                                 tGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGa      >  1:4078206/1‑150 (MQ=255)
                                                                                                                                                   gAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAAt    <  1:997824/150‑1 (MQ=255)
                                                                                                                                                   gAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAAt    <  1:3829121/150‑1 (MQ=255)
                                                                                                                                                    aaGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATc   >  2:495338/1‑150 (MQ=255)
                                                                                                                                                    aaGGATTTTTCCATGCTCAGGATCGTGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATc   <  2:3003725/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGTCGGTGGTAAAAGCGAACGCTTACGGGCATGGTATTGAGCGTATCTGGAGCGCGATCGGGGCCACCGATGGCTTTGCATTGCTTAACCTGGAAGAGGCAATAACGTTACGTGAGCGCGGCTGGAAAGGACCGATCCTGATGCTGGAAGGATTTTTCCATGCTCAGGATCTGGAGATTTATGACCAGCACCGCCTGACCACCTGCGTACACAGCAACTGGCAGCTCAAAGCACTGCAAAATGCGCGGCTAAAAGCACCGTTGGATATTTATCTTAAAGTGAACAGTGGGATGAATCG  >  CP009273/1234425‑1234723

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: