Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,499,887 G→C 5.5% pseudogene (905/1199 nt) yjgX ← KpLE2 phage‑like element; putative phosphoethanolamine transferase YjgX, N‑terminal fragment

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,499,8870GC 5.5% 619.5 / 18.9 181pseudogene (905/1199 nt)yjgXKpLE2 phage‑like element; putative phosphoethanolamine transferase YjgX, N‑terminal fragment
Reads supporting (aligned to +/- strand):  ref base G (88/83);  new base C (4/6);  total (92/89)
Fisher's exact test for biased strand distribution p-value = 5.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.03e-01

Only 100 of 181 total aligned reads displayed.

TCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTA  >  NC_000913/4499739‑4500033
                                                                                                                                                    |                                                                                                                                                  
  cACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGAcc                                                                                                                                                 >  1:847543/1‑150 (MQ=255)
  cACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGAcc                                                                                                                                                 >  2:2476160/1‑150 (MQ=255)
    cTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGa                                                                                                                                               >  2:641995/1‑150 (MQ=255)
     tGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGAc                                                                                                                                              >  2:27352/1‑150 (MQ=255)
      gAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACa                                                                                                                                             >  2:304572/1‑150 (MQ=255)
      gAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACa                                                                                                                                             >  1:379407/1‑150 (MQ=255)
        aaaaTACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGc                                                                                                                                           <  1:472436/150‑1 (MQ=255)
        aaaaTACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGc                                                                                                                                           >  1:2400414/1‑150 (MQ=255)
         aaaTACCTTCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGGTCAGCGACAgtg                                                                                                                                          <  2:1911700/150‑3 (MQ=255)
         aaaTACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCATGGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGTCACTGTCGCTGACCGACAGcc                                                                                                                                          <  2:1234383/150‑1 (MQ=255)
         aaaTACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGGTCAGCGACAgtg                                                                                                                                          <  2:1036256/150‑3 (MQ=255)
         aaaTACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGcc                                                                                                                                          <  2:550623/150‑1 (MQ=255)
         aaaTACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGcc                                                                                                                                          <  1:1643054/150‑1 (MQ=255)
          aaTACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGccc                                                                                                                                         <  2:2471684/150‑1 (MQ=255)
             aCCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGGTCAGCGACAgtgccaa                                                                                                                                      <  1:2472020/150‑7 (MQ=255)
             aCCTGCGAAATGACGGGCACTGATGTTTCTACTGACGGGCTGATGAGACGCTATGGCACTGGCGGTAATAATCAGTGGCACACTGTTGCCTCCCTGCACATGACCATAATGATGCACAGGATTGGCACTGTCGGTCAGCGACAgtgccaa                                                                                                                                      <  2:1187491/150‑7 (MQ=255)
                tGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGcc                                                                                                                                   >  1:986283/1‑150 (MQ=255)
                 gCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCa                                                                                                                                  >  1:1252308/1‑150 (MQ=255)
                       tGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAg                                                                                                                            <  2:869390/150‑1 (MQ=255)
                          cGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTaa                                                                                                                         >  1:2293883/1‑150 (MQ=255)
                          cGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTaa                                                                                                                         >  1:434923/1‑150 (MQ=255)
                          cGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTGATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCAAGAAAAGTaa                                                                                                                         >  2:829458/1‑150 (MQ=255)
                            ggCACTGATTTTTCTACTGACGGGCTGATTAGACGTTATGTCACTGGCGTTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAg                                                                                                                       <  2:2059036/150‑1 (MQ=255)
                            ggCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGGCCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAg                                                                                                                       <  2:1640017/150‑1 (MQ=255)
                            ggCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAg                                                                                                                       <  1:181855/150‑1 (MQ=255)
                             gCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGc                                                                                                                      <  2:415794/150‑1 (MQ=255)
                                cTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATg                                                                                                                   >  2:1108168/1‑150 (MQ=255)
                                      ttgcTACTGACGGGCTGAGGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGAGGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGaaa                                                                                                             <  1:1283481/147‑1 (MQ=255)
                                        tCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATg                                                                                                           <  1:223042/150‑1 (MQ=255)
                                           aCTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACg                                                                                                        <  1:275795/150‑1 (MQ=255)
                                           aCTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCACTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACg                                                                                                        <  2:2230913/150‑1 (MQ=255)
                                             tGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATTGCGAAATGACGTg                                                                                                      <  2:1703863/150‑1 (MQ=255)
                                             tGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTg                                                                                                      <  1:148707/150‑1 (MQ=255)
                                             tGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTg                                                                                                      <  1:698767/150‑1 (MQ=255)
                                              gACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACAATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGa                                                                                                     >  1:1086102/1‑150 (MQ=255)
                                                       atgaGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCg                                                                                            >  1:1803523/1‑150 (MQ=255)
                                                        tgaGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGa                                                                                           >  1:1126547/1‑150 (MQ=255)
                                                        tgaGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGa                                                                                           >  2:129048/1‑150 (MQ=255)
                                                          agaCGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAg                                                                                         <  2:441452/150‑1 (MQ=255)
                                                           gaCGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACCGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAgg                                                                                        >  2:1925277/1‑150 (MQ=255)
                                                               ttATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATAc                                                                                    <  1:1386092/150‑1 (MQ=255)
                                                                tATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATAcc                                                                                   >  2:245071/1‑150 (MQ=255)
                                                                 aTGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCa                                                                                  <  2:1086102/150‑1 (MQ=255)
                                                                     cACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCaa                                                                              >  2:2265714/1‑150 (MQ=255)
                                                                     cACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCaa                                                                              <  2:2307370/150‑1 (MQ=255)
                                                                     cACTGGCGGGAATAATCAGTGGCACACTGTAGCCTCCCTGAACATGACCATAATGATGAACAGGATTGGCCCTGTCGCTGACCGCCAGCCCAAGGCCAGAAAAGTAAAGCATGGAGAAATGACGTGAGTACCGGCGAAGGATAcctgcaa                                                                              >  2:1813234/1‑145 (MQ=255)
                                                                       cTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGc                                                                            >  2:1430170/1‑150 (MQ=255)
                                                                         ggCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTg                                                                          >  1:1169764/1‑150 (MQ=255)
                                                                            ggTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATTGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGcc                                                                       <  2:1285844/150‑1 (MQ=255)
                                                                            ggTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGcc                                                                       <  1:739335/150‑1 (MQ=255)
                                                                             gTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCa                                                                      >  2:2363569/1‑150 (MQ=255)
                                                                                 taatCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGaa                                                                  <  1:641995/150‑1 (MQ=255)
                                                                                 taatCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGaa                                                                  >  1:2367504/1‑150 (MQ=255)
                                                                                   atCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAg                                                                <  1:188979/150‑1 (MQ=255)
                                                                                   atCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAACCATGGCGAAATGACGTGAGTACCCGCGAAGGATACCATCAAAAAGCCCAAGaaaa                                                                >  1:1064979/1‑149 (MQ=255)
                                                                                        tGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATc                                                           >  1:975672/1‑150 (MQ=255)
                                                                                         ggCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGCCCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATcc                                                          >  1:360164/1‑150 (MQ=255)
                                                                                         ggCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATcc                                                          >  1:1727066/1‑150 (MQ=255)
                                                                                         ggCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATcc                                                          >  1:2085996/1‑150 (MQ=255)
                                                                                          gCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCa                                                         >  1:2015836/1‑150 (MQ=255)
                                                                                              acTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGttt                                                     <  2:2136185/150‑1 (MQ=255)
                                                                                              acTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGttt                                                     <  1:1874982/150‑1 (MQ=255)
                                                                                                tGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTAc                                                   <  2:34287/150‑1 (MQ=255)
                                                                                                  tAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTg                                                 >  2:23939/1‑150 (MQ=255)
                                                                                                  tAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTg                                                 <  1:831613/150‑1 (MQ=255)
                                                                                                   aGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCAAAGGACAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCTAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGa                                                >  2:9402/1‑150 (MQ=255)
                                                                                                    gCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGAt                                               >  2:745810/1‑150 (MQ=255)
                                                                                                     ccTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATg                                              <  1:1212138/150‑1 (MQ=255)
                                                                                                           tGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCg                                        <  1:2109841/150‑1 (MQ=255)
                                                                                                           tGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCg                                        >  2:802933/1‑150 (MQ=255)
                                                                                                             cacaTGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAg                                      <  1:1563694/150‑1 (MQ=255)
                                                                                                               caTGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACATCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGt                                    <  1:1933110/150‑1 (MQ=255)
                                                                                                                 tGACCATAATGATGAACAGGATTGGCACTGTCGCTGACTGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTaa                                  >  1:1610393/1‑150 (MQ=255)
                                                                                                                   aCCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACa                                >  2:1793871/1‑150 (MQ=255)
                                                                                                                   aCCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCAAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACa                                <  1:2476160/150‑1 (MQ=255)
                                                                                                                    ccATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAg                               >  2:594721/1‑150 (MQ=255)
                                                                                                                      aTAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGc                             <  1:1039019/150‑1 (MQ=255)
                                                                                                                      aTAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGc                             <  2:1648744/150‑1 (MQ=255)
                                                                                                                      aTAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGc                             >  1:783478/1‑150 (MQ=255)
                                                                                                                          tgctgaACATGCTTGGGAATGTCGCTGCACGACATGCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAAtt                         <  2:15563/147‑1 (MQ=255)
                                                                                                                          tgatgaACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAAc                                 >  2:2365611/1‑142 (MQ=255)
                                                                                                                            atgaACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAAtttt                       <  1:1267263/150‑1 (MQ=255)
                                                                                                                              gaACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTc                     >  1:1587828/1‑150 (MQ=255)
                                                                                                                                aCAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCgg                   <  2:1867687/150‑1 (MQ=255)
                                                                                                                                aCAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCgg                   <  1:2363569/150‑1 (MQ=255)
                                                                                                                                aCAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCgg                   >  1:2372171/1‑150 (MQ=255)
                                                                                                                                 cAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCggg                  <  2:523319/150‑1 (MQ=255)
                                                                                                                                   ggATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGCAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGa                >  2:1872833/1‑150 (MQ=255)
                                                                                                                                   ggATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAAtt                         >  1:1411034/1‑141 (MQ=255)
                                                                                                                                   ggATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAAtt                         <  2:1411034/141‑1 (MQ=255)
                                                                                                                                   ggATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGa                >  1:1539966/1‑150 (MQ=255)
                                                                                                                                     ccttgggcTGTCGGTCAGCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATa              >  2:1466924/8‑150 (MQ=255)
                                                                                                                                      cttgggcTGTCGGTCAGCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTACCAGGCAATTTTTCGGGGATAc             >  2:1265469/7‑150 (MQ=255)
                                                                                                                                        ggCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTg           >  1:275318/1‑150 (MQ=255)
                                                                                                                                         gCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGc          >  1:1035778/1‑150 (MQ=255)
                                                                                                                                          ggcTGTCGGTCAGCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCt         >  1:1906521/3‑150 (MQ=255)
                                                                                                                                          ggcTGTCGGTCAGCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCt         >  1:586851/3‑150 (MQ=255)
                                                                                                                                           aCTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTc        >  2:2084898/1‑150 (MQ=255)
                                                                                                                                                cGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGt   >  2:993265/1‑150 (MQ=255)
                                                                                                                                                 gCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTa  >  2:476006/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TCCACTGAAAAATACCTGCGAAATGACGGGCACTGATTTTTCTACTGACGGGCTGATGAGACGTTATGTCACTGGCGGTAATAATCAGTGGCACACTGTAGCCTCCCTGCACATGACCATAATGATGAACAGGATTGGCACTGTCGCTGACCGACAGCCCAAGGCCAGAAAAGTAAAGCATGGCGAAATGACGTGAGTACCGGCGAAGGATACCATCAAGCTGGCCAAGAAAGTTATCCAGTTTACTGATGCTGGCGAGGTAACAGGCAATTTTTCGGGGATACTGCTCCAGGTA  >  NC_000913/4499739‑4500033

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: