Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,589,933 G→A 93.2% 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,9330GA93.2% 1027.5 / 48.4 296T55I (ACC→ATC) ypfJputative neutral zinc metallopeptidase
Reads supporting (aligned to +/- strand):  ref base G (6/14);  new base A (140/136);  total (146/150)
Fisher's exact test for biased strand distribution p-value = 1.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.96e-01

Only 100 of 297 total aligned reads displayed.

GTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCG  >  CP009273/2589786‑2590081
                                                                                                                                                   |                                                                                                                                                    
gTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGAtt                                                                                                                                                    >  2:3762702/1‑150 (MQ=255)
gTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGAtt                                                                                                                                                    <  1:3652777/150‑1 (MQ=255)
 tCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTa                                                                                                                                                   <  2:2059737/150‑1 (MQ=255)
  cTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTaa                                                                                                                                                  <  2:4782741/150‑1 (MQ=255)
        cATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAac                                                                                                                                            <  2:3351207/150‑1 (MQ=255)
         aTCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAaca                                                                                                                                           >  1:1670322/1‑150 (MQ=255)
          tcttctCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAacac                                                                                                                                          <  1:1175777/150‑1 (MQ=255)
             tctcGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAt                                                                                                                                       >  2:675812/1‑150 (MQ=255)
               tcGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtag                                                                                                                                     >  2:3510533/1‑150 (MQ=255)
                cGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagt                                                                                                                                    >  2:879445/1‑150 (MQ=255)
                 gAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagta                                                                                                                                   >  1:787076/1‑150 (MQ=255)
                 gAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagta                                                                                                                                   <  2:3445441/150‑1 (MQ=255)
                  aaCTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGTTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCAtagtag                                                                                                                                  >  1:5135693/1‑150 (MQ=255)
                       ttgtCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGc                                                                                                                             >  2:4706235/1‑150 (MQ=255)
                            cccAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCa                                                                                                                        >  1:4478238/1‑150 (MQ=255)
                            cccAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCa                                                                                                                        >  1:2556625/1‑150 (MQ=255)
                             ccAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAg                                                                                                                       <  2:2919637/150‑1 (MQ=255)
                               aGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCa                                                                                                                     >  1:2823422/1‑150 (MQ=255)
                               aGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCa                                                                                                                     >  2:3836786/1‑150 (MQ=255)
                                ggTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCAc                                                                                                                    <  1:4429651/150‑1 (MQ=255)
                                ggTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCAc                                                                                                                    >  1:4982318/1‑150 (MQ=255)
                                ggTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCAc                                                                                                                    <  1:5107692/150‑1 (MQ=255)
                                ggTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCAc                                                                                                                    >  1:1945402/1‑150 (MQ=255)
                                 gtgtCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACg                                                                                                                   <  1:1244876/150‑1 (MQ=255)
                                 gtgtCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACg                                                                                                                   <  1:1802041/150‑1 (MQ=255)
                                  tgtCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGa                                                                                                                  >  2:1352888/1‑150 (MQ=255)
                                   gtCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGAc                                                                                                                 >  1:4758316/1‑150 (MQ=255)
                                    tCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACt                                                                                                                >  1:4319150/1‑150 (MQ=255)
                                      ttCCGTGGGTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTAt                                                                                                              <  2:2764556/150‑1 (MQ=255)
                                       tCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATc                                                                                                             >  2:2409527/1‑150 (MQ=255)
                                            ggTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTaa                                                                                                        <  2:1261617/150‑1 (MQ=255)
                                            ggTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTaa                                                                                                        <  1:1552201/150‑1 (MQ=255)
                                              ttGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTaaaa                                                                                                      <  2:2709084/150‑1 (MQ=255)
                                                gCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATc                                                                                                    >  2:4030735/1‑150 (MQ=255)
                                                 ccAGAATCACCGAGGTGAATTTTGCTGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCa                                                                                                   <  2:4743439/150‑1 (MQ=255)
                                                  aaGAATCACCGAGTTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACAGAGTTGATAGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAg                                                                                                  >  2:228806/2‑150 (MQ=255)
                                                     aaTCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCCGccc                                                                                               <  2:1416415/150‑1 (MQ=255)
                                                      aTCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCg                                                                                              >  2:1118707/1‑150 (MQ=255)
                                                      aTCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCg                                                                                              >  1:1481022/1‑150 (MQ=255)
                                                        cACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGcc                                                                                            >  1:2950422/1‑150 (MQ=255)
                                                         aCCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGcat                                                                                           <  2:1809127/150‑3 (MQ=255)
                                                         aCCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCt                                                                                           >  2:2529124/1‑150 (MQ=255)
                                                           cGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCttt                                                                                         >  1:2136644/1‑150 (MQ=255)
                                                                  aaTTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCtt                                                                                  <  1:4104560/150‑1 (MQ=255)
                                                                  aaTTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCtt                                                                                  <  2:2503490/150‑1 (MQ=255)
                                                                   aTTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTg                                                                                 <  1:4606456/150‑1 (MQ=255)
                                                                    ttttGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTgg                                                                                <  2:3069200/150‑1 (MQ=255)
                                                                       tGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAg                                                                             <  2:4352764/150‑1 (MQ=255)
                                                                       tGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAg                                                                             >  2:1517859/1‑150 (MQ=255)
                                                                           gCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGa                                                                         >  1:883650/1‑150 (MQ=255)
                                                                           gCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGa                                                                         >  2:1723614/1‑150 (MQ=255)
                                                                             ttCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGaaa                                                                       <  1:2085550/150‑1 (MQ=255)
                                                                             ttCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGaaa                                                                       >  2:2777973/1‑150 (MQ=255)
                                                                              tCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGaaaa                                                                      >  1:3877401/1‑150 (MQ=255)
                                                                              tCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGaaaa                                                                      <  2:536322/150‑1 (MQ=255)
                                                                               cGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAAc                                                                     <  1:1009108/150‑1 (MQ=255)
                                                                                gTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAAcc                                                                    >  2:5089529/1‑150 (MQ=255)
                                                                                gTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAAcc                                                                    <  2:772655/150‑1 (MQ=255)
                                                                                 tCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCg                                                                   <  2:4204070/150‑1 (MQ=255)
                                                                                 tCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCg                                                                   >  1:4976381/1‑150 (MQ=255)
                                                                                      aTTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGAcc                                                              >  1:4226251/1‑150 (MQ=255)
                                                                                      aTTTGGGCTAATTGACAGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGAcc                                                              >  1:808826/1‑150 (MQ=255)
                                                                                         tGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGcc                                                           >  1:492901/1‑150 (MQ=255)
                                                                                         tGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACGATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGcc                                                           <  2:460710/150‑1 (MQ=255)
                                                                                            gCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCAt                                                        <  1:3146229/150‑1 (MQ=255)
                                                                                             cTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATa                                                       <  1:3322808/150‑1 (MQ=255)
                                                                                             cTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATa                                                       <  1:2325073/150‑1 (MQ=255)
                                                                                             cTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATa                                                       <  2:3263341/150‑1 (MQ=255)
                                                                                                aTTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGAt                                                    >  1:3976374/1‑150 (MQ=255)
                                                                                                aTTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGAt                                                    <  2:4860226/150‑1 (MQ=255)
                                                                                                aTTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGAt                                                    >  1:2427992/1‑150 (MQ=255)
                                                                                                aTTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGAt                                                    >  1:4309549/1‑150 (MQ=255)
                                                                                                  tGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATgg                                                  >  2:186161/1‑150 (MQ=255)
                                                                                                   gACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGa                                                 <  1:4672638/150‑1 (MQ=255)
                                                                                                     ccGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGacc                                               >  2:6070/1‑150 (MQ=255)
                                                                                                     ccGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGacc                                               <  1:1292374/150‑1 (MQ=255)
                                                                                                      cgcgTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGacca                                              >  1:327179/1‑150 (MQ=255)
                                                                                                         gTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGAccacca                                           >  2:1213930/1‑150 (MQ=255)
                                                                                                           tGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCaga                                         >  1:1841145/1‑150 (MQ=255)
                                                                                                              ttgttgttgGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCt                                      <  2:1270816/150‑1 (MQ=255)
                                                                                                               tgttgttgGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTg                                     <  1:529869/150‑1 (MQ=255)
                                                                                                                  tgttgGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATCGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTg                                  <  2:3387060/150‑1 (MQ=255)
                                                                                                                   gttgGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTgc                                 >  1:4865307/1‑150 (MQ=255)
                                                                                                                    ttgGGAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCCTAGATGGACCACCAGAGCTGTTgcg                                >  1:1115060/1‑150 (MQ=255)
                                                                                                                      gggAAACCGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGcc                              <  2:4430943/150‑1 (MQ=255)
                                                                                                                             cGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTc                       >  1:2686029/1‑150 (MQ=255)
                                                                                                                             cGGCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTc                       >  1:3232505/1‑150 (MQ=255)
                                                                                                                              ggCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCa                      >  2:4539010/1‑150 (MQ=255)
                                                                                                                              ggCTGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCAGTCTTCa                      >  1:205279/1‑150 (MQ=255)
                                                                                                                                 tGCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACa                   <  2:159879/150‑1 (MQ=255)
                                                                                                                                  gCCCGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACAt                  <  1:3938334/150‑1 (MQ=255)
                                                                                                                                    ccGGTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTg                >  2:1926435/1‑150 (MQ=255)
                                                                                                                                       gTCATCAACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCa             >  1:3521996/1‑150 (MQ=255)
                                                                                                                                         catcaACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACt           >  2:4051749/1‑150 (MQ=255)
                                                                                                                                           tcaACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACttt         <  2:205061/150‑1 (MQ=255)
                                                                                                                                           tcaACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACttt         <  2:1444640/150‑1 (MQ=255)
                                                                                                                                            caACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTc        >  2:3919787/1‑150 (MQ=255)
                                                                                                                                            caACCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTc        <  1:4047233/150‑1 (MQ=255)
                                                                                                                                              aCCCGATTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCAc      <  1:250784/150‑1 (MQ=255)
                                                                                                                                              aCCCGATTAAATCAACACCAGAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCTACATTGTCACTTTCAc      <  2:3618247/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GTCTTACCCATCTTCTCGAACTGTTGTCCCCAGGTGTCTTCCGTGGTTGCCAGAATCACCGAGGTGAATTTTGCGGCTTCGTCTTCATTTGGGCTAATTGACCGCGTTGATTGTTGTTGGGAAACCGGCTGCCCGGTCATCAACCCGGTTAAATCAACACCATAGTAGCCTGCAACCAGCACGACTATCAGTAAAATCAGCCCGCCTTTACCGCTTGGCAGGCGAAAACCGGGACCGCCCATAGATGGACCACCAGAGCTGTTGCGCCTGTCTTCAACATTGTCACTTTCACGTCG  >  CP009273/2589786‑2590081

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: