Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,767,500 T→C intergenic (+145/‑210) rnlB → / → yfjP CP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,767,5000TC100.0% 236.8 / NA 69intergenic (+145/‑210)rnlB/yfjPCP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, 50S ribosome‑binding GTPase family protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (31/38);  total (31/38)

TTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCA  >  NC_000913/2767365‑2767634
                                                                                                                                        |                                                                                                                                      
ttctAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTcccc                                                                                                                                     <  1:997990/139‑1 (MQ=255)
        gAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTa                                                                                                                             <  2:36394/139‑1 (MQ=255)
        agacaggTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTa                                                                                                                             <  1:201326/133‑1 (MQ=255)
           ccAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAAc                                                                                                                          <  2:514236/139‑1 (MQ=255)
           ccAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCACTCACAACGTCCTGCATAGGGCACACGATTTCCACTCTGAACCCGCGACTAAATTCCCCCCCCTTTaacc                                                                                                                          >  2:341223/1‑137 (MQ=255)
            cAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACt                                                                                                                         <  1:571401/139‑1 (MQ=255)
                      ttCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACTCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCTTCCCATTAACATGCTCaaaaat                                                                                                               >  2:109795/1‑139 (MQ=255)
                      ttCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAAACCCCCCTCACAACGTCCTGCATAGCGAACACGTTTTCCCCTGTGCACCCACGACTAACTTCCCCCCCCTTTAAACTTCTCaataat                                                                                                               >  1:789941/1‑139 (MQ=255)
                        cGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCACCTGTGCACCCACGATTACATTCCCTCCCCTTTAAACTACTCCATAATGt                                                                                                             >  2:529433/1‑139 (MQ=255)
                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCGCCTCACAACGTCATGCATAGAAAACCAGTTTTACCCTGTGCGCCCCCGAATAAATTCACCCCCCTTTAAACTTCCCaatctgtct                                                                                                            >  2:201326/1‑133 (MQ=255)
                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACACCGCCCTACATATGGAAGTAGTTTTACCCTTAGAACCAACTACCAAATACCTCCCCCTTTTAACTTCTCAATAATGTc                                                                                                            >  1:36864/1‑139 (MQ=255)
                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGTGAACACGGTTTCCCCTGTACACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTc                                                                                                            >  2:1281238/1‑139 (MQ=255)
                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTcccc                                                                                                                                     <  1:1081469/114‑1 (MQ=255)
                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGATTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTc                                                                                                            >  1:719453/1‑139 (MQ=255)
                         gCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGACCTGCATAGGGCACACGTTTTCCCCTGTGCACCCACGACTAAATTcccc                                                                                                                                     >  2:1081469/1‑114 (MQ=255)
                                      gctatTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCATTAATGTCACCAATTTTCTgg                                                                                               <  2:886524/137‑1 (MQ=255)
                                                  tCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTatgat                                                                                   <  1:385393/139‑1 (MQ=255)
                                                  tCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCACAATAATGCAACCAATTTTCTGGAAATCTTatgct                                                                                   >  2:894052/1‑137 (MQ=255)
                                                     aaTATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGcc                                                                                <  1:529433/139‑1 (MQ=255)
                                                         tGATCAAATAGCATTCACCCCCCCTCGCAACCTTCTGTTTAGCGTAGACCTTCTCCCCAGTGCCCCCACGACGAAATTCCCCCCCCTTAACACTTCTCAATAATGTCGCTAATTTTCTGGAAATCCTATTTTGACTTTc                                                                            <  2:1491277/139‑1 (MQ=255)
                                                          gATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTATACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTcc                                                                           <  1:274758/139‑1 (MQ=255)
                                                                ttAGCATTAACCCCCCCTCACAACGTCCTGCATAGAGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCTGTaa                                                                     <  1:1157512/138‑1 (MQ=255)
                                                                    cATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTTCCCTGTGCACCCCCGACTAAATTCCCCCCACTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATgg                                                                 <  2:36864/139‑1 (MQ=255)
                                                                      ttAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAAt                                                                                           >  2:935468/1‑111 (MQ=255)
                                                                      ttAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAAt                                                                                           <  1:935468/111‑1 (MQ=255)
                                                                         aCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTc                                                                                                  <  2:801640/101‑1 (MQ=255)
                                                                         aCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtg                                                            >  2:1342719/1‑139 (MQ=255)
                                                                         aCCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGACACCAAtttta                                                                                                  >  1:801640/1‑100 (MQ=255)
                                                                          cccccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTTTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCtgt                                                           >  1:1284369/1‑139 (MQ=255)
                                                                             ccccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAg                                                        <  1:908247/139‑1 (MQ=255)
                                                                              cccTCACAACGTCCTTCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTgg                                                                                               <  1:801619/99‑1 (MQ=255)
                                                                              cccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTaaa                                                                                                                           >  1:1333856/1‑71 (MQ=255)
                                                                              cccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTaaa                                                                                                                           <  2:1333856/71‑1 (MQ=255)
                                                                              cccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTgg                                                                                               >  2:801619/1‑99 (MQ=255)
                                                                              cccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGt                                                                       <  2:44971/123‑1 (MQ=255)
                                                                              cccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGt                                                       >  1:513206/1‑139 (MQ=255)
                                                                              cccTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCCCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGt                                                                       >  1:44971/1‑123 (MQ=255)
                                                                                  cacaACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTaaa                                                                                                                           >  1:507711/1‑67 (MQ=255)
                                                                                  cacaACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTaaa                                                                                                                           <  2:507711/67‑1 (MQ=255)
                                                                                       cGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGa                                              >  2:1487050/1‑139 (MQ=255)
                                                                                       cGTCCTGCATAGGGAACACGTGTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGa                                              <  1:341223/139‑1 (MQ=255)
                                                                                         tCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGAtt                                            >  1:579049/1‑139 (MQ=255)
                                                                                          ccTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGa                                              <  1:1164786/136‑1 (MQ=255)
                                                                                          ccTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGa                                              >  2:1164786/1‑136 (MQ=255)
                                                                                            tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTg                                         <  1:1075620/139‑1 (MQ=255)
                                                                                            tGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTg                                         <  1:1487050/139‑1 (MQ=255)
                                                                                             gCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATTATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGc                                        <  2:1000897/139‑1 (MQ=255)
                                                                                              cATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGa                                              >  2:882960/1‑133 (MQ=255)
                                                                                              cATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGa                                              <  1:882960/133‑1 (MQ=255)
                                                                                               aTAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCaa                                      <  2:579049/139‑1 (MQ=255)
                                                                                                   ggAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGc                                        >  2:1135647/1‑133 (MQ=255)
                                                                                                   ggAACACGTTTTCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGc                                        <  1:1135647/133‑1 (MQ=255)
                                                                                                         cGTTTTCCGCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCa                            <  1:894052/139‑1 (MQ=255)
                                                                                                             ttCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGtt                                                      >  2:1385108/1‑109 (MQ=255)
                                                                                                             ttCCCCTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGtt                                                      <  1:1385108/109‑1 (MQ=255)
                                                                                                               ccccTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTg                         >  1:759779/1‑136 (MQ=255)
                                                                                                               ccccTGTGCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTg                         <  2:759779/136‑1 (MQ=255)
                                                                                                                    gtgCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt                                                                   >  2:968336/1‑89 (MQ=255)
                                                                                                                    gtgCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAAt                                                                   <  1:968336/89‑1 (MQ=255)
                                                                                                                    gtgCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATc                 >  1:217068/1‑139 (MQ=255)
                                                                                                                     tgCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGt                                                                       >  2:240241/1‑84 (MQ=255)
                                                                                                                     tgCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGt                                                                       <  1:240241/84‑1 (MQ=255)
                                                                                                                      gCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCtt               >  2:1075386/1‑139 (MQ=255)
                                                                                                                      gCACCCACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCtt               <  1:1342719/139‑1 (MQ=255)
                                                                                                                          ccACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGa                                                                                              <  2:1260015/56‑1 (MQ=255)
                                                                                                                          ccACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGa                                                                                              >  1:1260015/1‑56 (MQ=255)
                                                                                                                           cACGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagat          <  1:1075386/139‑1 (MQ=255)
                                                                                                                             cGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTgg                                                                                               <  2:1468781/52‑1 (MQ=255)
                                                                                                                             cGACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTgg                                                                                               >  1:1468781/1‑52 (MQ=255)
                                                                                                                              gACTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAAATTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCagataga       <  2:789941/139‑1 (MQ=255)
                                                                                                                               aCTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTg                                         <  1:1068086/104‑1 (MQ=255)
                                                                                                                               aCTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGc                                        <  1:1377802/105‑1 (MQ=255)
                                                                                                                               aCTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGc                                        >  2:1377802/1‑105 (MQ=255)
                                                                                                                               aCTAAATTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTGTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTg                                         >  2:1068086/1‑104 (MQ=255)
                                                                                                                                   aaTTCCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCa  <  1:1281238/139‑1 (MQ=255)
                                                                                                                                        |                                                                                                                                      
TTCTAAACGAACCAGTATAGAATTCGCCCTTTTGGCATGATATTAACTTGTCCAATATGATCAAATAGCATTAACCCCCCCTCACAACGTCCTGCATAGGGAACACGTTTTCCCCTGTGCACCCACGACTAAATTTCCCCCCCTTTAAACTTCTCAATAATGTCACCAATTTTCTGGAAATCCTATGATGCCTTTCCCGGTAAATGGTGCTGTGAGTTCAGAAAGATTCTGCAATCTATGTGCATTGAAAAAATCTTCAGATAGATATCA  >  NC_000913/2767365‑2767634

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