Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,472,033 G→A I494I (ATC→ATT fusA ← protein chain elongation factor EF‑G, GTP‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,472,0330GA98.7% 241.1 / ‑4.3 76I494I (ATC→ATTfusAprotein chain elongation factor EF‑G, GTP‑binding
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (39/36);  total (40/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.14e-01

CGTTGATGAACTCGTAGCCTTTCGGGTTTGAACCCGGCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGGATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTCTGGTTAGA  >  NC_000913/3471900‑3472158
                                                                                                                                     |                                                                                                                             
cGTTGATGAACTCGTAGCCTTTCGGGTTTGAACCCGGCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGGATAGt                                                                                                                          >  2:109972/1‑139 (MQ=255)
 gTTGATGAACTCGTAGCCTTTCGGGTTTGAACCCGGCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGtt                                                                                                                         >  1:586846/1‑139 (MQ=255)
 gTTGATGAACTCGTAGCCTTTCGGGTTTGAACCCGGCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGtt                                                                                                                         >  1:760356/1‑139 (MQ=255)
            cGTAGCCTTTCGGGTTTGAACCCGGCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGc                                                                                                              >  2:392818/1‑139 (MQ=255)
            cGTAGCCTTTCGGGTTTGAACCCGGCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGc                                                                                                              >  2:393334/1‑139 (MQ=255)
                  cTTTCGGGTTTGAACCCGGCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTg                                                                                                        <  1:536625/139‑1 (MQ=255)
                    ttCGGGTTTGAACCCGGCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCg                                                                                                      >  2:1217643/1‑139 (MQ=255)
                    ttCGGGTTTGAACCCGGCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCg                                                                                                      <  1:393334/139‑1 (MQ=255)
                    ttCGGGTTTGAACCCGGCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCg                                                                                                      >  1:133574/1‑139 (MQ=255)
                    ttCGGGTTTGAACCCGGCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCg                                                                                                      <  2:867533/139‑1 (MQ=255)
                    ttCGGGTTTGAACCCGGCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATAGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCg                                                                                                      <  1:448860/139‑1 (MQ=255)
                           ttGAACCCGGCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTAcc                                                                                               >  2:626248/1‑139 (MQ=255)
                           ttGAACCCGGCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTAcc                                                                                               >  2:651203/1‑139 (MQ=255)
                            tGAACCCGGCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCt                                                                                              <  2:570577/139‑1 (MQ=255)
                                 ccGGCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGtt                                                                                         <  1:392818/139‑1 (MQ=255)
                                   ggCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCAc                                                                                                                     >  2:329023/1‑109 (MQ=255)
                                   ggCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCAc                                                                                                                     <  1:329023/109‑1 (MQ=255)
                                   ggCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGtt                                                                                         <  1:1047806/137‑1 (MQ=255)
                                   ggCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGtt                                                                                         >  2:1047806/1‑137 (MQ=255)
                                   ggCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCg                                                                                       >  2:1086737/1‑139 (MQ=255)
                                                aCATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCAc                                                                                                                     <  2:435705/96‑1 (MQ=255)
                                                aCATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCAc                                                                                                                     >  1:435705/1‑96 (MQ=255)
                                                   tGTCGATAACAACATGACCATCCTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTc                                                                       <  1:651203/139‑1 (MQ=255)
                                                   tGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTc                                                                       <  1:1086737/139‑1 (MQ=255)
                                                       gATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGc                                                                   >  2:975729/1‑139 (MQ=255)
                                                         tAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTg                                                                                                        <  2:788340/100‑1 (MQ=255)
                                                         tAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTg                                                                                                        >  1:788340/1‑100 (MQ=255)
                                                            caacaTGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCAt                                                              >  2:1201151/1‑139 (MQ=255)
                                                             aacaTGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGtt                                                                                         <  2:1184611/111‑1 (MQ=255)
                                                             aacaTGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGtt                                                                                         >  1:1184611/1‑111 (MQ=255)
                                                             aacaTGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATa                                                             <  1:975729/139‑1 (MQ=255)
                                                              acaTGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATAc                                                            <  1:1217643/139‑1 (MQ=255)
                                                               caTGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTg                                                                                                        >  1:653606/1‑94 (MQ=255)
                                                               caTGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTg                                                                                                        <  2:653606/94‑1 (MQ=255)
                                                                       tACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAAcc                                                                                                          <  2:372341/84‑1 (MQ=255)
                                                                       tACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAAcc                                                                                                          >  1:372341/1‑84 (MQ=255)
                                                                       tACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACga                                                   >  2:962574/1‑139 (MQ=255)
                                                                       tACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACga                                                   >  2:381032/1‑139 (MQ=255)
                                                                            accacgaccacCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCt                                                                                     >  2:486699/1‑100 (MQ=255)
                                                                            accacgaccacCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCt                                                                                     <  1:486699/100‑1 (MQ=255)
                                                                              cacgaccacCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTc                                            >  2:839305/1‑139 (MQ=255)
                                                                               acgaccacCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCg                                           >  1:422484/1‑139 (MQ=255)
                                                                                 gaccacCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAAcc                                                                                                          >  1:221499/1‑74 (MQ=255)
                                                                                 gaccacCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAAcc                                                                                                          <  2:221499/74‑1 (MQ=255)
                                                                                      ccaGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAg                                                                                                               >  1:1131755/1‑64 (MQ=255)
                                                                                      ccaGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAg                                                                                                               <  2:1131755/64‑1 (MQ=255)
                                                                                        aGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGt                                  <  1:345133/139‑1 (MQ=255)
                                                                                              tttCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGt                                  >  2:1211702/1‑133 (MQ=255)
                                                                                              tttCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGt                                  <  1:1211702/133‑1 (MQ=255)
                                                                                              tttCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGccc                            >  1:491045/1‑139 (MQ=255)
                                                                                              tttCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGccc                            >  1:691267/1‑139 (MQ=255)
                                                                                              tttCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGccc                            >  1:925381/1‑139 (MQ=255)
                                                                                                 cgcgTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATa                         >  2:1169075/1‑139 (MQ=255)
                                                                                                 cgcgTGTTTACCTTCAACACGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATAc                        >  2:831286/1‑139 (MQ=255)
                                                                                                   cgTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATAcc                       >  2:295145/1‑139 (MQ=255)
                                                                                                   cgTGTTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATAcc                       >  1:699634/1‑139 (MQ=255)
                                                                                                     tgtTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCcg                     <  1:269521/139‑1 (MQ=255)
                                                                                                      gtTTACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCcgc                    >  1:99209/1‑139 (MQ=255)
                                                                                                        ttACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGtt                                                                             <  1:153583/80‑1 (MQ=255)
                                                                                                        ttACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGtt                                                                             >  2:153583/1‑80 (MQ=255)
                                                                                                         tACCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCgat                 >  2:1200145/1‑139 (MQ=255)
                                                                                                          aCCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCgatg                <  1:962574/139‑1 (MQ=255)
                                                                                                          aCCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCgatg                <  2:133574/139‑1 (MQ=255)
                                                                                                          aCCTTCAACATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCgatg                <  1:1201151/139‑1 (MQ=255)
                                                                                                             ttCAACACGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTCGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATgg            <  1:831286/139‑1 (MQ=255)
                                                                                                                aaCATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTc          <  1:1200145/139‑1 (MQ=255)
                                                                                                                aaCATCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTc          <  2:699634/139‑1 (MQ=255)
                                                                                                                    tCGGTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTCTGGt      <  1:295145/139‑1 (MQ=255)
                                                                                                                      ggTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACgatga                                                >  2:663922/1‑95 (MQ=255)
                                                                                                                      ggTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACgatga                                                <  1:663922/95‑1 (MQ=255)
                                                                                                                      ggTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTCTGGTTa    <  2:422484/139‑1 (MQ=255)
                                                                                                                       gTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTCTgg       <  1:378542/135‑1 (MQ=255)
                                                                                                                       gTAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTCTgg       >  2:378542/1‑135 (MQ=255)
                                                                                                                        tAACTTTCTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTCTGGTTAGa  <  2:691267/139‑1 (MQ=255)
                                                                                                                               cTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATgg            <  1:406160/122‑1 (MQ=255)
                                                                                                                               cTGGCGAATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATgg            >  2:406160/1‑122 (MQ=255)
                                                                                                                                     |                                                                                                                             
CGTTGATGAACTCGTAGCCTTTCGGGTTTGAACCCGGCTCCAGCGGGTACATGTCGATAACAACATGACCATACTGACCACGACCACCAGACTGTTTCGCGTGTTTACCTTCAACATCGGTAACTTTCTGGCGGATAGTTTCACGGTAAGCAACCTGCGGTTTACCTACGTTCGCTTCAACGTTGAATTCACGCTTCATACGGTCAACGATGATGTCGAGGTGCAGTTCGCCCATACCCGCGATGATGGTCTGGTTAGA  >  NC_000913/3471900‑3472158

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: