Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 581,723 C→T H64Y (CAT→TAT)  tfaD → DLP12 prophage; putative tail fiber assembly protein TfaD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913581,7230CT100.0% 231.0 / NA 65H64Y (CAT→TAT) tfaDDLP12 prophage; putative tail fiber assembly protein TfaD
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (31/34);  total (31/34)

CCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCC  >  NC_000913/581586‑581851
                                                                                                                                         |                                                                                                                                
cccGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGattat                                                                                                                                >  1:480047/1‑140 (MQ=255)
cccGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGattat                                                                                                                                >  2:938294/1‑140 (MQ=255)
      tCGGTTATTCCAAAATGCTGCTGGGTGGTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGt                                                                                                                          >  1:423505/1‑140 (MQ=255)
          ttATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTaaaa                                                                                                                      <  2:918588/140‑1 (MQ=255)
          ttATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTaaaa                                                                                                                      >  1:224744/1‑140 (MQ=255)
           tATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAAc                                                                                                                     <  2:51840/140‑1 (MQ=255)
                 aaaaTGCTGCTGGGTGGTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTa                                                                                                               <  2:423505/140‑1 (MQ=255)
                    aTGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGa                                                                                                            <  2:224744/140‑1 (MQ=255)
                       ctgctgGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGt                                                                                                         >  2:714120/1‑140 (MQ=255)
                          ctgGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATGTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGc                                                                                                      >  2:869652/1‑140 (MQ=255)
                            gggTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGt                                                                                                    <  2:287733/140‑1 (MQ=255)
                                  ttATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGt                                                                                                         <  2:746090/129‑1 (MQ=255)
                                  ttATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGt                                                                                                         >  1:746090/1‑129 (MQ=255)
                                  ttATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAgggg                                                                                              <  1:42645/140‑1 (MQ=255)
                                    aTGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGAc                                                                                            >  2:187520/1‑140 (MQ=255)
                                     tGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGAcg                                                                                           <  1:187520/140‑1 (MQ=255)
                                       ccTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGAcgcg                                                                                         >  1:528551/1‑140 (MQ=255)
                                       ccTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGAcgcg                                                                                         <  1:59531/140‑1 (MQ=255)
                                        cTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAgggg                                                                                              <  2:686904/134‑1 (MQ=255)
                                        cTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAgggg                                                                                              >  1:686904/1‑134 (MQ=255)
                                         tACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGtt                                                                                       >  1:364252/1‑140 (MQ=255)
                                          aCTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGtta                                                                                      >  1:489269/1‑140 (MQ=255)
                                           cTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGttat                                                                                     <  1:150644/140‑1 (MQ=255)
                                                  gagCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTc                                                                              <  1:151169/140‑1 (MQ=255)
                                                     cATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGa                                                                           >  1:647663/1‑140 (MQ=255)
                                                            agcGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGt                                                                    <  2:781830/140‑1 (MQ=255)
                                                             gcgcAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTc                                                                   >  2:832242/1‑140 (MQ=255)
                                                             gcgcAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTc                                                                   >  2:176487/1‑140 (MQ=255)
                                                              cgcAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTcc                                                                  >  2:196593/1‑140 (MQ=255)
                                                                 aaCACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGtt                                                               <  1:714120/140‑1 (MQ=255)
                                                                   cacCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTAc                                                             >  2:438940/1‑140 (MQ=255)
                                                                      ccTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCgg                                                          <  1:869652/140‑1 (MQ=255)
                                                                      ccTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCgg                                                          <  2:375547/140‑1 (MQ=255)
                                                                        ttATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGaa                                                        >  2:712893/1‑140 (MQ=255)
                                                                        ttATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGaa                                                        >  2:707465/1‑140 (MQ=255)
                                                                         tATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGaaa                                                       <  1:438940/140‑1 (MQ=255)
                                                                           ttCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAAt                                                     <  2:859870/140‑1 (MQ=255)
                                                                                gCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTAc                                                >  2:924965/1‑140 (MQ=255)
                                                                                gCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTAc                                                >  2:99846/1‑140 (MQ=255)
                                                                                gCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTAc                                                <  1:99846/140‑1 (MQ=255)
                                                                                gCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTAc                                                >  1:648930/1‑140 (MQ=255)
                                                                                gCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCAGTCCGTTACCGGAAAATGTTAc                                                >  1:696023/1‑140 (MQ=255)
                                                                                 cTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTAcc                                               >  1:40537/1‑140 (MQ=255)
                                                                                    gCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTgg                                            >  2:554778/1‑140 (MQ=255)
                                                                                         gTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATc                                       <  2:40537/140‑1 (MQ=255)
                                                                                              tGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCgg                                  <  2:417663/140‑1 (MQ=255)
                                                                                                  ttCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAgg                              <  1:938294/140‑1 (MQ=255)
                                                                                                  ttCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAgg                              <  1:691633/140‑1 (MQ=255)
                                                                                                  ttCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAgg                              <  1:695281/140‑1 (MQ=255)
                                                                                                    cAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAgggg                            <  1:554778/140‑1 (MQ=255)
                                                                                                    cAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAgggg                            <  2:648930/140‑1 (MQ=255)
                                                                                                          tgatgaGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTc                      >  1:377201/1‑140 (MQ=255)
                                                                                                            atgaGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAg                    >  1:188718/1‑140 (MQ=255)
                                                                                                            atgaGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAg                    >  1:67078/1‑140 (MQ=255)
                                                                                                               aGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAg                 <  1:176487/140‑1 (MQ=255)
                                                                                                               aGTCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAg                 <  2:67078/140‑1 (MQ=255)
                                                                                                                 tCATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTg               >  1:31552/1‑140 (MQ=255)
                                                                                                                  cATCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTgg              <  2:31552/140‑1 (MQ=255)
                                                                                                                   aTCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGa             <  2:364252/140‑1 (MQ=255)
                                                                                                                    tCGTGGCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGaa            >  1:695834/1‑140 (MQ=255)
                                                                                                                        ggCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGt        <  1:196593/140‑1 (MQ=255)
                                                                                                                         gCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTa       <  2:489269/140‑1 (MQ=255)
                                                                                                                         gCATCTCGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTa       >  1:201917/1‑140 (MQ=255)
                                                                                                                            tctcGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAg    >  2:137424/1‑140 (MQ=255)
                                                                                                                              tcGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGcc  <  1:712893/140‑1 (MQ=255)
                                                                                                                              tcGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGcc  <  2:201917/140‑1 (MQ=255)
                                                                                                                               cGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTa       <  2:418345/134‑1 (MQ=255)
                                                                                                                               cGTTGAAGATTATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTa       >  1:418345/1‑134 (MQ=255)
                                                                                                                                         |                                                                                                                                
CCCGTGTCGGTTATTCCAAAATGCTGCTGGGTGTTTATGCCTACTTTATAGAGCATAAGCAGCGCAACACCCTTATTCCAGCTGGCTTCGTGGCTGTTTTCAACAGTGATGAGTCATCGTGGCATCTCGTTGAAGATCATCGGGGTAAAACGGTTTATGACGTAGCGTCAGGGGACGCGTTATTTATTTCTGAACTCGGTCCGTTACCGGAAAATGTTACCTGGTTATCGCCGGAAGGGGAGTTTCAGAAGTGGAACGGTACAGCC  >  NC_000913/581586‑581851

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