Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,397,506 C→T 5.4% E3K (GAA→AAA)  abgB ← p‑aminobenzoyl‑glutamate hydrolase, B subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,397,5060CT 5.4% 1128.1 / 33.6 313E3K (GAA→AAA) abgBp‑aminobenzoyl‑glutamate hydrolase, B subunit
Reads supporting (aligned to +/- strand):  ref base C (165/130);  new base T (7/10);  total (173/140)
Fisher's exact test for biased strand distribution p-value = 3.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

Only 100 of 320 total aligned reads displayed.

GAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAACACCACGTAGGAGGCTTG  >  CP009273/1397357‑1397655
                                                                                                                                                     |                                                                                                                                                     
gAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTc                                                                                                                                                       >  1:3015630/1‑150 (MQ=255)
 aaGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTcc                                                                                                                                                      <  2:1380358/150‑1 (MQ=255)
 aaGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTcc                                                                                                                                                      >  2:4919350/1‑150 (MQ=255)
       gCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATc                                                                                                                                                <  1:2211191/150‑1 (MQ=255)
         aGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAg                                                                                                                                              >  1:2862656/1‑150 (MQ=255)
             tCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATAc                                                                                                                                          <  1:1071888/150‑1 (MQ=255)
                aGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTc                                                                                                                                       >  2:891405/1‑150 (MQ=255)
                  cgcCGAAGCCAGATGCTCCGCTGACCAGAACTCGTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTcgc                                                                                                                                     <  1:1235900/150‑1 (MQ=255)
                      gAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTcc                                                                                                                                 <  1:747633/150‑1 (MQ=255)
                       aaGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCa                                                                                                                                >  1:2388943/1‑150 (MQ=255)
                            aGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGaa                                                                                                                           >  2:2639433/1‑150 (MQ=255)
                             gATGCTCCGCTGACAAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGaaa                                                                                                                          <  2:602560/150‑1 (MQ=255)
                              aTGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGaaaa                                                                                                                         >  2:2766912/1‑150 (MQ=255)
                                gCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAAtt                                                                                                                       >  1:4254900/1‑150 (MQ=255)
                                  tCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGa                                                                                                                     <  1:3790769/150‑1 (MQ=255)
                                  tCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGa                                                                                                                     <  1:1174263/150‑1 (MQ=255)
                                       tGACCAGAACTCTTCAATACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGTTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGcgcg                                                                                                                <  2:4817278/150‑1 (MQ=255)
                                         aCCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGt                                                                                                              >  2:4186141/1‑150 (MQ=255)
                                           cAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGTGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTgc                                                                                                            >  2:3747389/1‑150 (MQ=255)
                                            aGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTgcg                                                                                                           >  2:361033/1‑150 (MQ=255)
                                             gAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTgcgc                                                                                                          >  2:4405048/1‑150 (MQ=255)
                                             gAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTgcgc                                                                                                          >  2:1797537/1‑150 (MQ=255)
                                               aCTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTgcgcgc                                                                                                        >  1:4685797/1‑150 (MQ=255)
                                                   ttCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGc                                                                                                    >  1:341714/1‑150 (MQ=255)
                                                      aaaaCGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGtt                                                                                                 <  1:2918391/150‑1 (MQ=255)
                                                        aaCGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTc                                                                                               <  2:1849557/150‑1 (MQ=255)
                                                            tgtTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACg                                                                                           >  1:792543/1‑150 (MQ=255)
                                                            tgtTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACg                                                                                           <  2:2904253/150‑1 (MQ=255)
                                                                 cTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAAt                                                                                      >  2:3449964/1‑150 (MQ=255)
                                                                   ggATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCACTAg                                                                                    >  2:4305173/1‑150 (MQ=255)
                                                                    gatgatCCCAGATTTGATCGGCAATATCGGTATAACTTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGc                                                                                   >  1:1452483/1‑150 (MQ=255)
                                                                     atgatCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCg                                                                                  >  1:4460933/1‑150 (MQ=255)
                                                                     atgatCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCg                                                                                  <  1:3872515/150‑1 (MQ=255)
                                                                          cccAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAAc                                                                             <  2:109190/150‑1 (MQ=255)
                                                                           ccAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAAcc                                                                            <  2:3366610/150‑1 (MQ=255)
                                                                            cAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCt                                                                           >  1:811096/1‑150 (MQ=255)
                                                                            cAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCt                                                                           <  2:2255428/150‑1 (MQ=255)
                                                                            cAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCt                                                                           <  1:4715161/150‑1 (MQ=255)
                                                                              gATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGc                                                                         >  2:5057010/1‑150 (MQ=255)
                                                                              gATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGc                                                                         >  2:4683044/1‑150 (MQ=255)
                                                                               aTTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCt                                                                        >  1:5128437/1‑150 (MQ=255)
                                                                                tttGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTc                                                                       >  2:2697108/1‑150 (MQ=255)
                                                                                tttGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTc                                                                       <  1:1300941/150‑1 (MQ=255)
                                                                                tttGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTc                                                                       <  1:450321/150‑1 (MQ=255)
                                                                                tttGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTc                                                                       <  1:3084854/150‑1 (MQ=255)
                                                                                  tGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGt                                                                     >  2:4371960/1‑150 (MQ=255)
                                                                                   gATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTc                                                                    <  2:2837480/150‑1 (MQ=255)
                                                                                    aTCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCa                                                                   <  1:2229885/150‑1 (MQ=255)
                                                                                       ggCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCaaaa                                                                >  1:4915003/1‑150 (MQ=255)
                                                                                        gCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTTCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAAt                                                               >  1:1920901/1‑150 (MQ=255)
                                                                                         cAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTTCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATc                                                              <  1:1915035/150‑1 (MQ=255)
                                                                                         cAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGACAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCACAATc                                                              >  2:4478417/1‑150 (MQ=255)
                                                                                         cAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATc                                                              <  2:271755/150‑1 (MQ=255)
                                                                                          aaTATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGGACCTGCTCGTCAAATTCg                                                             >  1:2697867/1‑150 (MQ=255)
                                                                                               cGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTTCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAAttt                                                        <  2:3089052/150‑1 (MQ=255)
                                                                                               cGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAAttt                                                        <  1:1083628/150‑1 (MQ=255)
                                                                                               cGGTATAACGTCGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGAGACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCATCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAAttt                                                        <  2:3575477/150‑1 (MQ=255)
                                                                                                ggTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAAtttt                                                       >  1:2120547/1‑150 (MQ=255)
                                                                                                 gTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAAttttt                                                      <  1:3970535/150‑1 (MQ=255)
                                                                                                  tataACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTc                                                     <  1:1990944/150‑1 (MQ=255)
                                                                                                   ataACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTGTTCg                                                    >  2:4722327/1‑150 (MQ=255)
                                                                                                     aaCGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGtt                                                  <  2:2466145/150‑1 (MQ=255)
                                                                                                       cGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTtgt                                                >  2:2594996/1‑150 (MQ=255)
                                                                                                        gTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTtgtg                                               >  1:1513687/1‑150 (MQ=255)
                                                                                                          tGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGtgat                                             >  2:668840/1‑150 (MQ=255)
                                                                                                          tGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGtgat                                             <  2:2022989/150‑1 (MQ=255)
                                                                                                             cGATCGGCTACACACGCATCGTCGATAAATCGATAGATTTCCTGCATCAGATACCACGCGACAAGGGAAAATTGAGCCCGCAGCGCGCCAGCGTTTCGACGCCAATAGAGAGAACCTGCTCGTCAAAATCGAATTTTTAGTTGtgattta                                          >  2:5014169/1‑147 (MQ=255)
                                                                                                              gATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGAcc                                         >  1:1343382/1‑150 (MQ=255)
                                                                                                               aTCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCt                                        >  2:2230803/1‑150 (MQ=255)
                                                                                                                 cGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTgc                                      <  1:767753/150‑1 (MQ=255)
                                                                                                                   gCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTgccg                                    <  2:243443/150‑1 (MQ=255)
                                                                                                                   gCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTgccg                                    >  2:1521432/1‑150 (MQ=255)
                                                                                                                   gCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTgccg                                    >  1:3741522/1‑150 (MQ=255)
                                                                                                                         aTCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCt                              <  2:1492666/150‑1 (MQ=255)
                                                                                                                          tCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCtg                             <  1:2755823/150‑1 (MQ=255)
                                                                                                                            gCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCtgtg                           <  1:2567025/150‑1 (MQ=255)
                                                                                                                             cATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCtgtgt                          <  2:2082024/150‑1 (MQ=255)
                                                                                                                              aTCGTCGATAAAACGATAGTTTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCtgtgtg                         <  2:394946/150‑1 (MQ=255)
                                                                                                                               tcgtcgATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGc                        >  1:3229321/1‑150 (MQ=255)
                                                                                                                                cgtcgATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGcc                       >  1:590891/1‑150 (MQ=255)
                                                                                                                                   cgATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCaaa                    >  2:3505630/1‑150 (MQ=255)
                                                                                                                                   cgATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCaaa                    <  1:920898/150‑1 (MQ=255)
                                                                                                                                     aTAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGTAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAaca                  >  2:3201916/1‑150 (MQ=255)
                                                                                                                                     aTAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAaca                  >  2:1158047/1‑150 (MQ=255)
                                                                                                                                      tAAAACGATAGATTTTCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAacac                 >  1:3492803/1‑150 (MQ=255)
                                                                                                                                      tAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAacac                 >  1:4699476/1‑150 (MQ=255)
                                                                                                                                      tAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAacac                 >  1:4699599/1‑150 (MQ=255)
                                                                                                                                       aaaaCGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAAcacc                >  1:2353819/1‑150 (MQ=255)
                                                                                                                                       aaaaCGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAAcacc                >  2:746731/1‑150 (MQ=255)
                                                                                                                                        aaaCGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAAcacca               >  1:5203809/1‑150 (MQ=255)
                                                                                                                                          aCGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAACACCACg             >  1:249182/1‑150 (MQ=255)
                                                                                                                                             aTAGATTTTCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAACACCACGTag          >  1:3400393/1‑150 (MQ=255)
                                                                                                                                             aTAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAACACCACGTag          >  1:830368/1‑150 (MQ=255)
                                                                                                                                               aGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAACACCACGTagga        >  1:3030038/1‑150 (MQ=255)
                                                                                                                                                gATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAACACCACGTaggag       >  2:408370/1‑150 (MQ=255)
                                                                                                                                                gATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAACACCACGTaggag       >  1:2545976/1‑150 (MQ=255)
                                                                                                                                                 aTTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAACACCACGTaggagg      <  1:1614825/150‑1 (MQ=255)
                                                                                                                                                   ttCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAACACCACGTAGGAGGCt    <  1:183508/150‑1 (MQ=255)
                                                                                                                                                    tCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAACACCACGTAGGAGGCtt   >  1:5055751/1‑150 (MQ=255)
                                                                                                                                                    tCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAACACCACGTAGGAGGCtt   <  1:1658650/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GAAGCCTGCAGATTCCAGCGCCGAAGCCAGATGCTCCGCTGACCAGAACTCTTCAAAACGTGTTTCTGGATGATCCCAGATTTGATCGGCAATATCGGTATAACGTTGGCGATCGGCTTCAATCGCATCGTCGATAAAACGATAGATTTCCTGCATCAGATACCTCGCGTCCAGGGAAAATTGAGCGCGGTGCGCGCCAGCGTTTCGACGGCAATAGCGAGAACCTGCTCGTCAAAATCGAATTTTTCGTTGTGATGACCTGCCGCCAGCTGTGTGCCAAACACCACGTAGGAGGCTTG  >  CP009273/1397357‑1397655

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