Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,447,844 T→C 11.8% S1245P (TCT→CCT)  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,8440TC11.8% 460.6 / 41.7 153S1245P (TCT→CCT) tamBtranslocation and assembly module for autotransporter export, inner membrane subunit
Reads supporting (aligned to +/- strand):  ref base T (68/67);  new base C (9/9);  total (77/76)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.58e-01

Only 100 of 161 total aligned reads displayed.

CGACTCCTCCCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAAC  >  NC_000913/4447708‑4447989
                                                                                                                                        |                                                                                                                                                 
       tccCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCAc                                                                                                                                 <  1:702346/148‑1 (MQ=255)
       tccCAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCAc                                                                                                                                 <  1:692861/148‑1 (MQ=255)
          ctggTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTg                                                                                                                                                >  4:572680/3‑130 (MQ=255)
          cAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacacg                                                                                                                                <  4:878056/146‑2 (MQ=255)
          cAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacacg                                                                                                                                >  3:878056/1‑145 (MQ=255)
          cAGGTAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGACa                                                                                                                                  >  2:323178/1‑144 (MQ=255)
              tAGTGGTCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacac                                                                                                                                 <  3:423507/141‑1 (MQ=255)
                    tCAGCGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTa                                                                                                                    <  3:914379/148‑1 (MQ=255)
                        cGGCTATGTATTGCCAGGTCTGCAAGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAgac                                                                                                                                   >  2:698403/1‑127 (MQ=255)
                                                aGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTAGACCAGACa                                                                                                                                  >  2:681916/1‑106 (MQ=255)
                                                aGTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGACa                                                                                                                                  <  1:681916/106‑1 (MQ=255)
                                                 gTGAAATACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCa                                                                                                                                      >  3:200026/1‑101 (MQ=255)
                                                       tACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                            <  1:403517/121‑1 (MQ=255)
                                                       tACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                            >  2:403517/1‑121 (MQ=255)
                                                       tACGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                 >  2:841085/1‑148 (MQ=255)
                                                        aCGGCGTGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGtt                                                                                                              >  4:908533/1‑118 (MQ=255)
                                                              tGGGTATATTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGt                                                                          >  3:736843/1‑148 (MQ=255)
                                                                  tataTTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGa                                                                                            >  1:844040/1‑126 (MQ=255)
                                                                  tataTTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACAAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTa                                                                                                     >  3:1055405/1‑117 (MQ=255)
                                                                  tataTTTGACTCTATAGCAACACTCACGTTACGTTAACGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTa                                                                                                     <  4:1055405/117‑1 (MQ=255)
                                                                   atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                      >  3:953931/1‑115 (MQ=255)
                                                                   atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                      <  2:956572/115‑1 (MQ=255)
                                                                   atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                      <  4:953931/115‑1 (MQ=255)
                                                                   atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacacggc                                                                                                                              >  1:1040807/1‑88 (MQ=38)
                                                                   atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacacggc                                                                                                                              <  2:1040807/91‑4 (MQ=38)
                                                                   atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGACa                                                                                                                                  >  2:84665/1‑87 (MQ=255)
                                                                   atatTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGACa                                                                                                                                  <  1:84665/87‑1 (MQ=255)
                                                                     atTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACt                                                                                                                                >  1:590754/1‑87 (MQ=255)
                                                                     atTTGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACt                                                                                                                                <  2:590754/87‑1 (MQ=255)
                                                                       ttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTa                                                                                                                    <  4:481736/97‑1 (MQ=255)
                                                                       ttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTa                                                                                                                    >  3:481736/1‑97 (MQ=255)
                                                                       ttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATg                                                                                               <  4:214155/118‑1 (MQ=255)
                                                                       ttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATg                                                                                               >  3:214155/1‑118 (MQ=255)
                                                                       ttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCcac                                                                 >  2:602983/1‑148 (MQ=255)
                                                                       ttGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCcac                                                                 >  2:48998/1‑148 (MQ=255)
                                                                        tGACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCAc                                                                                                                                 >  2:728201/1‑83 (MQ=255)
                                                                         gACTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaa                                                               <  4:736843/148‑1 (MQ=255)
                                                                          aCTCTATAGCAACCCTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCa                                                                                                                                  <  3:423734/80‑1 (MQ=255)
                                                                          aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCa                                                                                                                                  >  4:423734/1‑80 (MQ=255)
                                                                          aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCa                                                                                                                                  >  1:375631/1‑80 (MQ=255)
                                                                          aCTCTATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGCCTGGTCTACACCAGacac                                                                                                                                 >  1:386811/1‑81 (MQ=38)
                                                                             ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGc                                                                                                                                   <  2:912234/76‑1 (MQ=255)
                                                                             ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATc                                                                                                                  >  1:606869/1‑93 (MQ=255)
                                                                             ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGa                                                                                            <  1:948350/115‑1 (MQ=255)
                                                                             ctATAGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaacaa                                                           >  1:899195/1‑148 (MQ=255)
                                                                               ataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAat                                                                                          >  1:207944/1‑115 (MQ=255)
                                                                               ataGCAACACTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAat                                                                                          <  2:207944/115‑1 (MQ=255)
                                                                                        cTCACGTTACGTTATCGCCTGATGGCTAAGCTATATCTGGCAGCCGTGTCTGGTGGAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGATTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                >  3:112268/1‑148 (MQ=255)
                                                                                        cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACg                                                                               <  3:554530/117‑1 (MQ=255)
                                                                                        cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                            >  2:576791/1‑136 (MQ=255)
                                                                                        cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                >  3:1049410/1‑148 (MQ=255)
                                                                                        cTCACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                <  2:119894/148‑1 (MQ=255)
                                                                                          cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGct                                                                                                                       >  3:386530/1‑75 (MQ=255)
                                                                                          cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGct                                                                                                                       <  4:386530/75‑1 (MQ=255)
                                                                                          cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTc                                                 >  4:559198/1‑145 (MQ=255)
                                                                                          cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTc                                                 <  3:559198/145‑1 (MQ=255)
                                                                                          cACGTTACGTTATCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACt                                              >  1:735218/1‑148 (MQ=255)
                                                                                            cgttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAg                                                          <  2:944788/134‑1 (MQ=255)
                                                                                            cgttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAg                                                          >  1:944788/1‑134 (MQ=255)
                                                                                              ttacgttaTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATg                                                                                               >  2:385274/1‑95 (MQ=255)
                                                                                                    taTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTg                                                                                                                         >  1:571992/1‑63 (MQ=255)
                                                                                                     aTCGCCTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGTGAat                                                                                          <  3:772118/93‑1 (MQ=255)
                                                                                                        gCCGGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                       >  2:686188/1‑93 (MQ=255)
                                                                                                         ccTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                            <  3:300599/71‑1 (MQ=255)
                                                                                                         ccTGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                            >  4:300599/1‑71 (MQ=255)
                                                                                                           tGATGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCa                                                <  2:862535/129‑1 (MQ=255)
                                                                                                              tGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGa                                        >  1:82370/1‑134 (MQ=255)
                                                                                                                ccTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACt                        >  3:594819/1‑148 (MQ=255)
                                                                                                                  tAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGtttt                                                                                                      >  4:774241/1‑68 (MQ=255)
                                                                                                                  tAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGc                                                                             <  4:1066234/93‑1 (MQ=255)
                                                                                                                  tAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGc                                                                             >  3:1066234/1‑93 (MQ=255)
                                                                                                                   aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                       >  2:711779/1‑82 (MQ=255)
                                                                                                                   aaGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATAtt                                                                                       <  1:711779/82‑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                                        >  2:328877/1‑128 (MQ=255)
                                                                                                                    aGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGa                                        <  1:328877/128‑1 (MQ=255)
                                                                                                                     gCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCg                                                                                                            <  1:1112369/59‑1 (MQ=255)
                                                                                                                     gCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTg                                                                                     <  3:277368/82‑1 (MQ=255)
                                                                                                                     gCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTg                                                                                     >  4:277368/1‑82 (MQ=255)
                                                                                                                       tataTCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                 <  1:130657/84‑1 (MQ=255)
                                                                                                                       tataTCTGGAAGCCGTGTCTAGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCa                                                       <  1:998951/110‑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)
                                                                                                                              ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTa          <  4:1049410/148‑1 (MQ=255)
                                                                                                                               gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAat                                                                                          <  1:31925/67‑1 (MQ=255)
                                                                                                                               gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCacaaaca                                                            <  4:81599/97‑1 (MQ=255)
                                                                                                                               gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCg                   <  3:778222/138‑1 (MQ=255)
                                                                                                                               gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCg                   >  4:778222/1‑138 (MQ=255)
                                                                                                                               gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt         >  4:76422/1‑148 (MQ=255)
                                                                                                                                aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                 >  2:757407/1‑75 (MQ=255)
                                                                                                                                aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                 <  1:757407/75‑1 (MQ=255)
                                                                                                                                aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTa                                                                                 <  3:586787/75‑1 (MQ=255)
                                                                                                                                aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATc        <  3:76422/148‑1 (MQ=255)
                                                                                                                                     gtgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTAACTGGGca                   >  2:797122/1‑131 (MQ=255)
                                                                                                                                      tgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATa    >  1:593717/1‑146 (MQ=255)
                                                                                                                                        tCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACt                        <  4:512124/124‑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: