Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,217,066 G→A 5.5% I114I (ATC→ATT yehX ← putative ABC transporter ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,217,0660GA 5.5% 230.1 / 2.6 73I114I (ATC→ATTyehXputative ABC transporter ATPase
Reads supporting (aligned to +/- strand):  ref base G (39/30);  new base A (2/2);  total (41/32)
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.03e-01

CAAAAGGTTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGC  >  NC_000913/2216921‑2217207
                                                                                                                                                 |                                                                                                                                             
cAAAAGGTTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCAACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATa                                                                                                                                            >  1:213683/1‑149 (MQ=255)
cAAAAGGTTCATCAATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATa                                                                                                                                            <  7:312747/149‑1 (MQ=255)
             cATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATcc                                                                                                                               >  7:261918/1‑149 (MQ=255)
              aTTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGCCAACCGGAATGCTGCTGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCTTCGATACGATCTTCAATAcg                                                                                                                              >  4:292748/1‑149 (MQ=255)
              aTTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCAATACGATCGTCAATCcg                                                                                                                              >  7:362466/1‑149 (MQ=255)
               ttAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGCCTCCCGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCcgc                                                                                                                             >  5:154378/1‑149 (MQ=255)
                  gTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGcc                                                                                                                          >  4:314373/1‑149 (MQ=255)
                   tAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGccc                                                                                                                         >  1:69488/1‑149 (MQ=255)
                   tAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGccc                                                                                                                         >  5:127741/1‑149 (MQ=255)
                    aaGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCcg                                                                                                                        >  6:5539/1‑149 (MQ=255)
                        cTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATAc                                                                                                                                           >  6:29060/1‑126 (MQ=255)
                        cTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATAc                                                                                                                                           <  5:29060/126‑1 (MQ=255)
                         ttGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGAc                                                                                                                    >  4:261604/1‑148 (MQ=255)
                               aTCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCAttttt                                                                                                             <  5:5539/149‑1 (MQ=255)
                                tCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTg                                                                                                            <  8:261918/149‑1 (MQ=255)
                                  ggCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTa                                                                                                          <  3:314373/149‑1 (MQ=255)
                                   gCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTaa                                                                                                         >  7:91671/1‑149 (MQ=255)
                                     aGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAAtt                                                                                                       <  3:261604/149‑1 (MQ=255)
                                        cAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCg                                                                                                    >  1:263058/1‑149 (MQ=255)
                                        cAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCAATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCg                                                                                                    >  7:389617/1‑149 (MQ=255)
                                                gcgcCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATtcgtcg                                                                                                                                               <  1:380243/98‑1 (MQ=255)
                                                gcgcCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATtcgtcg                                                                                                                                               >  2:380243/1‑98 (MQ=255)
                                                    cACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCcgc                                                                                                                             >  1:272333/1‑112 (MQ=255)
                                                    cACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCcgc                                                                                                                             <  2:272333/112‑1 (MQ=255)
                                                       tCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCCGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGt                                                                                     >  7:315382/1‑149 (MQ=255)
                                                         ccACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGttt                                                                                   >  3:63385/1‑149 (MQ=255)
                                                            cacGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTgc                                                                                <  6:19336/149‑1 (MQ=255)
                                                            cacGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCAATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTgc                                                                                <  8:389617/149‑1 (MQ=255)
                                                             acGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGAGTCCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTgcg                                                                               <  6:127741/149‑1 (MQ=255)
                                                                         accaccGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCGCGCTCCa                                                                     <  2:213683/147‑1 (MQ=255)
                                                                           caccGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTgg                                                                 >  8:269393/1‑149 (MQ=255)
                                                                           caccGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTgg                                                                 >  7:51505/1‑149 (MQ=255)
                                                                           caccGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTgg                                                                 >  5:169938/1‑149 (MQ=255)
                                                                            accGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTggg                                                                <  7:269393/149‑1 (MQ=255)
                                                                            accGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGAGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCAGCGACCATTTTTGTAATTGCGGCACGGTAGCTATGTTTTTCGCCCCGCTCCAGTggg                                                                >  8:69549/1‑149 (MQ=255)
                                                                              cGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTgggg                                                               >  7:272877/1‑148 (MQ=255)
                                                                                gAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGaa                                                            >  3:244834/1‑149 (MQ=255)
                                                                                  aaGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGGCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACa                                                          <  8:91671/149‑1 (MQ=255)
                                                                                  aaGCTGATGCGGAAAACGCGCACGCAAATTTGCCTCCAGCCCCACTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACa                                                          <  2:69488/149‑1 (MQ=255)
                                                                                           cGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCTCGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATa                                                  <  7:323053/148‑1 (MQ=255)
                                                                                             gATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGAt                                               >  1:327396/1‑149 (MQ=255)
                                                                                                  cGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAAt                                          <  8:272877/149‑1 (MQ=255)
                                                                                                  cGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCAATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAAt                                          <  8:362466/149‑1 (MQ=255)
                                                                                                   gCTCACGCAAATTTGCCTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATg                                         >  6:357493/1‑149 (MQ=255)
                                                                                                   gCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACa                                                          <  3:315010/132‑1 (MQ=255)
                                                                                                   gCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATg                                         <  6:169938/149‑1 (MQ=255)
                                                                                                   gCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAAATCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTAGTAATTGCGGCACGCTAGCGAAGTTTTGCGCCACGCTCCAGTCGGGGAACa                                                          >  4:315010/1‑132 (MQ=255)
                                                                                                     tCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGaa                                           <  8:1692/145‑1 (MQ=255)
                                                                                                     tCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGaa                                           >  7:1692/1‑145 (MQ=255)
                                                                                                        cGCAAATTTGACTCCAGACCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGAGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCTCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATa                                    <  7:69549/149‑1 (MQ=255)
                                                                                                         gCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGa                                                                                                                                              >  1:294420/1‑42 (MQ=255)
                                                                                                         gCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGa                                                                                                                                              <  2:294420/42‑1 (MQ=255)
                                                                                                          cAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGACAGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGTAACAGTCCAATAGATTGAATGGCATAGc                                  <  2:260594/149‑1 (MQ=255)
                                                                                                          cAAATTTGACTCAAGCCCGAGTAGCGCGATTAATTCGTCGATACGATCGTCAATCCGCGCGCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGc                                  <  2:327396/149‑1 (MQ=255)
                                                                                                              tttGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCAt                              >  6:439227/1‑149 (MQ=255)
                                                                                                                 gACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATcc                            <  4:63385/148‑1 (MQ=255)
                                                                                                                    tCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCCCGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCcggcg                        >  7:171279/1‑149 (MQ=255)
                                                                                                                    tCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCcggcg                        >  4:137316/1‑149 (MQ=255)
                                                                                                                      ccgCCCCAGTAGCGACATTAAATCGTCGCTACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGGACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCggcggc                      <  3:292748/147‑1 (MQ=255)
                                                                                                                      cAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGa                                                                                        >  1:119641/1‑83 (MQ=255)
                                                                                                                      cAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGa                                                                                        <  2:119641/83‑1 (MQ=255)
                                                                                                                       aGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGgcggcg                     >  2:233019/1‑149 (MQ=255)
                                                                                                                         ccccAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGgcgc                    >  7:300037/1‑148 (MQ=255)
                                                                                                                         ccccAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGgcgc                    >  1:221297/1‑148 (MQ=255)
                                                                                                                             aGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTc               <  2:263058/149‑1 (MQ=255)
                                                                                                                               tAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCa             <  5:357493/149‑1 (MQ=255)
                                                                                                                                  cgcCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTa          <  8:171279/149‑1 (MQ=255)
                                                                                                                                     cATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGaa                                           <  4:266034/113‑1 (MQ=255)
                                                                                                                                     cATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGaa                                           >  3:266034/1‑113 (MQ=255)
                                                                                                                                     cATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTg       >  2:12150/1‑149 (MQ=255)
                                                                                                                                     cATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATAGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTg       >  6:220542/1‑149 (MQ=255)
                                                                                                                                      aTTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATAGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTgg      <  5:220542/149‑1 (MQ=255)
                                                                                                                                          aTTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCACCTCCAGTACTGGCagc  >  2:181287/1‑149 (MQ=255)
                                                                                                                                          aTTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCagc  >  1:407621/1‑149 (MQ=255)
                                                                                                                                           ttCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCagc  >  5:40350/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                             
CAAAAGGTTCATCCATTAGTAAGACTTGCGGATCGGCAGCCAGTGCACGCGCCACTCCCACACGTTGCTGCTGACCACCGGAAAGCTGATGCGGATAACGCTCACGCAAATTTGACTCCAGCCCCAGTAGCGCCATTAATTCGTCGATACGATCGTCAATCCGCGCCCGCGACCATTTTTGTAATTGCGGCACGGTAGCGATGTTTTGCGCCACGCTCCAGTGGGGGAACAGGCCAATAGATTGAATGGCATAGCCCATCCGGCGGCGCAACTCCAGTACTGGCAGC  >  NC_000913/2216921‑2217207

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