Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,904,586 T→C V231A (GTC→GCC)  manZ → mannose‑specific enzyme IID component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,904,5860TC100.0% 293.8 / NA 77V231A (GTC→GCC) manZmannose‑specific enzyme IID component of PTS
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (42/35);  total (42/35)

CAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCG  >  NC_000913/1904450‑1904714
                                                                                                                                        |                                                                                                                                
cAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGccc                                                                                                                                <  2:175751/139‑1 (MQ=255)
   aaaCTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGa                                                                                                                             <  2:3021/139‑1 (MQ=255)
        gACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTAtt                                                                                                                        <  2:1039016/139‑1 (MQ=255)
         aCGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTc                                                                                                                       >  1:636018/1‑139 (MQ=255)
              aGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGAc                                                                                                                  <  1:246124/139‑1 (MQ=255)
                 ggCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAg                                                                                                               >  1:1147259/1‑139 (MQ=255)
                         tCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGcc                                                                                                       <  1:1027297/139‑1 (MQ=255)
                         tCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGcc                                                                                                       >  2:198948/1‑139 (MQ=255)
                          ccTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCa                                                                                                      <  1:1211602/139‑1 (MQ=255)
                          ccTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCa                                                                                                      <  1:198948/139‑1 (MQ=255)
                           cTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAg                                                                                                     >  2:851331/1‑139 (MQ=255)
                                  tGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGt                                                                                              <  2:1292957/139‑1 (MQ=255)
                                   gTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTa                                                                                             >  1:786719/1‑139 (MQ=255)
                                   gTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTa                                                                                             >  1:353977/1‑139 (MQ=255)
                                              gggCATTTGTTAACAAGTGGACACGTGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGATCACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTCGTACCActgctgct                                                                                  >  1:922022/1‑139 (MQ=255)
                                              gggCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACATGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctgct                                                                                  >  2:914639/1‑139 (MQ=255)
                                              gggCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGcc                                                                                                  >  1:219415/1‑123 (MQ=255)
                                              gggCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGcc                                                                                                  <  2:219415/123‑1 (MQ=255)
                                              gggCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctgct                                                                                  >  1:1387357/1‑139 (MQ=255)
                                              gggCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctgct                                                                                  >  1:1324663/1‑139 (MQ=255)
                                                 cATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAc                                                                               <  1:1257904/139‑1 (MQ=255)
                                                 cATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAc                                                                               <  1:795953/139‑1 (MQ=255)
                                                 cATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAc                                                                               <  2:800055/139‑1 (MQ=255)
                                                  aTTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGAcc                                                                              >  2:1283790/1‑139 (MQ=255)
                                                   ttGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCt                                                                             <  1:298636/139‑1 (MQ=255)
                                                       ttAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGc                                                                         <  2:1387357/139‑1 (MQ=255)
                                                          aCAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTg                                                                      <  2:636018/139‑1 (MQ=255)
                                                              gTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTAtg                                                                  >  2:419386/1‑139 (MQ=255)
                                                               tGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTAtgt                                                                 <  2:1147259/139‑1 (MQ=255)
                                                               tGGACACATGTCAACATCCCGCTGGTTGTCTCTCCCATTACTGACCAGACGTGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTAtgt                                                                 >  2:635100/1‑139 (MQ=255)
                                                                ggACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTAtgtg                                                                >  1:701096/1‑139 (MQ=255)
                                                                 gACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGAc                                                                                                                  >  2:420664/1‑88 (MQ=255)
                                                                 gACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGAc                                                                                                                  <  1:420664/88‑1 (MQ=255)
                                                                 gACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTgg                                                               >  1:1369917/1‑139 (MQ=255)
                                                                        gTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTgc                                                        >  1:1426327/1‑139 (MQ=255)
                                                                        gTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTgc                                                        >  2:247192/1‑139 (MQ=255)
                                                                         tCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTgcg                                                       >  1:435551/1‑139 (MQ=255)
                                                                              aTCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGc                                                                                                   <  1:762837/90‑1 (MQ=255)
                                                                              aTCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGc                                                                                                   >  2:762837/1‑90 (MQ=255)
                                                                                cccGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCa                                                                                                      >  1:719263/1‑85 (MQ=255)
                                                                                cccGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCa                                                                                                      <  2:719263/85‑1 (MQ=255)
                                                                                  cGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCGCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGt                                              >  2:633153/1‑139 (MQ=255)
                                                                                  cGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGGCTATTCTGGAc                                                                                                                  <  2:316766/71‑1 (MQ=255)
                                                                                  cGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGAc                                                                                                                  >  1:316766/1‑71 (MQ=255)
                                                                                       gTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAcc                                         <  2:1324663/139‑1 (MQ=255)
                                                                                       gTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAcc                                         >  2:654325/1‑139 (MQ=255)
                                                                                        ttGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGcc                                                                                                  <  1:1212727/81‑1 (MQ=255)
                                                                                        ttGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGcc                                                                                                  >  2:1212727/1‑81 (MQ=255)
                                                                                          gTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGt                                                                                                              <  2:970480/67‑1 (MQ=255)
                                                                                          gTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGt                                                                                                              >  1:970480/1‑67 (MQ=255)
                                                                                          gTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTg                                                                      <  2:898702/107‑1 (MQ=255)
                                                                                          gTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTg                                                                      >  1:898702/1‑107 (MQ=255)
                                                                                          gTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGc                                      <  1:419386/139‑1 (MQ=255)
                                                                                           tctctcGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCt                                     <  2:1369917/139‑1 (MQ=255)
                                                                                                 gCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGAc                                                                                                                  >  1:801053/1‑56 (MQ=255)
                                                                                                 gCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGAc                                                                                                                  <  2:801053/56‑1 (MQ=255)
                                                                                                         gACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATgac                                                                                                       >  2:1468342/1‑57 (MQ=255)
                                                                                                         gACCAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGcc                                                                                                       <  1:1468342/59‑1 (MQ=255)
                                                                                                           ccAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCaa                                                    <  1:449105/108‑1 (MQ=255)
                                                                                                           ccAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCaa                                                    >  2:449105/1‑108 (MQ=255)
                                                                                                            cAGACGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGct                    >  1:17327/1‑139 (MQ=255)
                                                                                                               aCGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGAc                                                                                                                  <  2:356591/42‑1 (MQ=255)
                                                                                                               aCGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGAc                                                                                                                  >  1:356591/1‑42 (MQ=255)
                                                                                                               aCGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAc                                          >  2:74234/1‑114 (MQ=255)
                                                                                                               aCGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAc                                          <  1:74234/114‑1 (MQ=255)
                                                                                                               aCGGGCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGcttct                 >  1:329195/1‑139 (MQ=255)
                                                                                                                 gggCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGcc                                                                                                       <  1:1009794/51‑1 (MQ=255)
                                                                                                                 gggCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGcc                                                                                                       >  2:1009794/1‑51 (MQ=255)
                                                                                                                 gggCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCttcttc               >  1:889765/1‑139 (MQ=255)
                                                                                                                 gggCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCttcttc               >  1:1354129/1‑139 (MQ=255)
                                                                                                                 gggCAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCttcttc               >  1:1190813/1‑139 (MQ=255)
                                                                                                                    cAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAcc                                         <  1:1303639/110‑1 (MQ=255)
                                                                                                                    cAAAGAACACGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAAcc                                         >  2:1303639/1‑110 (MQ=255)
                                                                                                                             cGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATc   <  1:654325/139‑1 (MQ=255)
                                                                                                                             cGTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATc   >  1:131011/1‑139 (MQ=255)
                                                                                                                              gTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctgctg                                                                                 >  1:1457766/1‑60 (MQ=255)
                                                                                                                              gTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCActgctgctg                                                                                 <  2:1457766/60‑1 (MQ=255)
                                                                                                                              gTTACTACTGCCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCg  <  1:968285/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                
CAAAAACTGACGGAAGGGGCGTCTATCCTCGGCCTGTTTGTCATGGGGGCATTGGTTAACAAGTGGACACATGTCAACATCCCGCTGGTTGTCTCTCGCATTACTGACCAGACGGGCAAAGAACACGTTACTACTGTCCAGACTATTCTGGACCAGTTAATGCCAGGCCTGGTACCACTGCTGCTGACCTTTGCTTGTATGTGGCTACTGCGCAAAAAAGTTAACCCGCTGTGGATCATCGTTGGCTTCTTCGTCATCGGTATCG  >  NC_000913/1904450‑1904714

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