Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,987,928 C→A 15.3% 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,9280CA15.3% 824.6 / 91.3 261pseudogene (524/633 nt)ygeKpseudogene; putative 2‑component transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base C (121/100);  new base A (16/24);  total (137/124)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.27e-01

Only 100 of 263 total aligned reads displayed.

GTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTATCAAAATAAGTG  >  CP009273/2987781‑2988077
                                                                                                                                                   |                                                                                                                                                     
gTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCtt                                                                                                                                                     >  2:2128838/1‑150 (MQ=255)
gTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCtt                                                                                                                                                     <  2:3264175/150‑1 (MQ=255)
  aaGGATTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTaa                                                                                                                                                   <  1:1385107/150‑1 (MQ=255)
       ttCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCa                                                                                                                                              <  2:1585082/150‑1 (MQ=255)
       ttCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCa                                                                                                                                              <  2:3679358/150‑1 (MQ=255)
           gCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGaa                                                                                                                                          <  1:2090284/150‑1 (MQ=255)
            cAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGaaa                                                                                                                                         >  2:201797/1‑150 (MQ=255)
              aaaTGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGaaagt                                                                                                                                       >  1:2030359/1‑148 (MQ=255)
               aaTGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTa                                                                                                                                      >  2:2494864/1‑150 (MQ=255)
                 tGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGc                                                                                                                                    >  2:2950838/1‑150 (MQ=255)
                  gAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCg                                                                                                                                   >  1:2238766/1‑150 (MQ=255)
                    gCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGAt                                                                                                                                 <  2:501244/150‑1 (MQ=255)
                    gCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGAt                                                                                                                                 <  1:163885/150‑1 (MQ=255)
                       tAAAGCCTCTATTTTTTATATAGTGCAAACACCCATACGTAAAGCTGATTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTTCAAGAAATTAGCGATTNc                                                                                                                              <  2:74791/150‑1 (MQ=255)
                        aaaGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCt                                                                                                                             <  1:3428402/150‑1 (MQ=255)
                            ccactTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCtttgt                                                                                                                         >  1:4018124/4‑150 (MQ=255)
                             ctctTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCtttgtt                                                                                                                        >  1:407410/1‑150 (MQ=255)
                              tctTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTgtttgttt                                                                                                                       >  1:3298556/1‑150 (MQ=255)
                              tctTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCtttgttt                                                                                                                       >  2:450610/1‑150 (MQ=255)
                              tctTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCtttgttt                                                                                                                       <  2:1096583/150‑1 (MQ=255)
                                   tttttATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGtttt                                                                                                                  <  2:3608834/150‑1 (MQ=255)
                                   tttttATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGtttt                                                                                                                  <  2:4392799/150‑1 (MQ=255)
                                    ttttATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTa                                                                                                                 >  1:292810/1‑150 (MQ=255)
                                         tataGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATc                                                                                                            <  2:1741266/150‑1 (MQ=255)
                                          ataGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCa                                                                                                           <  1:567533/150‑1 (MQ=255)
                                           taGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAg                                                                                                          <  1:877646/150‑1 (MQ=255)
                                             gTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCt                                                                                                        <  2:614078/150‑1 (MQ=255)
                                              tGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTa                                                                                                       <  1:1017705/150‑1 (MQ=255)
                                                cAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATg                                                                                                     >  2:3552123/1‑150 (MQ=255)
                                                 aaaCACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAATAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGa                                                                                                    >  2:448292/1‑150 (MQ=255)
                                                 aaaCACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGa                                                                                                    <  1:1185199/150‑1 (MQ=255)
                                                  aaCACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAg                                                                                                   >  2:2769586/1‑150 (MQ=255)
                                                  aaCACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAAAATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACAATCAGCTATGAg                                                                                                   >  1:4020049/1‑150 (MQ=255)
                                                   acacCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGt                                                                                                  >  2:1805919/1‑150 (MQ=255)
                                                    cacCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGtt                                                                                                 >  2:3004745/1‑150 (MQ=255)
                                                     acCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTg                                                                                                <  1:3232768/150‑1 (MQ=255)
                                                      cccATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGt                                                                                               >  2:1119215/1‑150 (MQ=255)
                                                       ccATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTa                                                                                              <  2:1022595/150‑1 (MQ=255)
                                                        cATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTGTGTTTTACCATCAGCTATGAGTTGTaa                                                                                             >  2:339382/1‑150 (MQ=255)
                                                         aTACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAat                                                                                            >  2:1816430/1‑150 (MQ=255)
                                                         aTACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAat                                                                                            <  2:4624364/150‑1 (MQ=255)
                                                            cGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATAtt                                                                                         >  1:290947/1‑150 (MQ=255)
                                                            cGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATAtt                                                                                         <  2:285857/150‑1 (MQ=255)
                                                            cGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATAtt                                                                                         >  1:922326/1‑150 (MQ=255)
                                                             gTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATAttt                                                                                        >  1:4165901/1‑150 (MQ=255)
                                                             gTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATAttt                                                                                        >  1:3257895/1‑150 (MQ=255)
                                                               aaaGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCt                                                                                      >  1:3103175/1‑150 (MQ=255)
                                                                  gCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGAc                                                                                   >  1:464651/1‑150 (MQ=255)
                                                                  gCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGAc                                                                                   >  1:3521047/1‑150 (MQ=255)
                                                                      ttttGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCAc                                                                               <  2:2670623/150‑1 (MQ=255)
                                                                       tttGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACg                                                                              >  1:2714413/1‑150 (MQ=255)
                                                                                cTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAg                                                                     >  2:4294862/1‑150 (MQ=255)
                                                                                  cGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGtt                                                                   >  1:720957/1‑150 (MQ=255)
                                                                                   gATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTg                                                                  >  1:1012264/1‑150 (MQ=255)
                                                                                    aTACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGa                                                                 >  2:3339316/1‑150 (MQ=255)
                                                                                    aTACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGa                                                                 >  2:1369289/1‑150 (MQ=255)
                                                                                    aTACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGa                                                                 >  1:3795863/1‑150 (MQ=255)
                                                                                      aCCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGAtt                                                               >  2:2580151/1‑150 (MQ=255)
                                                                                      aCCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGAtt                                                               >  1:3568611/1‑150 (MQ=255)
                                                                                        cACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTAt                                                             >  2:3867907/1‑150 (MQ=255)
                                                                                        cACTGTGAACATATAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTAt                                                             <  2:3456407/150‑1 (MQ=255)
                                                                                          cTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCa                                                           <  1:2697478/150‑1 (MQ=255)
                                                                                           tttgAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCag                                                          <  2:2653277/148‑1 (MQ=255)
                                                                                            gtgAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCaga                                                         <  2:4022089/150‑1 (MQ=255)
                                                                                             tgAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCagag                                                        >  1:3593053/1‑150 (MQ=255)
                                                                                             tgAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCagag                                                        <  1:1422544/150‑1 (MQ=255)
                                                                                              gAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAgg                                                       >  1:3267685/1‑150 (MQ=255)
                                                                                              gAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAgg                                                       <  2:2847209/150‑1 (MQ=255)
                                                                                               aaCATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACAGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGa                                                      <  2:1059960/150‑1 (MQ=255)
                                                                                                aCATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGaa                                                     <  2:2527095/150‑1 (MQ=255)
                                                                                                aCATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGaa                                                     <  2:245044/150‑1 (MQ=255)
                                                                                                 cATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAAt                                                    <  2:114958/150‑1 (MQ=255)
                                                                                                 cATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAAt                                                    <  1:2128868/150‑1 (MQ=255)
                                                                                                     tAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCt                                                <  2:4338217/150‑1 (MQ=255)
                                                                                                     tAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCt                                                <  1:191423/150‑1 (MQ=255)
                                                                                                     tAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCt                                                <  1:1730083/150‑1 (MQ=255)
                                                                                                        tttCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGa                                             <  2:2946635/150‑1 (MQ=255)
                                                                                                         ttcttcATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGaa                                            >  2:842945/1‑150 (MQ=255)
                                                                                                            ttcATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGat                                         <  1:4237248/150‑1 (MQ=255)
                                                                                                              cATTATATTGAGTCTATGAGTTTCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGatat                                       <  1:3469482/150‑1 (MQ=255)
                                                                                                               aTGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATAtt                                      <  2:158138/150‑1 (MQ=255)
                                                                                                               aTGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATAtt                                      >  2:1306845/1‑150 (MQ=255)
                                                                                                               aTGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATAtt                                      >  2:1205330/1‑150 (MQ=255)
                                                                                                                tGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATAttt                                     <  1:2493961/150‑1 (MQ=255)
                                                                                                                tGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATAttt                                     >  1:146048/1‑150 (MQ=255)
                                                                                                                   tatTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATAtttttt                                  >  1:4202357/1‑150 (MQ=255)
                                                                                                                            tATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCa                         >  1:4321659/1‑150 (MQ=255)
                                                                                                                                aGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATg                     >  2:2152427/1‑150 (MQ=255)
                                                                                                                                   ttCGACTGTTTTCCTGCTTAACTGCAAGAAATTATCGATTTCTTTGTTTGTTTTACAATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTTTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTac                  <  1:2976022/150‑1 (MQ=255)
                                                                                                                                   ttCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTac                  <  2:1211035/150‑1 (MQ=255)
                                                                                                                                    tCGACTGTTTTCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTaca                 >  2:1402774/1‑150 (MQ=255)
                                                                                                                                      gACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACt               <  1:630392/150‑1 (MQ=255)
                                                                                                                                        cTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTAt             >  2:3939430/1‑150 (MQ=255)
                                                                                                                                        cTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTAt             <  1:253691/150‑1 (MQ=255)
                                                                                                                                        cTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTAt             >  2:2282780/1‑150 (MQ=255)
                                                                                                                                         tGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTATc            <  1:283166/150‑1 (MQ=255)
                                                                                                                                           ttttCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTATCaa          <  2:4018124/150‑1 (MQ=255)
                                                                                                                                              tCCTGATTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTATCAAAAt       <  2:675362/150‑1 (MQ=255)
                                                                                                                                                  gCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTATCAAAATAAgt   >  1:424706/1‑150 (MQ=255)
                                                                                                                                                   cTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTATCAAAATAAgtg  >  2:2356064/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GTAAGGCTTCAGCAAAATGAGCCTAAAGCCTCTTTTTTTTATATAGTGCAAACACCCATACGTAAAGCTGTTTTGATCAACTCGATACCACTGTGAACATCTAATTTCTTCATGATATTGAGTCTATGAGTTTCGACTGTTTTCCTGCTTAACTGCAAGAAATTAGCGATTTCTTTGTTTGTTTTACCATCAGCTATGAGTTGTAATATTTCTGACTCACGATTTGTAAGTTGATTATCAGAGGAATACCTGGAAGATATTTTTTCACAATCCATATGTACACTATCAAAATAAGTG  >  CP009273/2987781‑2988077

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: