Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,447,853 G→C 10.2% D1248H (GAC→CAC)  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,8530GC10.2% 563.5 / 47.6 166D1248H (GAC→CAC) tamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base G (76/73);  new base C (8/9);  total (84/82)
Fisher's exact test for biased strand distribution p-value = 8.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.80e-01

Only 100 of 168 total aligned reads displayed.

CGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAAC  >  NC_000913/4447708‑4447989
                                                                                                                                                 |                                                                                                                                        
cGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGAc                                                                                                                                        >  2:692861/1‑148 (MQ=255)
cGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGAc                                                                                                                                        >  2:702346/1‑148 (MQ=255)
       tccCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCAc                                                                                                                                 <  1:702346/148‑1 (MQ=255)
       tccCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCAc                                                                                                                                 <  1:692861/148‑1 (MQ=255)
          cAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacacg                                                                                                                                <  4:878056/146‑2 (MQ=255)
          cAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacacg                                                                                                                                >  3:878056/1‑145 (MQ=255)
           aGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGa                                                                                                                             >  4:914379/1‑148 (MQ=255)
              tAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacac                                                                                                                                 >  4:423507/1‑141 (MQ=255)
                    tCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTa                                                                                                                    <  3:914379/148‑1 (MQ=255)
                                   tGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                            <  4:480045/141‑1 (MQ=255)
                                   tGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGACa                                                                                                                                  <  3:1123768/119‑1 (MQ=255)
                                                aGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTAGACCAGACa                                                                                                                                  >  2:681916/1‑106 (MQ=255)
                                                aGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGACa                                                                                                                                  <  1:681916/106‑1 (MQ=255)
                                                 gTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCa                                                                                                                                      >  3:200026/1‑101 (MQ=255)
                                                 gTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                       >  3:218250/1‑148 (MQ=255)
                                                              tGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGt                                                                          >  3:736843/1‑148 (MQ=255)
                                                                  tataTTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGa                                                                                            >  1:844040/1‑126 (MQ=255)
                                                                  tataTTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGa                                                                                            <  2:844040/126‑1 (MQ=255)
                                                                  tataTTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACAAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTa                                                                                                     >  3:1055405/1‑117 (MQ=255)
                                                                  tataTTTGACTCTATAGCAACACTCACGTTACGTTAACGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTa                                                                                                     <  4:1055405/117‑1 (MQ=255)
                                                                   atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                      >  1:956572/1‑115 (MQ=255)
                                                                   atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                      <  2:956572/115‑1 (MQ=255)
                                                                   atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacacggc                                                                                                                              >  1:1040807/1‑88 (MQ=38)
                                                                   atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacacggc                                                                                                                              <  2:1040807/91‑4 (MQ=38)
                                                                       ttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTa                                                                                                                    >  3:481736/1‑97 (MQ=255)
                                                                       ttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATg                                                                                               <  4:214155/118‑1 (MQ=255)
                                                                       ttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCcac                                                                 >  2:48998/1‑148 (MQ=255)
                                                                          aCTCTATAGCAACCCTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCa                                                                                                                                  <  3:423734/80‑1 (MQ=255)
                                                                          aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacac                                                                                                                                 <  2:386811/81‑1 (MQ=38)
                                                                          aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacac                                                                                                                                 >  1:386811/1‑81 (MQ=38)
                                                                             ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATc                                                                                                                  >  1:606869/1‑93 (MQ=255)
                                                                             ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaacaa                                                           <  1:48998/148‑1 (MQ=255)
                                                                             ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaacaa                                                           <  2:899195/148‑1 (MQ=255)
                                                                             ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaacaa                                                           >  1:899195/1‑148 (MQ=255)
                                                                             ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGATCTATc                                                                                                                  <  2:606869/93‑1 (MQ=255)
                                                                                  gCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                            >  4:729935/1‑94 (MQ=255)
                                                                                        cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACg                                                                               >  4:554530/1‑117 (MQ=255)
                                                                                        cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                            >  2:576791/1‑136 (MQ=255)
                                                                                        cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                            <  1:576791/136‑1 (MQ=255)
                                                                                        cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                >  3:1049410/1‑148 (MQ=255)
                                                                                          cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGct                                                                                                                       >  3:386530/1‑75 (MQ=255)
                                                                                          cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTc                                                 <  3:559198/145‑1 (MQ=255)
                                                                                          cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACt                                              >  1:735218/1‑148 (MQ=255)
                                                                                            cgttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAg                                                          >  1:944788/1‑134 (MQ=255)
                                                                                            cgttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAg                                                          <  2:944788/134‑1 (MQ=255)
                                                                                              ttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATg                                                                                               <  1:385274/95‑1 (MQ=255)
                                                                                              ttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATg                                                                                               >  2:385274/1‑95 (MQ=255)
                                                                                                acgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGa                                        <  2:735218/148‑1 (MQ=255)
                                                                                                    taTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTg                                                                                                                         <  2:571992/63‑1 (MQ=255)
                                                                                                    taTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTg                                                                                                                         >  1:571992/1‑63 (MQ=255)
                                                                                                     aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGTGAat                                                                                          <  3:772118/93‑1 (MQ=255)
                                                                                                         ccTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                            >  4:300599/1‑71 (MQ=255)
                                                                                                           tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                >  1:862535/1‑129 (MQ=255)
                                                                                                           tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                <  2:862535/129‑1 (MQ=255)
                                                                                                              tGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGa                                        >  1:82370/1‑134 (MQ=255)
                                                                                                              tGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGa                                        <  2:82370/134‑1 (MQ=255)
                                                                                                                ccTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACt                        >  3:594819/1‑148 (MQ=255)
                                                                                                                  tAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                      >  4:774241/1‑68 (MQ=255)
                                                                                                                  tAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGc                                                                             >  3:1066234/1‑93 (MQ=255)
                                                                                                                  tAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGc                                                                             <  4:1066234/93‑1 (MQ=255)
                                                                                                                   aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                       <  1:711779/82‑1 (MQ=255)
                                                                                                                   aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                       >  2:711779/1‑82 (MQ=255)
                                                                                                                   aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAAt                                   <  4:801526/134‑1 (MQ=255)
                                                                                                                   aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAAt                                   >  3:801526/1‑134 (MQ=255)
                                                                                                                   aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTggg                     <  4:594819/148‑1 (MQ=255)
                                                                                                                   aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTggg                     >  2:272109/1‑148 (MQ=255)
                                                                                                                    aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGa                                        <  1:328877/128‑1 (MQ=255)
                                                                                                                    aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGa                                        >  2:328877/1‑128 (MQ=255)
                                                                                                                     gCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                            <  1:1112369/59‑1 (MQ=255)
                                                                                                                     gCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                            >  2:1112369/1‑59 (MQ=255)
                                                                                                                     gCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTg                                                                                     >  4:277368/1‑82 (MQ=255)
                                                                                                                     gCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTg                                                                                     <  3:277368/82‑1 (MQ=255)
                                                                                                                       tataTCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                 >  2:130657/1‑84 (MQ=255)
                                                                                                                       tataTCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                 <  1:130657/84‑1 (MQ=255)
                                                                                                                          atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                               <  3:422475/115‑1 (MQ=255)
                                                                                                                          atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCAc                                               >  4:422475/1‑115 (MQ=255)
                                                                                                                          atCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTAACTGGGcatt                 <  4:133464/145‑4 (MQ=255)
                                                                                                                            cTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                            <  3:120906/100‑1 (MQ=255)
                                                                                                                            cTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                            >  4:120906/1‑100 (MQ=255)
                                                                                                                              ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAAc                                                     <  3:321116/105‑1 (MQ=255)
                                                                                                                              ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTa          <  4:1049410/148‑1 (MQ=255)
                                                                                                                               gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAat                                                                                          <  1:31925/67‑1 (MQ=255)
                                                                                                                               gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAat                                                                                          >  2:31925/1‑67 (MQ=255)
                                                                                                                               gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                            <  4:81599/97‑1 (MQ=255)
                                                                                                                               gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCg                   <  3:778222/138‑1 (MQ=255)
                                                                                                                               gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt         >  4:76422/1‑148 (MQ=255)
                                                                                                                                aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                 >  4:586787/1‑75 (MQ=255)
                                                                                                                                    cGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa    <  1:272109/148‑1 (MQ=255)
                                                                                                                                     gtgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTAACTGGGca                   >  2:797122/1‑131 (MQ=255)
                                                                                                                                     gtgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTAACTGGGca                   <  1:797122/132‑2 (MQ=255)
                                                                                                                                      tgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCTACAAACAAGGCAACAGTCACTGGATGa                                        >  3:739957/1‑110 (MQ=255)
                                                                                                                                      tgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa    <  2:593717/146‑1 (MQ=255)
                                                                                                                                      tgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa    >  1:593717/1‑146 (MQ=255)
                                                                                                                                      tggcTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTTACTGGATGa                                        <  4:739957/107‑1 (MQ=255)
                                                                                                                                        tCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACt                        <  4:512124/124‑1 (MQ=255)
                                                                                                                                        tCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACt                        >  3:512124/1‑124 (MQ=255)
                                                                                                                                         cTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTc              >  4:59737/1‑133 (MQ=255)
                                                                                                                                          tGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTgg                      >  4:330547/1‑124 (MQ=255)
                                                                                                                                              gtAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAAc                                                     >  3:593575/1‑89 (MQ=255)
                                                                                                                                               tAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACt                        <  3:634058/117‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                        
CGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAAC  >  NC_000913/4447708‑4447989

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