Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,429,000 C→A 8.5% intergenic (‑519/+215) ynaE ← / ← ttcC cold shock protein, Rac prophage/pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,429,0000CA 8.5% 752.6 / 29.7 203intergenic (‑519/+215)ynaE/ttcCcold shock protein, Rac prophage/pseudogene, prophage Rac integration site ttcA duplication;Phage or Prophage Related
Reads supporting (aligned to +/- strand):  ref base C (74/111);  new base A (5/13);  total (79/124)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 212 total aligned reads displayed.

CGGATTGGGTTTAAAAGTATACATAGATAACAGTTTTTATCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATG  >  CP009273/1428854‑1429149
                                                                                                                                                  |                                                                                                                                                     
cGGATTGGGTTTAAAAGTATACATAGATAACAGTTTTTATCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTa                                                                                                                                                    <  1:4524315/150‑1 (MQ=39)
   tgatattggttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGc                                                                                                                                                 >  2:1849459‑M1/38‑150 (MQ=255)
     tGGGTTTAAAAGTATACATAGATAACAGTTTTTATCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCaa                                                                                                                                               >  2:2684125/1‑150 (MQ=255)
      tattggttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAGTGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAg                                                                                                                                              <  2:3351685‑M1/116‑1 (MQ=255)
       attggttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGa                                                                                                                                             <  1:4373437‑M1/117‑1 (MQ=255)
       attggttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGa                                                                                                                                             <  1:4482824‑M1/117‑1 (MQ=255)
         tggttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAAc                                                                                                                                           <  1:236422‑M1/119‑1 (MQ=255)
         tggttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAAc                                                                                                                                           <  1:2626513‑M1/119‑1 (MQ=255)
          ggttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACg                                                                                                                                          >  1:1644062‑M1/31‑150 (MQ=255)
           gttgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGt                                                                                                                                         <  1:1866076‑M1/121‑1 (MQ=255)
             tgtctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGtgc                                                                                                                                       <  2:3725716‑M1/123‑1 (MQ=255)
                ctggagattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGtgctgc                                                                                                                                    <  1:1986822‑M1/126‑1 (MQ=255)
                    agattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGCTGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCggtt                                                                                                                                >  2:3059861‑M1/21‑150 (MQ=255)
                    agattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCggtt                                                                                                                                <  1:1288557‑M1/130‑1 (MQ=255)
                     gattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCggttg                                                                                                                               <  1:3156017‑M1/131‑1 (MQ=255)
                     gattcagggggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCggttg                                                                                                                               >  2:3454874‑M1/20‑150 (MQ=255)
                     cATAGATAACAGTTTTTATCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCggttg                                                                                                                               <  2:681203/150‑1 (MQ=255)
                      aTAGATAACAGTTTTTATCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCggttgg                                                                                                                              >  1:3860372/1‑150 (MQ=255)
                             gggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACt                                                                                                                       >  1:3227809‑M1/12‑150 (MQ=255)
                             gggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACt                                                                                                                       >  2:520563‑M1/12‑150 (MQ=255)
                              ggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACtt                                                                                                                      >  1:4156785‑M1/11‑150 (MQ=255)
                              ggccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACtt                                                                                                                      <  2:64316‑M1/140‑1 (MQ=255)
                               gccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACttt                                                                                                                     <  2:3575135‑M1/141‑1 (MQ=255)
                                ccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTttg                                                                                                                    <  2:2407205‑M1/142‑1 (MQ=255)
                                ccagtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTttg                                                                                                                    <  2:3138695‑M1/142‑1 (MQ=255)
                                   tttATCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTttgatt                                                                                                                 <  2:2267943/150‑1 (MQ=255)
                                   gtctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTttgatt                                                                                                                 >  2:3742845‑M1/6‑150 (MQ=255)
                                     ctaCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTttgattga                                                                                                               <  2:4316316‑M1/147‑1 (MQ=255)
                                      taCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTgag                                                                                                              >  2:4332170‑M1/3‑150 (MQ=255)
                                      taCTGAAGAAGAAAAATATCAAGGTGAGATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTgag                                                                                                              <  1:1548555‑M1/148‑1 (MQ=255)
                                       aCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTgaga                                                                                                             >  2:2944928‑M1/2‑150 (MQ=255)
                                            agaagaaAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTTTCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGttt                                                                                                        >  1:2633712/1‑150 (MQ=255)
                                            agaagaaAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGttt                                                                                                        <  1:547265/150‑1 (MQ=255)
                                              aagaaAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTg                                                                                                      <  1:718381/150‑1 (MQ=255)
                                               agaaAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTgg                                                                                                     <  1:4192026/150‑1 (MQ=255)
                                               agaaAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTgg                                                                                                     >  1:1423110/1‑150 (MQ=255)
                                                 aaaaaTATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGaa                                                                                                   <  2:143565/150‑1 (MQ=255)
                                                    aaTATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAAttt                                                                                                <  2:3455531/150‑1 (MQ=255)
                                                    aaTATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAAttt                                                                                                <  1:1761988/150‑1 (MQ=255)
                                                    aaTATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAAttt                                                                                                <  2:577413/150‑1 (MQ=255)
                                                     atatCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAAtttt                                                                                               <  2:2799434/150‑1 (MQ=255)
                                                     atatCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAAtttt                                                                                               <  2:1147009/150‑1 (MQ=255)
                                                      tatCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAAtttta                                                                                              <  1:4034692/150‑2 (MQ=255)
                                                      tatCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGATGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAAttttt                                                                                              >  2:1260159/1‑150 (MQ=255)
                                                       atCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAAtttttt                                                                                             <  1:3791241/150‑1 (MQ=255)
                                                       atCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAAtttttt                                                                                             <  1:3764162/150‑1 (MQ=255)
                                                            ggTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGt                                                                                        <  2:2163485/150‑1 (MQ=255)
                                                              tGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTgg                                                                                      <  1:1146608/150‑1 (MQ=255)
                                                              tGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTgg                                                                                      <  1:1146716/150‑1 (MQ=255)
                                                                 tataGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCaa                                                                                   >  2:1095320/1‑150 (MQ=255)
                                                                  ataGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCaaa                                                                                  >  1:2420585/1‑150 (MQ=255)
                                                                    aGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCaaaaa                                                                                >  2:1898927/1‑150 (MQ=255)
                                                                          tatGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGc                                                                          <  2:3778965/150‑1 (MQ=255)
                                                                          tatGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGc                                                                          >  1:4628942/1‑150 (MQ=255)
                                                                           atGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCa                                                                         >  2:2907505/1‑150 (MQ=255)
                                                                               cTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAAc                                                                     >  1:1020525/1‑150 (MQ=255)
                                                                               cTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAAc                                                                     <  1:263217/150‑1 (MQ=255)
                                                                                tttGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACg                                                                    <  2:352813/150‑1 (MQ=255)
                                                                                tttGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACg                                                                    <  2:205777/150‑1 (MQ=255)
                                                                                  tGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCt                                                                  >  1:1221375/1‑150 (MQ=255)
                                                                                   gATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTg                                                                 <  1:880760/150‑1 (MQ=255)
                                                                                    aTGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATTGGGCAAAACGCTGc                                                                >  1:967404/1‑150 (MQ=255)
                                                                                    aTGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGc                                                                >  1:966889/1‑150 (MQ=255)
                                                                                    aTGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGc                                                                <  1:3549603/150‑1 (MQ=255)
                                                                                         gaggaATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaag                                                           <  2:504162/150‑1 (MQ=255)
                                                                                             aatgaatgTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaaggggc                                                       >  1:4228022/1‑150 (MQ=255)
                                                                                                 aatgTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaaggggcaaaa                                                   <  1:1839160/150‑1 (MQ=255)
                                                                                                   tgtgATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaaggggcaaaaaa                                                 >  2:4061746/1‑150 (MQ=255)
                                                                                                      gATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGATACGTTTTGGAATTTTTTTCGGTGGCAAAAATTGTGCAAAAAGCTGAaaaaggggcaaaaaaggg                                              >  1:4399298/1‑150 (MQ=38)
                                                                                                        tGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaaggggcaaaaaaggggc                                            <  1:1056894/150‑1 (MQ=255)
                                                                                                         gggAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaaggggcaaaaaaggggca                                           <  2:1334881/150‑1 (MQ=255)
                                                                                                          ggAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaaggggcaaaaaaggggcaa                                          <  1:1083386/150‑1 (MQ=255)
                                                                                                          ggAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCaaaaggggcaaaaaaggggcaa                                          <  2:1595352/150‑1 (MQ=255)
                                                                                                           gAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAaaaggggcaaaaaaggggcaaa                                         <  1:2102846/150‑1 (MQ=255)
                                                                                                              tGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAggggcaaaaaaggggcaaaaaa                                      >  2:2032133/1‑150 (MQ=255)
                                                                                                                   tATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGt                                 <  2:4239846/150‑1 (MQ=255)
                                                                                                                   tATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGt                                 >  2:1427289/1‑150 (MQ=255)
                                                                                                                   tATCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTTCTGCGGTTGGATTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAGGGGCAAAAAAAGAGTg                                <  1:1815256/150‑1 (MQ=255)
                                                                                                                    aTCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTg                                <  2:2626513/150‑1 (MQ=255)
                                                                                                                    aTCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTg                                <  2:236422/150‑1 (MQ=255)
                                                                                                                    aTCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTg                                >  2:2441739/1‑150 (MQ=255)
                                                                                                                    aTCTGAATAGTTGAAAAACCGCAGTCACGTAGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTg                                <  1:2749285/150‑1 (MQ=255)
                                                                                                                        gAATAGTTGAACAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGAtt                            <  2:4720288/150‑1 (MQ=255)
                                                                                                                           tAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATc                         <  1:3454263/150‑1 (MQ=255)
                                                                                                                            aGTTGAAAAACCGCCGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCg                        <  1:2557135/150‑1 (MQ=255)
                                                                                                                              ttGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTa                      <  1:4402100/150‑1 (MQ=255)
                                                                                                                                 aaaaaCCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCt                   >  2:4563276/1‑150 (MQ=255)
                                                                                                                                  aaaaCCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCtt                  >  1:2881591/1‑150 (MQ=255)
                                                                                                                                   aaaCCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTa                 >  1:1689907/1‑150 (MQ=255)
                                                                                                                                   aaaCCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTa                 >  1:69314/1‑150 (MQ=255)
                                                                                                                                      ccGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTAttg              <  2:3983867/150‑1 (MQ=255)
                                                                                                                                      ccGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTAttg              <  2:3879345/150‑1 (MQ=255)
                                                                                                                                        gCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTAttgtt            <  2:1267060/150‑1 (MQ=255)
                                                                                                                                             cACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGGGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGc       <  1:2778507/150‑1 (MQ=255)
                                                                                                                                                gtcgtATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCtga    <  2:3307324/150‑1 (MQ=255)
                                                                                                                                                gtcgtATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCtga    <  2:913634/150‑1 (MQ=255)
                                                                                                                                                 tcgtATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAACATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCtgat   <  2:2632823/150‑1 (MQ=255)
                                                                                                                                                 tcgtATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCtgat   <  2:1548555/150‑1 (MQ=255)
                                                                                                                                                  cgtATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCtgatg  <  2:4482824/150‑1 (MQ=255)
                                                                                                                                                  cgtATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCtgatg  <  2:4373437/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
CGGATTGGGTTTAAAAGTATACATAGATAACAGTTTTTATCTGAAGAAGAAAAATATCAAGGTGATATAGCCTATATGCCTTTGATGCGGAGGAATGAATGTGATGGGAGTGATGTATCTGAATAGTTGAAAAACCGCAGTCACGTCGTATGCAAGAACGTGCTGCGGTTGGTTTGACTTTGATTGAGACGTTTTGGAATTTTTTTCGGTGGCAAAAATGGGGCAAAACGCTGCAAAAGGGGCAAAAAAGGGGCAAAAAAAGAGTGGATTATCGTAGCTTATTGTTGTCGCTGATG  >  CP009273/1428854‑1429149

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: