Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 322,133 G→A 22.0% intergenic (+420/+38) aguA → / ← PP_0267 agmatine deiminase/ferric siderophore receptor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947322,1330GA22.0% 109.3 / 27.4 59intergenic (+420/+38)aguA/PP_0267agmatine deiminase/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (25/21);  new base A (6/7);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 7.55e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.02e-01

AACACCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAAT  >  NC_002947/321985‑322279
                                                                                                                                                         |                                                                                                                                                  
aaCACCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCgg                                                                                                                                                    <  2:197638/149‑1 (MQ=255)
   acCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCACTACTGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCgg                                                                                                                                                    <  1:299370/149‑1 (MQ=255)
   acCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCa                                                                                                                                                 <  3:310058/149‑1 (MQ=255)
      aTGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGc                                                                                                                                              >  4:352512/1‑149 (MQ=255)
           aGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTcgcg                                                                                                                                         <  2:179748/149‑1 (MQ=255)
             cATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGaa                                                                                                                                       <  2:235830/149‑1 (MQ=255)
              aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGc                                                                                                                                              <  8:123257/141‑1 (MQ=255)
              aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGc                                                                                                                                              >  7:123257/1‑141 (MQ=255)
              aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTcgcg                                                                                                                                         >  2:209550/1‑146 (MQ=255)
               tCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCAGTGTGGCACTGATCGATTGCAGCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAaga                                                                                                                                     <  7:9361/149‑1 (MQ=255)
                  gggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGGCTCTTCGCGGGCATGCTCGCGAAGAGGc                                                                                                                                  <  1:209550/149‑1 (MQ=255)
                  gggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGc                                                                                                                                  <  3:352512/149‑1 (MQ=255)
                  gggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGa                                                                                                                                        <  5:177589/143‑1 (MQ=255)
                  gggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGa                                                                                                                                        >  6:177589/1‑143 (MQ=255)
                                 accCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACt                                                                                                                    >  3:35261/1‑148 (MQ=255)
                                   ccACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTa                                                                                                                  >  3:363119/1‑148 (MQ=255)
                                     acgacgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTaaaa                                                                                                               <  5:103727/149‑1 (MQ=255)
                                     acgacgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCcgcg                                                                                                                                         >  6:371043/1‑123 (MQ=255)
                                     acgacgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCcgcg                                                                                                                                         <  5:371043/123‑1 (MQ=255)
                                     acgacgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCcgcg                                                                                                                                         >  6:366033/1‑123 (MQ=255)
                                     acgacgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCcgcg                                                                                                                                         <  5:366033/123‑1 (MQ=255)
                                           tCTGCTTGCCGGTCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGaa                                                                                                         <  4:35261/149‑1 (MQ=255)
                                               cTTGCCGGGCGATACCCCCAATGGTTGATCA‑‑GCGAAGCAGAGATGCCTCAGAGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGaa                                                                                                         <  1:304256/147‑1 (MQ=38)
                                               cTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTg                                                                                                     >  2:43994/1‑149 (MQ=255)
                                                  gCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTa                                                                                                   <  4:363119/148‑1 (MQ=255)
                                                       gCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                     <  2:201224/109‑1 (MQ=255)
                                                       gCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                     >  1:201224/1‑109 (MQ=255)
                                                              ccccAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTc                                                                                                                                            <  1:131392/95‑1 (MQ=255)
                                                              ccccAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTc                                                                                                                                            >  2:131392/1‑95 (MQ=255)
                                                              ccccAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGaaa                                                                                                        <  6:190054/131‑1 (MQ=255)
                                                              ccccAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGaaa                                                                                                        >  5:190054/1‑131 (MQ=255)
                                                               cccAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTcacca                                                                                     >  4:32937/1‑149 (MQ=255)
                                                                     ggTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGaa                                                                                                                                       <  5:49452/93‑1 (MQ=255)
                                                                     ggTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGaa                                                                                                                                       >  6:49452/1‑93 (MQ=255)
                                                                       ttGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTa                                                                             >  2:356191/1‑149 (MQ=255)
                                                                          aTCA‑‑GCGAAGCGGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcgc                                                                          >  5:215883/1‑149 (MQ=255)
                                                                          aTCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcgc                                                                          >  6:167301/1‑149 (MQ=255)
                                                                                         gagaTGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCAc                                                               >  8:61134/1‑148 (MQ=255)
                                                                                             tGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGCCACCACAAGGCTACGCGGATCGCCCACCTGGa                                                          <  1:43994/149‑1 (MQ=255)
                                                                                                        tGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGcgccg                                                 >  7:307448/1‑149 (MQ=255)
                                                                                                                       aCTGATCGATTGCACCTGTACCTGCCTCTTCGCGGTCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGTATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGc                                  <  3:32937/149‑1 (MQ=255)
                                                                                                                       aCTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGc                                  >  3:15962/1‑149 (MQ=255)
                                                                                                                           aTCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTg                                >  2:191736/1‑147 (MQ=255)
                                                                                                                               aTTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGtagtag                          >  5:252995/1‑149 (MQ=255)
                                                                                                                                     ccTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGCGGCAGGGACAGGCAACTTAAACTCAGCAACTGTAGCTCCCCGTCCCCCCACGGCTCCGCGGAGCGCCCCCCTGGCTCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCtgg                    >  3:211722/1‑147 (MQ=37)
                                                                                                                                         tACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCa                 <  6:215883/148‑1 (MQ=255)
                                                                                                                                         tACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCa                 <  5:167301/148‑1 (MQ=255)
                                                                                                                                            cTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTgg                                           <  8:103980/119‑1 (MQ=255)
                                                                                                                                            cTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTgg                                           >  7:103980/1‑119 (MQ=255)
                                                                                                                                            cTGCCGCGTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAAGCAGAAACTGTAGCTACCCGGCCCCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGGCGCCGAGCTGTAGTAGGTCTTGTCAGCaa             <  4:211722/149‑1 (MQ=255)
                                                                                                                                             tGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAAt            >  5:189411/1‑149 (MQ=255)
                                                                                                                                                ctctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGtt         <  6:189411/149‑1 (MQ=255)
                                                                                                                                                ctctTCGCGGGCATGCCCGCGAAGAGGCAGGTACAGGCACCTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCAACCTGGATCTGAGCCGCACTGGTCGCCGAGCTGTTGTAGGGCTTGTCAGCAATGtt         >  2:53314/1‑149 (MQ=255)
                                                                                                                                                 tctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTg        >  4:91170/1‑149 (MQ=255)
                                                                                                                                                 tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGTAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCaa             >  5:92920/1‑144 (MQ=255)
                                                                                                                                                 tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGTAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCaa             <  6:92920/144‑1 (MQ=255)
                                                                                                                                                 tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTg        >  6:335461/1‑149 (MQ=255)
                                                                                                                                                 tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTg        >  2:197695/1‑149 (MQ=255)
                                                                                                                                                      gcgGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGt          <  2:200646/142‑1 (MQ=255)
                                                                                                                                                      gcgGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGt          >  1:200646/1‑142 (MQ=255)
                                                                                                                                                       cgGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAAt  <  1:356191/149‑1 (MQ=255)
                                                                                                                                                         |                                                                                                                                                  
AACACCATGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCA‑‑GCGAAGCAGAGATGC‑‑AGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAAT  >  NC_002947/321985‑322279

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