Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,137,229 A→G S143S (AGT→AGC)  yijF ← DUF1287 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,137,2290AG100.0% 322.8 / NA 87S143S (AGT→AGC) yijFDUF1287 family protein
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base G (46/41);  total (46/41)

CTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACTCAGCAAAATTCTT  >  NC_000913/4137096‑4137363
                                                                                                                                     |                                                                                                                                      
cttcCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTggg                                                                                                                                   >  2:1680311/1‑139 (MQ=255)
  tcCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGt                                                                                                                                 >  2:1162833/1‑139 (MQ=255)
     tGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCt                                                                                                                              <  2:1307651/139‑1 (MQ=255)
       cgcgcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTg                                                                                                                            >  2:992816/1‑139 (MQ=255)
       cgcgcCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTg                                                                                                                            >  2:104663/1‑139 (MQ=255)
                gcCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTgggg                                                                                                                   <  1:1680311/139‑1 (MQ=255)
                 ccAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                  <  1:104663/139‑1 (MQ=255)
                 ccAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                  <  1:14285/139‑1 (MQ=255)
                    aTGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTgggg                                                                                                                   >  1:1469497/1‑135 (MQ=255)
                    aTGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTgggg                                                                                                                   <  2:1469497/135‑1 (MQ=255)
                         aTGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCg                                                                                                                 >  2:1581516/1‑132 (MQ=255)
                         aTGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCg                                                                                                                 <  1:1581516/132‑1 (MQ=255)
                              tCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATg                                                                                                     <  2:62581/139‑1 (MQ=255)
                              tCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATg                                                                                                     <  1:992816/139‑1 (MQ=255)
                               cACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATg                                                                                                     >  2:485071/1‑138 (MQ=255)
                               cACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATg                                                                                                     <  1:485071/138‑1 (MQ=255)
                                aCCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTc                                                                                                                               >  2:907973/1‑111 (MQ=255)
                                aCCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTc                                                                                                                               <  1:907973/111‑1 (MQ=255)
                                 ccAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCt                                                                                                                                      >  2:797975/1‑103 (MQ=255)
                                 ccAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTg                                                                                                                                     <  2:307078/104‑1 (MQ=255)
                                 ccAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTg                                                                                                                                     >  1:307078/1‑104 (MQ=255)
                                 ccAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAAAATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATggcg                                                                                                  >  2:997030/1‑139 (MQ=255)
                                    gcggcgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCtggt                                                                                                                       <  2:1491765/115‑1 (MQ=255)
                                    gcggcgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCtggt                                                                                                                       >  1:1491765/1‑115 (MQ=255)
                                    gcggcgTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCt                                                                                               >  2:827735/1‑139 (MQ=255)
                                          tACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCa                                                                                         >  1:504695/1‑139 (MQ=255)
                                               tCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGttt                                                                                    >  2:112565/1‑139 (MQ=255)
                                               tCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGACAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGttt                                                                                    >  1:1571235/1‑139 (MQ=255)
                                                cgcgGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTc                                                                                   >  1:1461977/1‑139 (MQ=255)
                                                  cgGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCa                                                                                 <  1:997030/139‑1 (MQ=255)
                                                    ggCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                  >  2:1442826/1‑104 (MQ=255)
                                                    ggCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                  <  1:1442826/104‑1 (MQ=255)
                                                             aTCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCAcc                                                                      <  2:379860/139‑1 (MQ=255)
                                                              tCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATg                                                                                                     <  1:436007/107‑1 (MQ=255)
                                                              tCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATg                                                                                                     >  2:436007/1‑107 (MQ=255)
                                                               cTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                  <  1:295869/93‑1 (MQ=255)
                                                               cTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                  >  2:295869/1‑93 (MQ=255)
                                                               cTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAc                                                                                           >  2:1292378/1‑116 (MQ=255)
                                                               cTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAAc                                                                                           <  1:1292378/116‑1 (MQ=255)
                                                               cTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCg                                                                    <  1:827735/139‑1 (MQ=255)
                                                                tGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGa                                                                   <  2:1461977/139‑1 (MQ=255)
                                                                  aaaCCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACg                                                                 <  2:733123/139‑1 (MQ=255)
                                                                        ccccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                  >  2:1309081/1‑84 (MQ=255)
                                                                        ccccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                  <  1:1309081/84‑1 (MQ=255)
                                                                         cccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                  <  1:746626/83‑1 (MQ=255)
                                                                         cccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                  >  2:746626/1‑83 (MQ=255)
                                                                         cccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCg                                                          >  1:770716/1‑139 (MQ=255)
                                                                         cccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCg                                                          >  1:411690/1‑139 (MQ=255)
                                                                          ccAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGa                                                         <  2:770716/139‑1 (MQ=255)
                                                                            aaTATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                  <  1:1002613/80‑1 (MQ=255)
                                                                            aaTATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGc                                                                                                                  >  2:1002613/1‑80 (MQ=255)
                                                                                tgtgCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGc                                                   >  2:97085/1‑139 (MQ=255)
                                                                                  tgCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTg                                                 >  2:1050507/1‑139 (MQ=255)
                                                                                   gCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTgggg                                                                                                                                  <  1:1113311/57‑1 (MQ=255)
                                                                                   gCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTgggg                                                                                                                                  >  2:1113311/1‑57 (MQ=255)
                                                                                   gCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGt                                                >  2:369840/1‑139 (MQ=255)
                                                                                    ccAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTc                                               <  2:1571235/139‑1 (MQ=255)
                                                                                      aGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCg                                             >  2:694760/1‑139 (MQ=255)
                                                                                      aGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCg                                             <  1:175943/139‑1 (MQ=255)
                                                                                       gCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCtggt                                                                                                                       <  2:539802/64‑1 (MQ=255)
                                                                                       gCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCtggt                                                                                                                       >  1:539802/1‑64 (MQ=255)
                                                                                        cccGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCa                                                                                         <  1:114453/93‑1 (MQ=255)
                                                                                        cccGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCa                                                                                         >  2:114453/1‑93 (MQ=255)
                                                                                            ttATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGtt                                                                                                                                <  2:790272/50‑1 (MQ=255)
                                                                                            ttATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGtt                                                                                                                                >  1:790272/1‑50 (MQ=255)
                                                                                            ttATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCt                                       >  2:903940/1‑139 (MQ=255)
                                                                                              aTCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCttt                                     <  1:97085/139‑1 (MQ=255)
                                                                                                  aGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTgggg                                                                                                                                  >  2:1580377/1‑42 (MQ=255)
                                                                                                  aGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTgggg                                                                                                                                  <  1:1580377/42‑1 (MQ=255)
                                                                                                   gTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCgt                                                                                                                >  2:133111/1‑59 (MQ=255)
                                                                                                   gTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCgt                                                                                                                <  1:133111/59‑1 (MQ=255)
                                                                                                   gTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAg                                                                           >  2:709546/1‑96 (MQ=255)
                                                                                                   gTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAg                                                                           <  1:751921/96‑1 (MQ=255)
                                                                                                   gTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAg                                                                           <  1:709546/96‑1 (MQ=255)
                                                                                                   gTCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAg                                                                           >  2:751921/1‑96 (MQ=255)
                                                                                                    tCGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTc                               >  1:1658421/1‑139 (MQ=255)
                                                                                                     cGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCtggt                                                                                                                       <  1:1500661/50‑1 (MQ=255)
                                                                                                     cGCCAGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCtggt                                                                                                                       >  2:1500661/1‑50 (MQ=255)
                                                                                                         aGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTggg                                                                                                                    <  1:1710325/49‑1 (MQ=255)
                                                                                                         aGGAGACAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTggg                                                                                                                    >  2:1710325/1‑49 (MQ=255)
                                                                                                           gagaCAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCtggtg                                                                                                                      <  1:1463436/45‑1 (MQ=255)
                                                                                                           gagaCAATATCGCCCGCTTGATAGTCGCTGGGGTTCTTGCtggtg                                                                                                                      >  2:1463436/1‑45 (MQ=255)
                                                                                                                      gCCCGCTTGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACTCa             >  1:444065/1‑139 (MQ=255)
                                                                                                                            ttGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCAtt                           <  2:835625/119‑1 (MQ=255)
                                                                                                                            ttGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCAtt                           >  1:835625/1‑119 (MQ=255)
                                                                                                                            ttGATAGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACTCAGCaaaa       >  2:473889/1‑139 (MQ=255)
                                                                                                                                 aGTCGCTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACTCAGCAAAAttctt  <  1:473889/139‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                      
CTTCCTGCGCGCCTGCGCCAATGTTATGTATCACCAGCGGCGTACCGTCGCGGGCGAAGCCATCTGAAACCACCCCAATATGTGCCAGCCCGTTATCCAGTCGCCAGGAGACAATATCGCCCGCTTGATAGTCACTGGGGTTCTTGCTGGTGGGGCGTGTTTTATCATGGCGGCTAAACCAGGTTTCCAGATTAGGCACCCGACGGTGATCGATGTTGCTGTCCGGGCGCTTTAACTTCCATTTTTGCGGGTACTCAGCAAAATTCTT  >  NC_000913/4137096‑4137363

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