Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,522,021 G→A 5.3% R153* (CGA→TGA)  xapR ← transcriptional activator of xapAB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,522,0210GA 5.3% 249.9 / 3.0 76R153* (CGA→TGA) xapRtranscriptional activator of xapAB
Reads supporting (aligned to +/- strand):  ref base G (32/40);  new base A (2/2);  total (34/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

ATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGAATCGCCGCATAACCGGGCGCATCCGTCCCCACATTGCCGTTCCCA  >  NC_000913/2521876‑2522167
                                                                                                                                                 |                                                                                                                                                  
aTTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCa                                                                                                                                                 <  2:632673/149‑1 (MQ=255)
aTTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCa                                                                                                                                                 <  3:196652/149‑1 (MQ=255)
   gTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGAt                                                                                                                                              <  1:251933/149‑1 (MQ=255)
      aCAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCACCAGATccc                                                                                                                                           >  2:21569/1‑149 (MQ=255)
      aCAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGACGGTGGTTCTGTCGCCATTCACCAGATccc                                                                                                                                           >  2:92284/1‑149 (MQ=255)
       cAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCg                                                                                                                                          <  1:134669/149‑1 (MQ=255)
        aaaaTACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCgg                                                                                                                                         >  4:216911/1‑149 (MQ=255)
           aTACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCAt                                                                                                                                      <  1:405715/149‑1 (MQ=255)
                cGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGt                                                                                                                                 <  2:642044/149‑1 (MQ=255)
                          cgcTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGtt                                                                                                                       >  2:692242/1‑149 (MQ=255)
                           gcTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTc                                                                                                                      >  4:491063/1‑149 (MQ=255)
                           gcTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTc                                                                                                                      <  3:143532/149‑1 (MQ=255)
                             ttCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCa                                                                                                                    >  2:206842/1‑149 (MQ=255)
                              tCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAg                                                                                                                   <  1:353093/149‑1 (MQ=255)
                                 aGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCaa                                                                                                                <  3:31564/149‑1 (MQ=255)
                                    gggACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGc                                                                                                             >  2:454157/1‑149 (MQ=255)
                                         ggTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCAttt                                                                                                        <  2:321799/149‑1 (MQ=255)
                                            ggAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCaa                                                                                                                                   <  4:465911/119‑1 (MQ=255)
                                            ggAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCaa                                                                                                                                   >  3:465911/1‑119 (MQ=255)
                                                 atgatgAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATcgcg                                                                                                <  2:330600/149‑1 (MQ=255)
                                                  tgatgaGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCgg                                                                                               <  4:443994/149‑1 (MQ=255)
                                                   gatgaGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCggg                                                                                              <  1:692242/149‑1 (MQ=255)
                                                        gagaTGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCt                                                                                         >  2:18445/1‑149 (MQ=255)
                                                         agaTGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCtt                                                                                        >  1:371469/1‑149 (MQ=255)
                                                            tgatgCTCTTCAGGCATCGCCACCAGAAACGTCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGAGCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTc                                                                                     <  4:26425/149‑1 (MQ=255)
                                                            tgatgCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTc                                                                                     <  2:143892/149‑1 (MQ=255)
                                                            tgatgCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTc                                                                                     <  3:216911/149‑1 (MQ=255)
                                                            tgatgCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTc                                                                                     <  4:581254/149‑1 (MQ=255)
                                                             gatgCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTcg                                                                                    >  3:253485/1‑149 (MQ=255)
                                                              atgCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTcgc                                                                                   >  1:540940/1‑149 (MQ=255)
                                                                  tctTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGaaa                                                                               >  3:232305/1‑149 (MQ=255)
                                                                  tctTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGaaa                                                                               >  3:630109/1‑149 (MQ=255)
                                                                     tCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAg                                                                            >  1:674690/1‑149 (MQ=255)
                                                                        ggCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCACCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAAc                                                                         <  1:92284/149‑1 (MQ=255)
                                                                        ggCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCACCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAAc                                                                         <  1:21569/149‑1 (MQ=255)
                                                                         gCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACg                                                                        >  1:237165/1‑149 (MQ=255)
                                                                          cATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGt                                                                       <  1:206842/149‑1 (MQ=255)
                                                                            tCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGaa                                                                          <  3:435196/144‑1 (MQ=255)
                                                                            tCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCa                                                                     >  2:299031/1‑149 (MQ=255)
                                                                            tCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGGTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGaa                                                                          >  4:435196/1‑144 (MQ=255)
                                                                              gCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAAc                                                                   >  1:212815/1‑149 (MQ=255)
                                                                               ccaccaGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACg                                                                  >  1:303767/1‑149 (MQ=255)
                                                                                  ccaGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTa                                                               <  2:303767/149‑1 (MQ=255)
                                                                                  ccaGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTACAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTa                                                               <  4:187066/149‑1 (MQ=255)
                                                                                   caGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAg                                                              >  2:489955/1‑149 (MQ=255)
                                                                                   caGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAg                                                              >  2:455754/1‑149 (MQ=255)
                                                                                          gCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTcgcg                                                                                                                            <  1:266532/80‑1 (MQ=255)
                                                                                          gCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTcgcg                                                                                                                            >  2:266532/1‑80 (MQ=255)
                                                                                          gCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTc                                                       <  1:455754/149‑1 (MQ=255)
                                                                                           ccGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTcc                                                      >  4:165632/1‑149 (MQ=255)
                                                                                            cGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTccc                                                     <  1:454157/149‑1 (MQ=255)
                                                                                             gATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCt                                                    <  4:253485/149‑1 (MQ=255)
                                                                                               ttCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGa                                                  <  2:674690/149‑1 (MQ=255)
                                                                                               ttCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGa                                                  <  2:40733/149‑1 (MQ=255)
                                                                                                tCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAg                                                 <  2:371469/149‑1 (MQ=255)
                                                                                                tCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAg                                                 <  2:237165/149‑1 (MQ=255)
                                                                                                    gCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGaa                                              <  3:226790/148‑1 (MQ=255)
                                                                                                    gCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGaa                                              >  4:226790/1‑148 (MQ=255)
                                                                                                          gTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGAATcgccgc                                       <  4:232305/149‑1 (MQ=255)
                                                                                                           tAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGt                                                                                                                        <  4:675118/67‑1 (MQ=255)
                                                                                                           tAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGt                                                                                                                        >  3:675118/1‑67 (MQ=255)
                                                                                                                  gTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGAATCGCCGCATAACCgg                               >  2:131222/1‑149 (MQ=255)
                                                                                                                  gTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGAATCGCCGCATAACCgg                               >  2:618178/1‑149 (MQ=255)
                                                                                                                      aaCCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTAGCGGCGTTCCAGAAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGAATCGCCGCATAACCGGgcgc                           >  4:507790/1‑149 (MQ=255)
                                                                                                                           gTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGAATCGCCGCATAACCGGGCGCATCCg                      >  1:488485/1‑149 (MQ=255)
                                                                                                                            tCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGAATCGCCGCATAACCGGGCGCATCCGt                     >  3:260886/1‑149 (MQ=255)
                                                                                                                            tCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGAATCGCCGCATAACCGGGCGCATCCGt                     >  2:533838/1‑149 (MQ=255)
                                                                                                                                tggtTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGAATCGCCGCATAACCGGGCGCATCCGTcccc                 <  1:131222/149‑1 (MQ=255)
                                                                                                                                       gTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGAATCGCCGCATAACCGGGCGCAt                         >  2:195573/1‑134 (MQ=255)
                                                                                                                                       gTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGAATCGCCGCATAACCGGGCGCAt                         <  1:195573/134‑1 (MQ=255)
                                                                                                                                       gTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGAATCGCCGCATAACCGGGCGCATCCGTCCCCACATTGc          <  3:165632/149‑1 (MQ=255)
                                                                                                                                        tCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGAATCGCCGCATAACCGGGCGCATCCGTCCCCACATTGcc         <  2:69389/149‑1 (MQ=255)
                                                                                                                                          gCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGAATCGCCGCATAACCGGGCGCATCCGTc                    >  2:264436/1‑136 (MQ=255)
                                                                                                                                          gCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGAATCGCCGCATAACCGGGCGCATCCGTc                    <  1:264436/136‑1 (MQ=255)
                                                                                                                                           ccATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGAATCGCCGCATAACCGGGCGCATCCGTCCCCACATTGCCGtt      <  3:526741/149‑1 (MQ=255)
                                                                                                                                               tCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGAATCGCCGCATAACCGGGCGCATCCGTCCCCACATTGCCGTTCCCa  <  4:260886/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
ATTGTAACAAAATACTCGTCACGTAGCGCTTCCAGCGGGACGGTGGAAAATGATGAGAGATGATGCTCTTCAGGCATCGCCACCAGAAACGCCGATTCATGCAACCGTAAGCTGGTAAAACCAGTCGGTGGTTCTGTCGCCATTCGCCAGATCCCGGCATCAAGTTCGCGGCGTTCCAGCAAGGCCATTTGCATCGCGGGCATCTTTTCGCGAAAAAGAACGTCAACGTTAGGATTTTCCCTGAGGAATCGCCGCATAACCGGGCGCATCCGTCCCCACATTGCCGTTCCCA  >  NC_000913/2521876‑2522167

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