Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 918,118 G→A 6.0% G549D (GGT→GAT) ‡ ybjD → putative OLD family ATP‑dependent endonuclease; DUF2813 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913918,1180GA 6.0% 263.7 / 5.4 83G549D (GGT→GAT) ‡ybjDputative OLD family ATP‑dependent endonuclease; DUF2813 family protein
Reads supporting (aligned to +/- strand):  ref base G (43/35);  new base A (2/3);  total (45/38)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

ACGCAAAATTATCTCGAAAGCGATCCATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCTTTTCGTGCTATCTGTGCTGGAAGGCGGTAATG  >  NC_000913/917972‑918262
                                                                                                                                                  |                                                                                                                                                
aCGCAAAATTATCTCGAAAGCGATCCATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTc                                                                                                                                                <  2:263652/149‑1 (MQ=255)
aCGCAAAATTATCTCGAAAGCGATCCATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTc                                                                                                                                                <  2:182331/149‑1 (MQ=255)
  gCAAAATTATCTCGAAAGCGATCCATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTcgc                                                                                                                                              >  2:215061/1‑149 (MQ=255)
  gCAAAATTATCTCGAAAGCGATCCATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTcgc                                                                                                                                              >  1:742143/1‑149 (MQ=255)
    aaaaTTATCTCGAAAGCGATCCATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTcgcgc                                                                                                                                            >  4:838044/1‑149 (MQ=255)
      aaTTATCTCGAAAGCGATCCATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCgg                                                                                                                                          <  4:22675/149‑1 (MQ=255)
       aTTATCTCGAAAGCGATCCATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGa                                                                                                                                         <  3:108144/149‑1 (MQ=255)
       aTTATCTCGAAAGCGATCCATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGa                                                                                                                                         <  2:669775/149‑1 (MQ=255)
       aTTATCTCGAAAGCGATCCATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGACcgcgcgcc                                                                                                                                         <  1:787162/149‑3 (MQ=255)
        ttATCTCGAAAGCGATCCATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGACcgcgcgcca                                                                                                                                        <  3:378495/149‑4 (MQ=255)
        ttATCTCGAAAGCGATCCATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGACcgcgcgcca                                                                                                                                        <  4:671954/149‑4 (MQ=255)
            ctcGAAAGCGATCCATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAAc                                                                                                                                    <  4:376344/149‑1 (MQ=255)
               gAAAGCGATCCATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACcgc                                                                                                                                 >  2:231578/1‑149 (MQ=255)
               gAAAGCGATCCATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACcgc                                                                                                                                 <  2:430992/149‑1 (MQ=255)
                      aTCCATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACAt                                                                                                                          >  1:432886/1‑149 (MQ=255)
                       tCCATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATc                                                                                                                         >  1:188985/1‑149 (MQ=255)
                        ccATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCg                                                                                                                        <  4:686310/149‑1 (MQ=255)
                         cATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGt                                                                                                                       <  1:270887/149‑1 (MQ=255)
                          aTCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGt                                                                                                                       <  4:659999/148‑1 (MQ=255)
                          aTCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGt                                                                                                                       >  3:659999/1‑148 (MQ=255)
                          aTCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTg                                                                                                                      <  4:495541/149‑1 (MQ=255)
                          aTCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTg                                                                                                                      <  3:229649/149‑1 (MQ=255)
                          aTCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACACTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTg                                                                                                                      <  1:191651/149‑1 (MQ=255)
                          aTCGCTCTTCCAAACCCGATCTGGCAATAGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGAGTAACCGCGAAACATCGTg                                                                                                                      >  2:42026/1‑149 (MQ=255)
                              ctctTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCa                                                                                                                  >  3:239731/1‑149 (MQ=255)
                                 ttCCATACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGAGCGGATTAACCGCGAAACATCGTGGCCAttt                                                                                                               >  1:255616/1‑149 (MQ=255)
                                   ccAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTtgt                                                                                                             <  2:188985/149‑1 (MQ=255)
                                   ccAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTtgt                                                                                                             >  1:799740/1‑149 (MQ=255)
                                    cAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTtgtg                                                                                                            >  3:54384/1‑149 (MQ=255)
                                    cAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTtgtg                                                                                                            >  1:566074/1‑149 (MQ=255)
                                    cAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTtgtg                                                                                                            <  3:212063/149‑1 (MQ=255)
                                       aCCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCt                                                                                                         <  2:566074/149‑1 (MQ=255)
                                          cGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGaa                                                                                                      <  4:786802/149‑1 (MQ=255)
                                           gATCTTGCCATTGAAGTGGCAATGGAGACAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAAt                                                                                                     <  2:215533/149‑1 (MQ=255)
                                               ttGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCg                                                                                                 >  3:660518/1‑149 (MQ=255)
                                                 gCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTc                                                                                               <  4:696107/149‑1 (MQ=255)
                                                   cATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGa                                                                                             >  3:480642/1‑149 (MQ=255)
                                                     ttGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGc                                                                                           >  4:333101/1‑149 (MQ=255)
                                                        aaGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCAtc                                                                                        <  3:430292/149‑1 (MQ=255)
                                                           tGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGCCGAGCATCTCa                                                                                     <  1:466624/149‑1 (MQ=255)
                                                               aaTGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAgc                                                                                 <  4:783921/149‑1 (MQ=255)
                                                                      gCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTg                                                                                                                      >  2:824768/1‑105 (MQ=255)
                                                                      gCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTg                                                                                                                      <  1:824768/105‑1 (MQ=255)
                                                                      gCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACg                                                                          <  1:660875/149‑1 (MQ=255)
                                                                      gCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACg                                                                          <  4:480642/149‑1 (MQ=255)
                                                                             gtcgtGGTGTGGACTCCGTACCGACGCTGCTGAAAAAATTGTTCTCACGCGTGCTGTGGCTGCCGCGCGGTCGCGCGGTTTAACGGCGAAACATGGTGCCCTTTGGTGGCTGAATAGGCTCGAGCATCTCATTGCGCCGACGGTATTcc                                                                   >  2:546602/1‑148 (MQ=255)
                                                                               cgtGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGc                                                                 <  2:584663/149‑1 (MQ=255)
                                                                                 tggtgTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGccc                                                               >  3:597964/1‑149 (MQ=255)
                                                                                   gtgtGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGt                                                             >  2:854025/1‑149 (MQ=255)
                                                                                         aCTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTa                                                       >  1:744067/1‑149 (MQ=255)
                                                                                          cTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTAc                                                      >  3:143907/1‑149 (MQ=255)
                                                                                           tCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACt                                                     >  4:331916/1‑149 (MQ=255)
                                                                                           tCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACt                                                     >  4:845978/1‑149 (MQ=255)
                                                                                               tACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCg                                                 >  4:568521/1‑149 (MQ=255)
                                                                                                aCCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGa                                                >  1:332900/1‑149 (MQ=255)
                                                                                                 ccGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGAt                                               >  1:861399/1‑149 (MQ=255)
                                                                                                 ccGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGAt                                               >  3:841329/1‑149 (MQ=255)
                                                                                                       ctgctgAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGc                                         >  4:658659/1‑149 (MQ=255)
                                                                                                        tgctgAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCa                                        >  2:442077/1‑149 (MQ=255)
                                                                                                        tgctgAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCa                                        >  3:482460/1‑149 (MQ=255)
                                                                                                             aaaaaaaTGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAAtctc                                   >  2:473924/1‑149 (MQ=255)
                                                                                                              aaaaaaTGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATctct                                  >  4:643958/1‑149 (MQ=255)
                                                                                                               aaaaaTGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCtt                                 >  2:231734/1‑149 (MQ=255)
                                                                                                               aaaaaTGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCtt                                 >  1:287102/1‑149 (MQ=255)
                                                                                                               aaaaaTGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCtt                                 >  2:613625/1‑149 (MQ=255)
                                                                                                                  aaTGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGTGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCGCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCTTTTc                              >  4:859483/1‑149 (MQ=255)
                                                                                                                  aaTGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCTTTTc                              >  3:277827/1‑149 (MQ=255)
                                                                                                                    tGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGcc                                                                              <  1:499232/99‑1 (MQ=255)
                                                                                                                    tGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGcc                                                                              >  2:499232/1‑99 (MQ=255)
                                                                                                                    tGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCTTTTCGt                            <  1:231578/149‑1 (MQ=255)
                                                                                                                     gTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCTTTTCGTg                           <  3:838044/149‑1 (MQ=255)
                                                                                                                      ttCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCTTTTCGTGc                          >  3:33568/1‑149 (MQ=255)
                                                                                                                             gcgTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCTTTTCGTGCTATCtgt                   <  1:353464/149‑1 (MQ=255)
                                                                                                                                 gCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCTTTTCGTGCTATCTGTGCTg               <  1:802112/149‑1 (MQ=255)
                                                                                                                                   tgtgGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCTTTTCGTGCTATCTGTGCTGGa             <  2:799740/149‑1 (MQ=255)
                                                                                                                                      ggctggcGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCTTTTCGTGCTATCTGTGCTGGAAgg          <  3:188209/149‑1 (MQ=255)
                                                                                                                                       gctggcGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCTTTTCGTGCTATCTGTGCTGGAAGGc         <  3:333101/149‑1 (MQ=255)
                                                                                                                                         tggcGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCTTTTCGTGCTATCTGTGCTGGAAGGCgg       >  1:398037/1‑149 (MQ=255)
                                                                                                                                          tccgcgcgACCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCTTTTCGTGCTATCTGTGCTGGAAGGCGGt      >  4:321628/3‑149 (MQ=255)
                                                                                                                                          tccgcgcgACCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCTTTTCGTGCTATCTGTGCTGGAAGGCGGt      >  3:708374/3‑149 (MQ=255)
                                                                                                                                          ggcGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCTTTTCGTGCTATCTGTGCTGGAAGGCGGt      >  3:756327/1‑149 (MQ=255)
                                                                                                                                           gcgcgcGGTCGCGCGGATTAACCGCGAGACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCTTTTCGTGCTATCTGTGCTGGAAGGCGGTa     >  2:105461/1‑149 (MQ=255)
                                                                                                                                              cgcgGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCTTTTCGTGCTATCTGTGCTGGAAGGCGGTAATg  <  1:127879/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
ACGCAAAATTATCTCGAAAGCGATCCATCGCTCTTCCAAACCCGATCTTGCCATTGAAGTGGCAATGGAGGCAGGACGTCGTGGTGTGGACTCCGTACCGACGCTGCTGAAAAAAATGTTCTCACGCGTGCTGTGGCTGGCGCGCGGTCGCGCGGATTAACCGCGAAACATCGTGGCCATTTGTGGCTGAATAGCGTCGAGCATCTCATAGCGCCGACGGTATTCAGCCCGTTTTTTACTGGCGATTTCGGCAATCTCTTTTCGTGCTATCTGTGCTGGAAGGCGGTAATG  >  NC_000913/917972‑918262

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: