Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,076,526 C→A 78.8% D211Y (GAC→TAC)  speB ← agmatinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,076,5260CA78.8% 617.2 / 213.7 293D211Y (GAC→TAC) speBagmatinase
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base A (127/104);  total (157/136)
Fisher's exact test for biased strand distribution p-value = 3.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

Only 100 of 294 total aligned reads displayed.

TACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCA  >  CP009273/3076380‑3076674
                                                                                                                                                  |                                                                                                                                                    
tACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCAt                                                                                                                                                    <  1:1037724/149‑1 (MQ=255)
 aCCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATc                                                                                                                                                   <  1:3314555/149‑1 (MQ=255)
   cAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCa                                                                                                                                                 <  1:4796767/149‑1 (MQ=255)
    aGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCAc                                                                                                                                                <  1:1000345/149‑1 (MQ=255)
      tttAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGc                                                                                                                                              >  1:818610/1‑149 (MQ=255)
      tttAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGc                                                                                                                                              <  2:1243938/149‑1 (MQ=255)
       ttAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCt                                                                                                                                             <  2:4230893/149‑1 (MQ=255)
           tAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGa                                                                                                                                         >  1:1888695/1‑149 (MQ=255)
           tAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGa                                                                                                                                         <  1:1267922/149‑1 (MQ=255)
               gcgATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGt                                                                                                                                     <  1:1349216/149‑1 (MQ=255)
                  aTCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCa                                                                                                                                  >  2:446545/1‑149 (MQ=255)
                    cGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCAcc                                                                                                                                <  2:2469646/149‑1 (MQ=255)
                    cGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCAcc                                                                                                                                <  2:289704/149‑1 (MQ=255)
                    cGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCAcc                                                                                                                                <  1:2620649/149‑1 (MQ=255)
                        ggTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGc                                                                                                                            >  2:4889080/1‑149 (MQ=255)
                         gTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCa                                                                                                                           <  1:306907/149‑1 (MQ=255)
                                cgccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTcc                                                                                                                    <  2:2987895/149‑1 (MQ=255)
                                cgccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTcc                                                                                                                    <  1:1825482/149‑1 (MQ=255)
                                cgccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTcc                                                                                                                    <  1:4847641/149‑1 (MQ=255)
                                 gccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGAGATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCa                                                                                                                   >  1:919062/1‑149 (MQ=255)
                                  ccAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAg                                                                                                                  >  1:3168946/1‑149 (MQ=255)
                                    aaTCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCa                                                                                                                <  2:3715690/149‑1 (MQ=255)
                                        aCTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGt                                                                                                            >  2:2105283/1‑149 (MQ=255)
                                         cTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTa                                                                                                           >  1:5165789/1‑149 (MQ=255)
                                                aCCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGt                                                                                                    <  2:5135786/149‑1 (MQ=255)
                                                     tGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCt                                                                                               >  2:3760571/1‑149 (MQ=255)
                                                     tGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCt                                                                                               >  2:533562/1‑149 (MQ=255)
                                                        cTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCACAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTa                                                                                            <  2:4296417/149‑1 (MQ=255)
                                                         tGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTAt                                                                                           >  2:2934313/1‑149 (MQ=255)
                                                           gTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCa                                                                                         <  2:4589587/149‑1 (MQ=255)
                                                             gCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCaaa                                                                                       >  2:3967284/1‑149 (MQ=255)
                                                             gCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCaaa                                                                                       >  1:53475/1‑149 (MQ=255)
                                                              cAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAAc                                                                                      <  2:4695473/149‑1 (MQ=255)
                                                               aaaaGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACt                                                                                     <  2:1055232/149‑1 (MQ=255)
                                                                aaaGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTc                                                                                    >  2:2159346/1‑149 (MQ=255)
                                                                       gATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGa                                                                             >  1:1054946/1‑149 (MQ=255)
                                                                         tCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAAt                                                                           >  2:3492252/1‑149 (MQ=255)
                                                                          ccAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATa                                                                          <  1:3736335/149‑1 (MQ=255)
                                                                             ggCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCa                                                                       <  1:1655/149‑1 (MQ=255)
                                                                               cAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAAt                                                                     <  1:5040590/149‑1 (MQ=255)
                                                                               cAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAAt                                                                     >  2:1525847/1‑149 (MQ=255)
                                                                               cAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAAt                                                                     <  1:221010/149‑1 (MQ=255)
                                                                               cAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAAt                                                                     <  1:702649/149‑1 (MQ=255)
                                                                                aGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATc                                                                    <  2:1798067/149‑1 (MQ=255)
                                                                                 gTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCt                                                                   >  2:166938/1‑149 (MQ=255)
                                                                                 gTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGACTACCAATCt                                                                   >  1:4072742/1‑149 (MQ=255)
                                                                                 gTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCt                                                                   <  1:3513954/149‑1 (MQ=255)
                                                                                    gATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCa                                                                >  1:2185166/1‑149 (MQ=255)
                                                                                         cAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGaa                                                           <  2:2807892/149‑1 (MQ=255)
                                                                                         cAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGaa                                                           <  1:582477/149‑1 (MQ=255)
                                                                                          aaaaGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAAt                                                          <  2:4847339/149‑1 (MQ=255)
                                                                                           aaaGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATg                                                         <  1:1705862/149‑1 (MQ=255)
                                                                                             aGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGAt                                                       >  2:5064108/1‑149 (MQ=255)
                                                                                             aGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGAt                                                       >  1:4755944/1‑149 (MQ=255)
                                                                                              gTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGAtt                                                      <  2:903947/149‑1 (MQ=255)
                                                                                               tCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTc                                                     <  1:3535967/149‑1 (MQ=255)
                                                                                                   gTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGt                                                 <  2:369444/149‑1 (MQ=255)
                                                                                                      aaCCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGa                                              >  1:2672296/1‑149 (MQ=255)
                                                                                                      aaCCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGa                                              >  1:5022024/1‑149 (MQ=255)
                                                                                                      aaCCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGa                                              >  2:2307383/1‑149 (MQ=255)
                                                                                                       aCCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGAt                                             >  1:4659027/1‑149 (MQ=255)
                                                                                                        ccGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATc                                            <  2:1888695/149‑1 (MQ=255)
                                                                                                          ggCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAg                                          >  1:1937262/1‑149 (MQ=255)
                                                                                                              tatCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACct                                      <  2:306705/149‑1 (MQ=255)
                                                                                                               atCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACctt                                     >  1:3543075/1‑149 (MQ=255)
                                                                                                               atCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACctt                                     >  1:3785024/1‑149 (MQ=255)
                                                                                                                tCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACcttc                                    >  2:2426860/1‑149 (MQ=255)
                                                                                                                tCACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACcttc                                    <  1:2943282/149‑1 (MQ=255)
                                                                                                                 cACCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACcatct                                   >  1:4391682/1‑149 (MQ=255)
                                                                                                                  aCCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACcttctt                                  >  1:2462512/1‑149 (MQ=255)
                                                                                                                  aCCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACcttctt                                  <  2:4477132/149‑1 (MQ=255)
                                                                                                                  aCCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACcttctt                                  >  1:3871027/1‑149 (MQ=255)
                                                                                                                  aCCCACAATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACcttctt                                  <  2:1346034/149‑1 (MQ=255)
                                                                                                                    ccACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTc                                >  2:2745091/1‑149 (MQ=255)
                                                                                                                     cacaATCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCg                               >  2:1174140/1‑149 (MQ=255)
                                                                                                                      acaATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCgg                              >  2:2662773/1‑149 (MQ=255)
                                                                                                                         aTCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGgcgc                           <  1:2126751/149‑1 (MQ=255)
                                                                                                                         aTCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGgcgc                           >  2:2634689/1‑149 (MQ=255)
                                                                                                                         aTCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGgcgc                           <  2:2005669/149‑1 (MQ=255)
                                                                                                                          tCTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGcgcg                          >  1:538462/1‑149 (MQ=255)
                                                                                                                           cTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCgg                         <  2:3786712/149‑1 (MQ=255)
                                                                                                                           cTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCgg                         <  2:3786807/149‑1 (MQ=255)
                                                                                                                           cTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCgg                         >  1:199454/1‑149 (MQ=255)
                                                                                                                           cTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCgg                         >  1:4805977/1‑149 (MQ=255)
                                                                                                                           cTGTTTCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCgg                         >  1:3410855/1‑149 (MQ=255)
                                                                                                                              tttCACTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGtat                      >  1:64828/1‑149 (MQ=255)
                                                                                                                                  aCTTGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGaa                  <  1:1169450/149‑1 (MQ=255)
                                                                                                                                    ttGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACa                >  2:2188057/1‑149 (MQ=255)
                                                                                                                                    ttGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACa                >  2:196381/1‑149 (MQ=255)
                                                                                                                                     tGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACAt               <  1:955980/149‑1 (MQ=255)
                                                                                                                                     tGGGCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACAt               <  2:3199027/149‑1 (MQ=255)
                                                                                                                                      gggCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGACTGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATa              >  2:585306/1‑149 (MQ=255)
                                                                                                                                       ggCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAg             >  1:2873296/1‑149 (MQ=255)
                                                                                                                                       ggCGATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAg             >  1:3945693/1‑149 (MQ=255)
                                                                                                                                          gATAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGc          <  1:504968/149‑1 (MQ=255)
                                                                                                                                           aTAACGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGcc         <  1:309066/149‑1 (MQ=255)
                                                                                                                                             aaCGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGt       <  1:3685617/149‑1 (MQ=255)
                                                                                                                                              aCGTAATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTg      >  2:4077311/1‑149 (MQ=255)
                                                                                                                                                 tCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTc   >  2:2839122/1‑149 (MQ=255)
                                                                                                                                                  cATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCa  >  1:4228317/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TACCAGTTTAATAGCGCGATCGGAGGTCAGGCCGCCAATCACTGGCGTACCGGTGCCTGGTGCAAAAGCAGGATCCAGGCAGTCGATATCAAAAGTCAGGTAAACCGGCATATCACCCACAATCTGTTTCACTTGGGCGATAACGTCATCCACGCTGCGATCGTTCACCTGGCAGGCGTCCAGCACGGTAAAGCCGTTGTCTTTATCAAACTCGGTACGAATACCAATCTGCACGGAATGATTCGGGTCGATCAGACCTTCTTTCGGCGCGGTATAGAACATAGTGCCGTGGTCA  >  CP009273/3076380‑3076674

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: ‑