Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,987,928 C→A 23.0% pseudogene (524/633 nt) ygeK ← pseudogene; putative 2‑component transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,987,9280CA23.0% 661.3 / 154.2 263pseudogene (524/633 nt)ygeKpseudogene; putative 2‑component transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (92/110);  new base A (29/32);  total (121/142)
Fisher's exact test for biased strand distribution p-value = 8.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 263 total aligned reads displayed.

CGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTATCAAAATAAGT  >  CP009273/2987780‑2988076
                                                                                                                                                    |                                                                                                                                                    
cGTAAGGCTTAAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGATGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCt                                                                                                                                                     >  1:3451551/1‑150 (MQ=255)
 gTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTACTTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCtt                                                                                                                                                    <  1:2678665/150‑1 (MQ=255)
 gTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCtt                                                                                                                                                    <  2:716175/150‑1 (MQ=255)
 gTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCtt                                                                                                                                                    <  2:4692008/150‑1 (MQ=255)
         tCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCaa                                                                                                                                            <  1:1645157/150‑1 (MQ=255)
           aGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGa                                                                                                                                          >  1:2421957/1‑150 (MQ=255)
            gCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGaa                                                                                                                                         <  1:916907/150‑1 (MQ=255)
            gCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGaa                                                                                                                                         <  2:1795577/150‑1 (MQ=255)
             cAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGaaa                                                                                                                                        <  2:3106719/150‑1 (MQ=255)
               aaaTGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAAtt                                                                                                                                      >  2:3058864/1‑150 (MQ=255)
               aaaTGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAAtt                                                                                                                                      >  1:1530905/1‑150 (MQ=255)
                   gAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCg                                                                                                                                  <  2:382473/150‑1 (MQ=255)
                   gAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCg                                                                                                                                  >  2:544291/1‑150 (MQ=255)
                    aGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGa                                                                                                                                 <  1:2695189/150‑1 (MQ=255)
                      ccTAAATCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGAtt                                                                                                                               <  2:4360309/150‑1 (MQ=255)
                       cTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGAttt                                                                                                                              <  2:2185276/150‑1 (MQ=255)
                       cTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGAttt                                                                                                                              <  1:1173734/150‑1 (MQ=255)
                       cTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCAATGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGAttt                                                                                                                              <  1:2363377/150‑1 (MQ=255)
                         aaaGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCt                                                                                                                            >  1:552875/1‑150 (MQ=255)
                              ctctTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCtttgtt                                                                                                                       <  2:284264/150‑1 (MQ=255)
                                     ttttATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTa                                                                                                                <  1:1669581/150‑1 (MQ=255)
                                     ttttATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTa                                                                                                                <  1:136984/150‑1 (MQ=255)
                                       ttATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTAcc                                                                                                              <  2:1130761/150‑1 (MQ=255)
                                       ttATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTAcc                                                                                                              <  2:5051479/150‑1 (MQ=255)
                                        tatataGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCa                                                                                                             >  2:1344775/1‑150 (MQ=255)
                                           ataGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCa                                                                                                          <  1:2161248/150‑1 (MQ=255)
                                           ataGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCa                                                                                                          <  2:2305561/150‑1 (MQ=255)
                                           ataGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCa                                                                                                          >  2:1084877/1‑150 (MQ=255)
                                           ataGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCa                                                                                                          >  2:94302/1‑150 (MQ=255)
                                           ataGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCa                                                                                                          <  2:3327339/150‑1 (MQ=255)
                                           ataGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCa                                                                                                          >  2:722694/1‑150 (MQ=255)
                                              gTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCt                                                                                                       <  2:841180/150‑1 (MQ=255)
                                              gTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCt                                                                                                       <  2:4320482/150‑1 (MQ=255)
                                              gTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCt                                                                                                       >  1:4316469/1‑150 (MQ=255)
                                              gTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCt                                                                                                       <  2:4320397/150‑1 (MQ=255)
                                                gCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTAt                                                                                                     <  1:417723/150‑1 (MQ=255)
                                                gCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTAt                                                                                                     <  1:4800897/150‑1 (MQ=255)
                                                 cAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATg                                                                                                    <  1:3306792/150‑1 (MQ=255)
                                                 cAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATg                                                                                                    >  1:3945655/1‑150 (MQ=255)
                                                  aaaCACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGa                                                                                                   >  1:3585709/1‑150 (MQ=255)
                                                  aaaCACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGa                                                                                                   >  2:2562165/1‑150 (MQ=255)
                                                    acacCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGt                                                                                                 >  1:441948/1‑150 (MQ=255)
                                                    acacCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGt                                                                                                 <  1:880180/150‑1 (MQ=255)
                                                    acacCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGt                                                                                                 >  2:257109/1‑150 (MQ=255)
                                                     cacCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGtt                                                                                                >  2:2925322/1‑150 (MQ=255)
                                                     cacCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGtt                                                                                                >  1:3529068/1‑150 (MQ=255)
                                                     cacCCATACGTAAAGCTGTTTTGATCAACTCGAGACCACTGTGAACATCTACTTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGtt                                                                                                <  1:2047277/150‑1 (MQ=255)
                                                         cATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTaa                                                                                            >  2:1954665/1‑150 (MQ=255)
                                                           tACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAata                                                                                          <  2:5042876/150‑1 (MQ=255)
                                                            aCGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAatat                                                                                         <  1:660661/150‑1 (MQ=255)
                                                                aaaGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCt                                                                                     >  1:5008486/1‑150 (MQ=255)
                                                                 aaGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTg                                                                                    <  1:2477504/150‑1 (MQ=255)
                                                                   gCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGAc                                                                                  >  2:1517878/1‑150 (MQ=255)
                                                                    cTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACt                                                                                 >  2:163237/1‑150 (MQ=255)
                                                                    cTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACt                                                                                 >  2:3126061/1‑150 (MQ=255)
                                                                        tttGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACg                                                                             <  2:2251650/150‑1 (MQ=255)
                                                                               aaCTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTa                                                                      >  2:44577/1‑150 (MQ=255)
                                                                                aCTCGATACCACTGTGAACATCTAATTTTTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTaa                                                                     <  1:2615969/150‑1 (MQ=255)
                                                                                aCTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTaa                                                                     >  1:3151075/1‑150 (MQ=255)
                                                                                 cTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAg                                                                    >  2:4049278/1‑150 (MQ=255)
                                                                                       aCCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGAtt                                                              <  1:50788/150‑1 (MQ=255)
                                                                                       aCCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGAtt                                                              <  1:1023779/150‑1 (MQ=255)
                                                                                            tgtgAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCag                                                         <  1:1099373/150‑1 (MQ=255)
                                                                                            tgtgAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCag                                                         <  1:26891/150‑1 (MQ=255)
                                                                                              tgAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCagag                                                       >  2:2979703/1‑150 (MQ=255)
                                                                                               gAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAgg                                                      <  1:1736127/150‑1 (MQ=255)
                                                                                                  cATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAAt                                                   <  1:4955470/150‑1 (MQ=255)
                                                                                                  cATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAAt                                                   >  2:4894656/1‑150 (MQ=255)
                                                                                                   aTCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATa                                                  >  2:2371250/1‑150 (MQ=255)
                                                                                                    tCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATAc                                                 >  2:5115993/1‑150 (MQ=255)
                                                                                                      tAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCt                                               <  1:2777602/150‑1 (MQ=255)
                                                                                                      tAATGTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCt                                               >  1:1751206/1‑150 (MQ=255)
                                                                                                          ttcttcATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGaa                                           >  2:379574/1‑150 (MQ=255)
                                                                                                          ttcttcATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGaa                                           <  1:207667/150‑1 (MQ=255)
                                                                                                              tcATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGata                                       <  1:243901/150‑1 (MQ=255)
                                                                                                                   atatTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATAttttt                                  <  1:2224369/150‑1 (MQ=255)
                                                                                                                     atTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTc                                <  2:3617995/150‑1 (MQ=255)
                                                                                                                       tGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTcac                              >  1:3880133/1‑150 (MQ=255)
                                                                                                                        gAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTcaca                             <  2:4206572/150‑1 (MQ=255)
                                                                                                                           tCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATc                          <  1:4673138/150‑1 (MQ=255)
                                                                                                                           tCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATc                          >  2:1180333/1‑150 (MQ=255)
                                                                                                                            cTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATcc                         >  1:298848/1‑150 (MQ=255)
                                                                                                                            cTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCATTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATcc                         <  2:2694678/150‑1 (MQ=255)
                                                                                                                            cTATGAGTTTCGACTGTTTTCCTGATTAACTGCATGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATcc                         <  1:815200/150‑1 (MQ=255)
                                                                                                                              aTGTGTTTAGACTGTTTTACTGATTTACTGCTATTTTTTAGCGTTTTGTTTGTTTGTTTTACCATCAGCTTTGAGTTGTAATATTTCTGACTGACGATTTGTAAGTTGTTTATCAGAGGAATACCTGGAAGATATTTTTTCACTATGCat                       <  2:1416947/150‑1 (MQ=255)
                                                                                                                                 aGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGCTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATg                    >  1:2312058/1‑150 (MQ=255)
                                                                                                                                 aGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATg                    <  2:424549/150‑1 (MQ=255)
                                                                                                                                    ttCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTac                 <  1:1955440/150‑1 (MQ=255)
                                                                                                                                     tCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTaca                <  1:1867088/150‑1 (MQ=255)
                                                                                                                                      cGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTacac               >  1:1322153/1‑150 (MQ=255)
                                                                                                                                      cGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTacac               >  1:1381126/1‑150 (MQ=255)
                                                                                                                                       gACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACt              >  2:81890/1‑150 (MQ=255)
                                                                                                                                         cTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTAt            <  1:3927605/150‑1 (MQ=255)
                                                                                                                                          tGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTATc           <  2:770937/150‑1 (MQ=255)
                                                                                                                                             tttCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTATCaaa        <  1:3615968/150‑1 (MQ=255)
                                                                                                                                              ttCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTATCaaaa       >  1:47629/1‑150 (MQ=255)
                                                                                                                                              ttCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTATCaaaa       >  2:2142240/1‑150 (MQ=255)
                                                                                                                                               tCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTATCAAAAt      <  1:3762266/150‑1 (MQ=255)
                                                                                                                                                ccTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTATCAAAATa     <  2:5092305/150‑1 (MQ=255)
                                                                                                                                                   gCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTATCAAAATAAgt  <  1:1212329/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTATCAAAATAAGT  >  CP009273/2987780‑2988076

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: