Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,589,933 G→A 100% T55I (ACC→ATC)  ypfJ ← putative neutral zinc metallopeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,589,9330GA96.9% 912.5 / 11.1 261T55I (ACC→ATC) ypfJputative neutral zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (1/7);  new base A (133/119);  total (135/126)
Fisher's exact test for biased strand distribution p-value = 3.15e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 261 total aligned reads displayed.

GTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGT  >  CP009273/2589786‑2590079
                                                                                                                                                   |                                                                                                                                                  
gTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGAtt                                                                                                                                                  <  1:2343997/150‑1 (MQ=255)
gTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGAtt                                                                                                                                                  <  1:3787994/150‑1 (MQ=255)
gTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGAtt                                                                                                                                                  >  1:1670290/1‑150 (MQ=255)
 tCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTa                                                                                                                                                 >  2:2199591/1‑150 (MQ=255)
   ttACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTaaa                                                                                                                                               <  2:3228438/150‑1 (MQ=255)
    tACCCATCTTCTCGAACTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCTCGTTGATTGTTGTTTGGAAACCGGCTGCCCGGTCATCAACCCGATTAAAt                                                                                                                                              >  2:1966455/1‑150 (MQ=255)
    tACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAAt                                                                                                                                              >  2:4092967/1‑150 (MQ=255)
    tACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAAt                                                                                                                                              <  2:3518598/150‑1 (MQ=255)
    tACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAAt                                                                                                                                              <  2:1885518/150‑1 (MQ=255)
         aTCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAaca                                                                                                                                         >  2:3792131/1‑150 (MQ=255)
          tcttctCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAacac                                                                                                                                        >  2:2531366/1‑150 (MQ=255)
          tcttctCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAacac                                                                                                                                        <  2:4603692/150‑1 (MQ=255)
          tcttctCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAacac                                                                                                                                        >  1:2702542/1‑150 (MQ=255)
              ctcGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAta                                                                                                                                    <  2:1297918/150‑1 (MQ=255)
               tcGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtag                                                                                                                                   <  1:566287/150‑1 (MQ=255)
               tcGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtag                                                                                                                                   >  1:3735659/1‑150 (MQ=255)
                cGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagt                                                                                                                                  >  2:3537010/1‑150 (MQ=255)
                  aaCTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagtag                                                                                                                                >  2:1186761/1‑150 (MQ=255)
                  aaCTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagtag                                                                                                                                >  1:126745/1‑150 (MQ=255)
                   aCTTTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGc                                                                                                                               <  1:3820854/150‑1 (MQ=255)
                          tCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAAc                                                                                                                        <  1:4427015/150‑1 (MQ=255)
                            cccAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCa                                                                                                                      <  1:1316608/150‑1 (MQ=255)
                             ccAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAg                                                                                                                     <  1:3551731/150‑1 (MQ=255)
                              cAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGc                                                                                                                    <  1:4462988/150‑1 (MQ=255)
                                 gtgtCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACg                                                                                                                 <  1:1509013/150‑1 (MQ=255)
                                  tgtCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGa                                                                                                                >  1:1255298/1‑150 (MQ=255)
                                   gtCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGAc                                                                                                               >  2:3443762/1‑150 (MQ=255)
                                      ttCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGCTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTAt                                                                                                            >  2:63831/1‑150 (MQ=255)
                                      ttCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTAt                                                                                                            <  2:2069230/150‑1 (MQ=255)
                                        ccGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCa                                                                                                          >  1:1964322/1‑150 (MQ=255)
                                          gTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGt                                                                                                        >  1:1028887/1‑150 (MQ=255)
                                              ttGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTaaaa                                                                                                    <  1:1831017/150‑1 (MQ=255)
                                                 ccAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCa                                                                                                 <  2:4283142/150‑1 (MQ=255)
                                                   aGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGc                                                                                               >  1:4211683/1‑150 (MQ=255)
                                                     aaTCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGccc                                                                                             <  2:1063601/150‑1 (MQ=255)
                                                      aTCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCg                                                                                            <  2:2905273/150‑1 (MQ=255)
                                                        cACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGcc                                                                                          >  2:3870208/1‑150 (MQ=255)
                                                         aCCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCt                                                                                         >  1:3016204/1‑150 (MQ=255)
                                                          ccGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCtt                                                                                        >  1:959545/1‑150 (MQ=255)
                                                            gAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTa                                                                                      >  2:546552/1‑150 (MQ=255)
                                                             aGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTAc                                                                                     >  2:4632120/1‑150 (MQ=255)
                                                              ggTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTAcc                                                                                    >  2:2916019/1‑150 (MQ=255)
                                                              ggTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTAcc                                                                                    <  2:2461166/150‑1 (MQ=255)
                                                              ggTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTAcc                                                                                    <  1:2508043/150‑1 (MQ=255)
                                                                 gAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCt                                                                                 <  2:3894562/150‑1 (MQ=255)
                                                                  aaTTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCtt                                                                                <  1:33508/150‑1 (MQ=255)
                                                                   aTTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTg                                                                               >  1:599995/1‑150 (MQ=255)
                                                                    ttttGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTgg                                                                              <  1:1449387/150‑1 (MQ=255)
                                                                    ttgtGCGGCTTCGTCTGCATTTGGGCTAATTGACCGCGTTGATTGTTGTGGGGAACCCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTgg                                                                              <  1:2899963/147‑1 (MQ=255)
                                                                     tttGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGc                                                                             <  2:4421821/150‑1 (MQ=255)
                                                                         cGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGc                                                                         <  1:887932/150‑1 (MQ=255)
                                                                         cGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGc                                                                         >  2:3191492/1‑150 (MQ=255)
                                                                          ggCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCg                                                                        <  1:2294043/150‑1 (MQ=255)
                                                                            cTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGaa                                                                      >  1:1697749/1‑150 (MQ=255)
                                                                                gTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAAcc                                                                  >  1:1072585/1‑150 (MQ=255)
                                                                                gTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAAcc                                                                  <  2:4348316/150‑1 (MQ=255)
                                                                                gTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAAcc                                                                  >  2:80253/1‑150 (MQ=255)
                                                                                gTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAAcc                                                                  >  1:3058198/1‑150 (MQ=255)
                                                                                gTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAAcc                                                                  >  1:4016268/1‑150 (MQ=255)
                                                                                cgtctCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAAcc                                                                  <  2:2799478/146‑1 (MQ=255)
                                                                                   ttCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCggg                                                               <  2:2078561/150‑1 (MQ=255)
                                                                                    tCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGa                                                              <  1:2107435/150‑1 (MQ=255)
                                                                                      aTTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGAcc                                                            <  2:2755686/150‑1 (MQ=255)
                                                                                          gggCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGccc                                                        <  1:633780/150‑1 (MQ=255)
                                                                                             cTACTTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATa                                                     <  1:199703/150‑1 (MQ=255)
                                                                                             cTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATa                                                     <  1:1584378/150‑1 (MQ=255)
                                                                                             cTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGc                                                                         >  2:3691940/1‑130 (MQ=255)
                                                                                                 ttGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATg                                                 <  2:879426/150‑1 (MQ=255)
                                                                                                  tGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATgg                                                >  1:478835/1‑150 (MQ=255)
                                                                                                  tGACCGCGTTGAAAGTTGTTGGGAAACCGGCTGACCGGTCATAAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTAACAGTAAAATCAGCCCGCCTGTACCTCTTGTCAGGCGAAAACAGGGACCGCCCATAGATgg                                                >  2:2009334/1‑150 (MQ=255)
                                                                                                     ccGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGacc                                             <  2:3588058/150‑1 (MQ=255)
                                                                                                        cgTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGaccacc                                          >  2:2399543/1‑150 (MQ=255)
                                                                                                        cgTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGaccacc                                          >  1:3668668/1‑150 (MQ=255)
                                                                                                        cgTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGaccacc                                          >  1:852597/1‑150 (MQ=255)
                                                                                                         gTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGAccacca                                         <  1:4411037/150‑1 (MQ=255)
                                                                                                         gTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGAccacca                                         <  1:894079/150‑1 (MQ=255)
                                                                                                         gTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGAccacca                                         >  1:2417503/1‑150 (MQ=255)
                                                                                                          ttGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCag                                        >  2:68233/1‑150 (MQ=255)
                                                                                                           tGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCaga                                       >  1:3602153/1‑150 (MQ=255)
                                                                                                               tgttgttgGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTg                                   <  2:2288910/150‑1 (MQ=255)
                                                                                                                    ttgGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTgcg                              >  1:3454903/1‑150 (MQ=255)
                                                                                                                     tgGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTgcgc                             >  2:519291/1‑150 (MQ=255)
                                                                                                                      gggAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGcc                            >  1:1122490/1‑150 (MQ=255)
                                                                                                                      gggAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGcc                            >  2:3871082/1‑150 (MQ=255)
                                                                                                                       ggAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCt                           <  2:1980382/150‑1 (MQ=255)
                                                                                                                         aaaCCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGt                         >  1:777844/1‑150 (MQ=255)
                                                                                                                         aaaCCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGt                         >  2:501821/1‑150 (MQ=255)
                                                                                                                             cGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTc                     >  2:175671/1‑150 (MQ=255)
                                                                                                                              ggCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCa                    <  2:1538825/150‑1 (MQ=255)
                                                                                                                                 tGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACa                 >  1:559893/1‑150 (MQ=255)
                                                                                                                                   cccGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACAtt               >  2:3976986/1‑150 (MQ=255)
                                                                                                                                   cccGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACAtt               <  1:3703962/150‑1 (MQ=255)
                                                                                                                                   cccGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACAtt               >  2:2667946/1‑150 (MQ=255)
                                                                                                                                     cGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGt             >  1:321636/1‑150 (MQ=255)
                                                                                                                                     cGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGt             <  2:3027389/150‑1 (MQ=255)
                                                                                                                                     cGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCATCACCGCACGGCAGGCGACACCCGGGACCGCCCAGAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGt             <  2:1938758/150‑1 (MQ=255)
                                                                                                                                         catcaACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACt         <  2:4429440/150‑1 (MQ=255)
                                                                                                                                         catcaACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACt         >  2:1393281/1‑150 (MQ=255)
                                                                                                                                           tcaACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACttt       >  1:2721934/1‑150 (MQ=255)
                                                                                                                                             aaCCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCa     <  1:2610405/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGT  >  CP009273/2589786‑2590079

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: