Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,447,891 G→A *1260* (TAG→TAA tamB → translocation and assembly module for autotransporter export, inner membrane subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,447,8910GA97.5% 244.2 / ‑4.7 79*1260* (TAG→TAAtamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (1/1);  new base A (41/36);  total (42/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCG  >  NC_000913/4447754‑4448023
                                                                                                                                         |                                                                                                                                    
cAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACa                                                                                                                                    >  2:644738/1‑139 (MQ=255)
               gTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTc                                                                                                                     >  1:664880/1‑139 (MQ=255)
                          tGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGttt                                                                                                          <  2:100868/139‑1 (MQ=255)
                            aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCg                                                                                                                               <  1:368292/116‑1 (MQ=255)
                            aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCg                                                                                                                               >  2:368292/1‑116 (MQ=255)
                             ctctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACg                                                                                                       <  2:291349/139‑1 (MQ=255)
                                 ataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCcac                                                                                                   >  1:582498/1‑139 (MQ=255)
                                          cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                                          >  1:480118/1‑139 (MQ=255)
                                          cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAgg                                                                                           <  2:93023/138‑1 (MQ=255)
                                          cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAgg                                                                                           >  1:93023/1‑138 (MQ=255)
                                          cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                                          >  2:265666/1‑139 (MQ=255)
                                                ttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                                                    >  1:60401/1‑139 (MQ=255)
                                                  acgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACgg                                                                                                                >  2:499014/1‑109 (MQ=255)
                                                  acgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACgg                                                                                                                <  1:499014/109‑1 (MQ=255)
                                                       aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGa                                                                             >  2:200357/1‑139 (MQ=255)
                                                              gATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                                                      <  1:615334/139‑1 (MQ=255)
                                                               aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                                                      >  1:9595/1‑138 (MQ=255)
                                                               aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                                                      <  2:9595/138‑1 (MQ=255)
                                                                  ccTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGcc                                                                  >  2:609614/1‑139 (MQ=255)
                                                                    tAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACACCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCa                                                                <  2:664880/139‑1 (MQ=255)
                                                                           tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                                            >  1:235175/1‑120 (MQ=255)
                                                                           tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                                            <  2:235175/120‑1 (MQ=255)
                                                                           tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTg                                                         >  1:271194/1‑139 (MQ=255)
                                                                            atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTgg                                                        >  2:357795/1‑139 (MQ=255)
                                                                              cTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                                                      <  2:582498/139‑1 (MQ=255)
                                                                              cTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                                                      <  1:609614/139‑1 (MQ=255)
                                                                                ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGa                                                    >  2:799504/1‑139 (MQ=255)
                                                                                 gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGc                                                                                                      >  2:152778/1‑88 (MQ=255)
                                                                                 gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                                                 >  2:442393/1‑109 (MQ=255)
                                                                                 gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                                                                 <  1:442393/109‑1 (MQ=255)
                                                                                 gAAGCCGTGTCTGGTGTAGACCAGGAACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGc                                                                                                      <  1:152778/88‑1 (MQ=255)
                                                                                  aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTa                                                                                                                   <  2:827234/74‑1 (MQ=255)
                                                                                  aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTa                                                                                                                   >  1:827234/1‑74 (MQ=255)
                                                                                  aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        >  2:194384/1‑101 (MQ=255)
                                                                                  aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        <  1:194384/101‑1 (MQ=255)
                                                                                      cGTGTCTGGTGTAGACCAGGCACTGGATTGTTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTc                                                                                                                     >  1:758476/1‑68 (MQ=255)
                                                                                      cGTGTCTGGTGTAGACCAGGCACTGGATTGTTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTc                                                                                                                     <  2:758476/68‑1 (MQ=255)
                                                                                        tgtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                                               <  1:209563/136‑1 (MQ=255)
                                                                                        tgtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                                               >  2:209563/1‑136 (MQ=255)
                                                                                                    aCCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACg                                                                                                                 <  2:243289/58‑1 (MQ=255)
                                                                                                    aCCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACg                                                                                                                 >  1:243289/1‑58 (MQ=255)
                                                                                                    aCCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAcc                                >  1:443348/1‑139 (MQ=255)
                                                                                                     ccAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCa                               <  1:357795/139‑1 (MQ=255)
                                                                                                      cAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAg                              <  2:686747/139‑1 (MQ=255)
                                                                                                        ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAatat                                                                                                                          <  1:397152/45‑1 (MQ=255)
                                                                                                        ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAatat                                                                                                                          >  2:397152/1‑45 (MQ=255)
                                                                                                        ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGc                                                                                                      >  1:235459/1‑65 (MQ=255)
                                                                                                        ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGc                                                                                                      <  2:235459/65‑1 (MQ=255)
                                                                                                        ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGcta                                                                                        >  1:837696/1‑77 (MQ=255)
                                                                                                        ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        <  1:563794/79‑1 (MQ=255)
                                                                                                        ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        >  2:563794/1‑79 (MQ=255)
                                                                                                        ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        <  2:837696/79‑1 (MQ=255)
                                                                                                        ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTggg                                                       >  2:743405/1‑112 (MQ=255)
                                                                                                        ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTggg                                                       <  1:743405/112‑1 (MQ=255)
                                                                                                        ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCa                               >  1:381289/1‑136 (MQ=255)
                                                                                                        ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCa                               <  2:381289/136‑1 (MQ=255)
                                                                                                         gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTa                                                                                                                   <  1:268842/51‑1 (MQ=255)
                                                                                                         gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTa                                                                                                                   >  2:268842/1‑51 (MQ=255)
                                                                                                         gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                                                    >  1:130113/1‑82 (MQ=255)
                                                                                                         gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                                                    <  2:130113/82‑1 (MQ=255)
                                                                                                         gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGa                                        <  2:148895/126‑1 (MQ=255)
                                                                                                         gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGa                                        >  1:148895/1‑126 (MQ=255)
                                                                                                              tggattGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa                                      >  1:776652/5‑124 (MQ=255)
                                                                                                              tggattGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa                                      <  2:776652/120‑1 (MQ=255)
                                                                                                              tggattGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGtata                       >  1:539664/5‑139 (MQ=255)
                                                                                                              tggattGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGtata                       >  2:526258/5‑139 (MQ=255)
                                                                                                               ggattGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        <  1:24254/70‑1 (MQ=255)
                                                                                                               ggattGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        <  2:786664/70‑1 (MQ=255)
                                                                                                               ggattGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        >  1:786664/4‑73 (MQ=255)
                                                                                                               ggattGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        >  2:24254/4‑73 (MQ=255)
                                                                                                                     ctctATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                                                      <  2:275965/101‑1 (MQ=255)
                                                                                                                     ctctATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                                                      >  1:275965/1‑101 (MQ=255)
                                                                                                                     ctctATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAg                              >  2:479774/1‑125 (MQ=255)
                                                                                                                     ctctATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAg                              <  1:479774/125‑1 (MQ=255)
                                                                                                                         aTCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCAc            <  1:200357/139‑1 (MQ=255)
                                                                                                                               tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCa      >  1:795847/1‑139 (MQ=255)
                                                                                                                               tCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCa      >  2:837437/1‑139 (MQ=255)
                                                                                                                                   gTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgcg  <  2:480118/139‑1 (MQ=255)
                                                                                                                                   gTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgcg  <  1:265666/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                    
CAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCG  >  NC_000913/4447754‑4448023

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