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,8910GA100.0% 313.5 / NA 95*1260* (TAG→TAAtamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (53/42);  total (53/42)

GTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCG  >  NC_000913/4447769‑4448023
                                                                                                                          |                                                                                                                                    
gTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTc                                                                                                                     <  1:227043/139‑1 (MQ=255)
 tGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCt                                                                                                                    >  1:348900/1‑139 (MQ=255)
             aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTAc                                                                                                        >  1:439395/1‑139 (MQ=255)
             aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGATCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTAc                                                                                                        >  1:1786914/1‑139 (MQ=255)
              ctctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACg                                                                                                       >  2:814787/1‑139 (MQ=255)
               tctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGc                                                                                                      >  1:28870/1‑139 (MQ=255)
                ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGcc                                                                                                     >  2:896358/1‑139 (MQ=255)
                ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGcc                                                                                                     <  2:28870/139‑1 (MQ=255)
                   taGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCcaca                                                                                                  >  2:252329/1‑139 (MQ=255)
                       aaCACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATAttt                                                                                                                        >  1:1689994/1‑113 (MQ=255)
                       aaCACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATAttt                                                                                                                        <  2:1689994/113‑1 (MQ=255)
                         cacTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAAt                                                                                                                                  <  2:749716/101‑1 (MQ=255)
                         cacTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAAt                                                                                                                                  >  1:749716/1‑101 (MQ=255)
                           cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGc                                                                                          >  1:1506486/1‑139 (MQ=255)
                           cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTCGCCACAAACAAGGCa                                                                                         <  2:348900/139‑1 (MQ=255)
                              acgttacgttATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAAc                                                                                       >  1:1805052/1‑139 (MQ=255)
                                gttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAg                                                                                     >  1:1857906/1‑139 (MQ=255)
                                        aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGa                                                                                                                              >  1:221105/1‑90 (MQ=255)
                                        aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGa                                                                                                                              <  2:221105/90‑1 (MQ=255)
                                        aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        >  1:81591/1‑127 (MQ=255)
                                        aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        <  2:81591/127‑1 (MQ=255)
                                         tCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGAt                                                                            <  2:1864496/139‑1 (MQ=255)
                                          cGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATg                                                                           >  1:1268738/1‑139 (MQ=255)
                                              tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAGTGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                                                    >  2:1297696/1‑126 (MQ=255)
                                              tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGt                                                                                    <  1:1297696/126‑1 (MQ=255)
                                              tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCa                                                                       >  2:830629/1‑139 (MQ=255)
                                              tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCa                                                                       >  1:871990/1‑139 (MQ=255)
                                               gATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                                                      <  2:1805052/139‑1 (MQ=255)
                                                aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCATAAACAAGGCAACAGTCACTGGATGACCaa                                                                      <  1:1086194/138‑1 (MQ=255)
                                                aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                                                      >  2:1086194/1‑138 (MQ=255)
                                                aTGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAAt                                                                     <  2:871990/139‑1 (MQ=255)
                                                            tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATg                                                                    <  1:880234/128‑1 (MQ=255)
                                                            tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATg                                                                    >  2:880234/1‑128 (MQ=255)
                                                             atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAAt                                                                                                                                  >  1:489188/1‑65 (MQ=255)
                                                             atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAAt                                                                                                                                  <  2:489188/65‑1 (MQ=255)
                                                               cTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGc                                                      <  2:1636987/139‑1 (MQ=255)
                                                                  gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAt                                                   >  2:1668088/1‑139 (MQ=255)
                                                                  gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAt                                                   >  2:1791326/1‑139 (MQ=255)
                                                                  gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAt                                                   >  2:1324640/1‑139 (MQ=255)
                                                                   aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaa                                                                                                >  2:337634/1‑93 (MQ=255)
                                                                   aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaa                                                                                                <  1:337634/93‑1 (MQ=255)
                                                                   aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        >  1:67319/1‑101 (MQ=255)
                                                                   aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        <  2:67319/101‑1 (MQ=255)
                                                                   aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAtt                                                  >  2:1003529/1‑139 (MQ=255)
                                                                   aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGAtt                                                  >  1:610857/1‑139 (MQ=255)
                                                                    aGCCCTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGGATATTTGTCTACGGCAGTTTACGCCCCAAAAAAGGCCACAATCACTGGATGACCAATGCCCAGTTACTCGGCGAttt                                                 >  1:584356/1‑139 (MQ=255)
                                                                      ccGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACa                                                                                                                                    >  2:1646444/1‑54 (MQ=255)
                                                                      ccGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACa                                                                                                                                    <  1:1646444/54‑1 (MQ=255)
                                                                      ccGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTa                                                                                                         >  2:145460/1‑81 (MQ=255)
                                                                      ccGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTa                                                                                                         <  1:145460/81‑1 (MQ=255)
                                                                      ccGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                                               <  1:1668088/139‑1 (MQ=255)
                                                                      ccGTGTCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                                               >  1:943427/1‑139 (MQ=255)
                                                                        gtgtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGttt                                                                                                          <  1:410058/78‑1 (MQ=255)
                                                                        gtgtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGttt                                                                                                          >  2:410058/1‑78 (MQ=255)
                                                                        gtgtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        >  2:273465/1‑96 (MQ=255)
                                                                        gtgtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        <  1:273465/96‑1 (MQ=255)
                                                                         tgtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAAt                                                                                                                                  >  2:1311370/1‑53 (MQ=255)
                                                                         tgtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAAt                                                                                                                                  <  1:1311370/53‑1 (MQ=255)
                                                                         tgtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTa                                            <  1:252329/139‑1 (MQ=255)
                                                                          gtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAatat                                                                                                                          >  1:217137/1‑60 (MQ=255)
                                                                          gtCTGGTGTAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAatat                                                                                                                          <  2:217137/60‑1 (MQ=255)
                                                                               gtgtAGACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa                                      <  2:610857/139‑1 (MQ=255)
                                                                                    gACCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTAc                                 >  1:1499718/1‑139 (MQ=255)
                                                                                     aCCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACa                                                                                                                                    >  2:720988/1‑39 (MQ=255)
                                                                                     aCCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACa                                                                                                                                    <  1:720988/39‑1 (MQ=255)
                                                                                     aCCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAatat                                                                                                                          >  2:1560089/1‑49 (MQ=255)
                                                                                     aCCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAatat                                                                                                                          <  1:1560089/49‑1 (MQ=255)
                                                                                     aCCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                                                      <  1:357020/101‑1 (MQ=255)
                                                                                     aCCAGGCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCaa                                                                      >  2:357020/1‑101 (MQ=255)
                                                                                         ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATg                                                                    <  2:772977/99‑1 (MQ=255)
                                                                                         ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATg                                                                    >  1:772977/1‑99 (MQ=255)
                                                                                         ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                                               <  1:629668/120‑1 (MQ=255)
                                                                                         ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCa                                               >  2:629668/1‑120 (MQ=255)
                                                                                         ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGtt                            >  2:1024934/1‑139 (MQ=255)
                                                                                         ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGtt                            >  1:1717931/1‑139 (MQ=255)
                                                                                         ggCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGtt                            >  1:1462430/1‑139 (MQ=255)
                                                                                          gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATAttt                                                                                                                        <  2:1485856/46‑1 (MQ=255)
                                                                                          gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATAttt                                                                                                                        >  1:1485856/1‑46 (MQ=255)
                                                                                          gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        <  2:487849/78‑1 (MQ=255)
                                                                                          gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCaa                                                                                        >  1:487849/1‑78 (MQ=255)
                                                                                          gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTc                                                                                   >  1:1231280/1‑83 (MQ=255)
                                                                                          gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTc                                                                                   <  2:1231280/83‑1 (MQ=255)
                                                                                          gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                                                  <  2:815422/84‑1 (MQ=255)
                                                                                          gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                                                  >  1:815422/1‑84 (MQ=255)
                                                                                          gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTg                           >  2:508311/1‑139 (MQ=255)
                                                                                          gCACTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTg                           >  1:231144/1‑139 (MQ=255)
                                                                                            aCTGGATTGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGta                         <  1:830629/139‑1 (MQ=255)
                                                                                                   ttGCTCTATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCt                   <  2:943427/139‑1 (MQ=255)
                                                                                                         tATCAGTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCa             <  2:1857906/139‑1 (MQ=255)
                                                                                                              gTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGt                          <  1:846025/121‑1 (MQ=255)
                                                                                                              gTTCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGt                          >  2:846025/1‑121 (MQ=255)
                                                                                                               ttCGAGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATcc       <  2:1499718/139‑1 (MQ=255)
                                                                                                                   aGTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgc   >  1:1703312/1‑139 (MQ=255)
                                                                                                                    gTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgcg  <  2:231144/139‑1 (MQ=255)
                                                                                                                    gTTTTAACAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgcg  <  1:1324640/139‑1 (MQ=255)
                                                                                                                          |                                                                                                                                    
GTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCG  >  NC_000913/4447769‑4448023

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