Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,448,009 G→A G39D (GGC→GAC) 
G79G (GGC→GGT) 		 ytfP →
yzfA ←		 gamma‑glutamylamine cyclotransferase family protein YtfP
putative YjfA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,448,0090GA100.0% 222.6 / NA 68G39D (GGC→GAC) 
G79G (GGC→GGT) 		ytfP
yzfA		gamma‑glutamylamine cyclotransferase family protein YtfP
putative YjfA
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (37/31);  total (37/31)

CTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACGCCG  >  NC_000913/4447873‑4448145
                                                                                                                                        |                                                                                                                                        
ctATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCa                                                                                                                                       >  1:173870/1‑140 (MQ=255)
 tATCAGTTCGAGTTTTCACAATTCGAATATTGGTCAACGTCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCAc                                                                                                                                      <  2:464729/140‑1 (MQ=255)
  aTCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACt                                                                                                                                     >  2:72915/1‑140 (MQ=255)
      gTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATcc                                                                                                                                 <  1:828012/140‑1 (MQ=255)
      gTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATcc                                                                                                                                 <  1:565473/140‑1 (MQ=255)
       ttCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCa                                                                                                                                >  2:876921/1‑140 (MQ=255)
         cGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCCGTTGTATAGCCTGGACCACTATCCAgg                                                                                                                              >  2:1014928/1‑140 (MQ=255)
           aGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGgcg                                                                                                                            <  1:778655/140‑1 (MQ=255)
           aGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCCATAACAACCAGTTGTACAGCCTAGACCACTATCCAGgcg                                                                                                                            >  2:457070/1‑140 (MQ=255)
                  gCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTcc                                                                                                                     <  2:251613/140‑1 (MQ=255)
                      tGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCgggg                                                                                                                 >  1:939700/1‑140 (MQ=255)
                          aaTATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGggaacg                                                                                                             <  1:958176/140‑1 (MQ=255)
                             atTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGggaacggaa                                                                                                          <  1:704991/140‑1 (MQ=255)
                             atTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGggaacggaa                                                                                                          <  2:173870/140‑1 (MQ=255)
                             atTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGggaacggaa                                                                                                          <  2:58845/140‑1 (MQ=255)
                              tttGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGggaacggaac                                                                                                         <  2:171964/140‑1 (MQ=255)
                                    tACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTaca                                                                                                   <  1:242829/140‑1 (MQ=255)
                                    tACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCAATATCCAGGCGCAGTTCCGTGGAACGGAACGGTaca                                                                                                   >  2:442245/1‑140 (MQ=255)
                                     aCGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTAAAGCCTGGAACACTATCCAGGCGCAGTTCCGTGGAACGGAACGGTacac                                                                                                  >  2:322955/1‑140 (MQ=255)
                                      cGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACg                                                                                                 <  2:661690/140‑1 (MQ=255)
                                      cGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACg                                                                                                 >  1:788864/1‑140 (MQ=255)
                                       ggCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACgg                                                                                                >  2:698729/1‑140 (MQ=255)
                                         cAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTg                                                                                              <  2:759598/140‑1 (MQ=255)
                                         cAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTg                                                                                              <  2:779128/140‑1 (MQ=255)
                                                cGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTa                                                                                       <  2:108597/140‑1 (MQ=255)
                                                cGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTa                                                                                       <  2:248592/140‑1 (MQ=255)
                                                    acaaacaaGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGt                                                                                   <  2:751436/140‑1 (MQ=255)
                                                      aaacaaGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTAt                                                                                 >  1:121649/1‑140 (MQ=255)
                                                                aCAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGcc                                                                       >  1:814689/1‑140 (MQ=255)
                                                                aCAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGcc                                                                       >  1:31250/1‑140 (MQ=255)
                                                                 cAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCa                                                                      >  2:870678/1‑140 (MQ=255)
                                                                 cAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCa                                                                      >  1:770484/1‑140 (MQ=255)
                                                                  aGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCAc                                                                     >  1:556896/1‑140 (MQ=255)
                                                                       cTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTgg                                                                <  2:34480/140‑1 (MQ=255)
                                                                           aTGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGa                                                            <  2:125979/140‑1 (MQ=255)
                                                                           aTGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTGTATGGTCTTGACATCGCCACGCTGGCCGa                                                            >  1:490793/1‑140 (MQ=255)
                                                                            tGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGaa                                                           >  1:202320/1‑140 (MQ=255)
                                                                              aCCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACt                                                         >  2:95020/1‑140 (MQ=255)
                                                                              aCCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACt                                                         >  1:930581/1‑140 (MQ=255)
                                                                              aCCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACt                                                         >  1:731557/1‑140 (MQ=255)
                                                                                  aTGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGAt                                                     >  2:765214/1‑140 (MQ=255)
                                                                                  aTGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGAt                                                     >  2:150345/1‑140 (MQ=255)
                                                                                        aGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTg                                               <  1:72915/140‑1 (MQ=255)
                                                                                               gggCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCa                                        <  2:788864/140‑1 (MQ=255)
                                                                                                ggCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAg                                       >  1:595308/1‑140 (MQ=255)
                                                                                                ggCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAg                                       >  2:824536/1‑140 (MQ=255)
                                                                                                 gagATTTCAGTATCGCTAACTACCATTTGCATAGCCTCGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAgg                                      <  1:322955/138‑1 (MQ=255)
                                                                                                   gATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAgggg                                    <  1:698729/140‑1 (MQ=255)
                                                                                                   gATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAgggg                                    <  1:938653/140‑1 (MQ=255)
                                                                                                   gATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAgggg                                    <  1:525173/140‑1 (MQ=255)
                                                                                                    aTTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGc                                   >  2:923581/1‑140 (MQ=255)
                                                                                                    aTTTCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCAGGCCGAACTTGATGCCTTGCGCACCAGGGGc                                   >  2:437475/1‑140 (MQ=255)
                                                                                                    aTTTCAGTATCAATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGc                                   >  2:245440/1‑140 (MQ=255)
                                                                                                     tttCAGTATCGATAACTACCAGTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCg                                  <  2:515610/140‑1 (MQ=255)
                                                                                                     tttCAGTATCGATAAATACCAGTTGCATAGCCTGGACTACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCg                                  <  1:442245/140‑1 (MQ=255)
                                                                                                                  aaCTACCAGTTGTATAGCCTGGACCCCTACCCAGGCGGAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATAcgcgcg                     <  1:353275/140‑1 (MQ=255)
                                                                                                                          gTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGa             >  1:640603/1‑140 (MQ=255)
                                                                                                                          gTTGTATAGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCAAACTTGATACCTCGCGAACCAGGGGGGGTGAATACGCGCGCCAGTTGa             >  2:426951/1‑140 (MQ=255)
                                                                                                                               ataGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAg        >  1:542170/1‑140 (MQ=255)
                                                                                                                               ataGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAg        >  2:696054/1‑140 (MQ=255)
                                                                                                                               ataGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAg        >  2:229302/1‑140 (MQ=255)
                                                                                                                               ataGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCACACCAGGGGCGCTGCATACGCTCGCCAGTTGATTCAg        >  1:482824/1‑140 (MQ=255)
                                                                                                                                taGCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGa       <  1:718306/140‑1 (MQ=255)
                                                                                                                                taGCCTGGACCACTATCCAGGCGCAGCTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGa       >  2:707575/1‑140 (MQ=255)
                                                                                                                                  gCCTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACg     <  1:923581/140‑1 (MQ=255)
                                                                                                                                   ccTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACg     <  2:850749/139‑1 (MQ=255)
                                                                                                                                   ccTGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACg     >  1:850749/1‑139 (MQ=255)
                                                                                                                                     tGGACCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACGCCg  >  2:164134/1‑140 (MQ=255)
                                                                                                                                        |                                                                                                                                        
CTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGGGGCGGTGAATACGCGCGCCAGTTGATTCAGACGCCG  >  NC_000913/4447873‑4448145

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑